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ABSTRACT 


APPLICATION  OP  WELL-STIRRED  REACTOR  THEORY 
TO  THE  PREDICTION  OP  COMBUSTOR  PERFORMANCE 

By 

Allan  Hardy  Bonne 11 


Submitted  to  the  Department  of  Chemical  Engineering  on 
May  19.  1958  In  partial  fulfillment  of  the  requirements 
for  the  degree  of  Doctor  of  Science. 


Analysis  and  prediction  of  the  performance  character- 
istics of  complex  processes  has  often  been  successfully 
accomplished  by  use  of  single  conceptual  models.  In  the 
models  such  processes  are  visualized  in  terms  of  important 
basic  mechanisms  of  the  process  which  can  interact  in  an 
idealized  manner. 

Exemplary  in  the  field  of  combustion  has  been  the  use 
of  wrinkled  flame  front  models  to  analyze  the  behavior  of 
low-intensity  turbulent  flames.  In  high -output  combustion 
systems,  however,  substantial  portions  of  the  combustion 
\ chamber  may  be  devoted  to  homogeneous  combustion,  and  in 
\such  situations  wrinkled  flame  front  models  break  down. 

In  the  work  reported  here  conceptual  models  of  high- 
output  systems  are  considered  in  which  the  combustion 
chamber  is  visualized  as  consisting  only  of  a small  number 
of  interacting  elements  composed  of  homogeneous  reaction 
volumes.  Pressure  interaction  and  recirculation  are 
considered.  Performance  curves  developed  for  four  typical 
models  show  qualitatively  flow  interactions  and  residual 
flame  phenomena  typically  observed  in  real  combustion 
chambers . 

Predicted  performance  characteristics  based  on  the 
above  models  were  compared  with  experimental  data  obtained 
in  a small  combustion  chamber  burning  premixed  propane-air 
mixtures.  Stability  limits,  combustion  efficiency-loading 
rate  data,  and  visual  observation  of  flame  and  flow  were 
typical  measurements  made  on  chamber  performance. 


V 


The  performance  characteristics  of  combustion  chamoers 
predicted  by  these  models,  - including  the  existence  of 
primary  and  secondary  blowout,  the  effects  of  air-fuel  ratio 
on  blowout,  the  effects  of  recirculation  and  of  non-uniform 
fuel-air  ratio  on  combustion  efficiency,  the  existence  of 
conditions  which  can  lead  to  pulsing  combustion,  the 
magnitude  of  the  combustion  rate,  - all  these  are  in 
encouraging  agreement  with  measured  performance  of  actual 
chambers. 

Additional  experimental  work  is  required  to  establish 
more  quantitatively  the  ability  of  simple  conceptual  models 
to  predict  actual  combustor  performance,  but  the  experimental 
work  discussed  here  provides  reasonable  assurance  that  such 
models  have  practical  utility.  It  is  recommended  that  reasons 
for  the  deviations  between  the  model  predictions  and  the 
experimental  performance  be  investigated  further. 
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CHAPTER  I 

SUMMARY 

; j 

The  importance  of  mixing  as  a rate-limiting  step  in 
combustion  processes  has  long  been  recognized.  Prior  to 
the  development  of  high-output  combustors  for  Jet  propulsion, 
however,  mixing  rates  were  sufficiently  slow  in  most  practical 
combustors  to  make  chemical  reaction  rates  appear  infinitely 
j fast  by  comparison. 

The  use  of  a simple  conceptual  model  to  aid  in  visual- 
ization of  the  processes  involved  in  typical  high-output 
combustion  chambers  could  simplify  the  experimental  develop- 
ment of  new  chamber  designs.  Since  in  the  chambers  under  con- 
sideration both  chemical  reaction  rates  and  mixing  rates  have 

1 

I important  roles,  it  would  be  necessary  to  incorporate  both 

phenomena  into  such  a model.  The  model  should  predict 
behavior  characteristics  of  combustion  chambers  (such  as 
combustion  efficiency,  stability  and  tendency  to  resonate) 
resulting  not  only  from  limitations  set  by  chemical  reaction 
rate  or  mixing  rates,  but  also  from  interactions  of  two  or 
more  phenomena.  Although  a sufficiently  complex  model  could 
completely  describe  a chamber  behavior  pattern,  the  results 
of  application  of  the  model  concepts  to  prediction  of  per- 
formance would  be  Impressive  only  if  the  model  were  simple. 


In  the  models  discussed  here,  the  combustion  chamber 
is  considered  to  be  divided  into  a small  number  of  zones 
each  of  which  is  assumed  to  be  uniform  in  composition.  Flow 
between  these  zones  is  postulated  as  representing  the  fluid- 
flow  and  mixing  pattern.  Within  each  zone,  since  the  com- 
position is  uniform,  it  is  postulated  that  the  chemical  reaction 
rate  is  determined  uniquely  by  the  composition  and  temperature 
of  that  zone.  In  the  present  discussion  heat  transferors 
neglected;  however,  adiabaticity  is  not  a necessary  assumption. 
The  relation  between  burning  rate  and  composition  is  taken 
from  stirred  reactor  data  on  hydrocarbon  combustion  and  is 

V 

limited  to  a burning  rate-composition  relation  involving  only 
one  composition  variable.  More  complex  reaction  schemes  could 
be  used,  but  the  lengthy  computations  required  for  models 
involving  these  severely  limit  the  number  of  composition 
variables  which  may  be  considered.  Models  similar  in  some 
respects  to  those  discussed  here  have  appeared  in  the  recent 
literature  (2,  6,  18,  24);  however,  little  has  been  done  to 
relate  their  predicted  performance  quantitatively  to  the 
performance  of  practical  combustion  chambers. 

The  purposes  of  this  investigation  were:  (l)  to  consider 
some  possible  conceptual  models  of  combustion  chambers  which 
allow  for  interaction  between  flow  phenomena  and  reaction  rate 
phenomena,  (2)  to  study  the  operation  of  a simple  experi- 
mental combustion  chamber  and  (3)  to  examine  the  question 


of  whether  the  performance  predicted  by  simple  conceptual 
models  can  be  quantitatively  fitted  to  the  measured 
performance  of  the  experimental  combustion  chamber. 

Oata  with  which  the  predicted  performance  character- 
istics could  be  compared  were  obtained  from  a small  experi- 
mental combustion  chamber.  The  chamber  design  was  chosen 
to  be  simple  and  yet  to  give  characteristics  typical  of 
practical  chambers.  In  the  chamber  used  a small,  rectangular 
two-dimensional  combustion  2one  was  fired  with  either  one  or 
three  feed  jets.  Opposing  the  feed  jets  were  corresponding 

i 

exit  flues.  The  combustion  zone  was  formed  by  kaolin  fire- 
bricks whose  arrangement  could  be  readily  altered  depending 
on  the  number  of  feed  Jets  used.  The  front  and  back  sides 
of  the  combustion  zone  were  closed  with  Vycor  windows  which 
permitted  direct  observation  of  the  flame.  The  feed  Jets 
were  each  0.030"  wide  by  2"  deep  in  the  2-dimensional 
direction,  and  mounted  in  the  floor  of  the  chamber  with  their 
axes  accurately  aligned,  '■‘■'he  opposing  exit  flues  in  the 
chamber  roof  were  1/2"  wide.  When  three  Jets  were  used  the 
chamber  was  3"  wide  with  1"  spacing  between  jet  centerlines. 

With  one  feed  Jet  a 2"  chamber  width  was  used.  A pre-mlxed 
propane-air  mixture  was  prepared  and  metered  individually 
to  each  jet  to  insure  uniform  operation. 

Operating  characteristics  of  these  chambers  demonstrated 
two  general  phenomena.  In  using  the  three- Jet  chamber  it  was 

l . a 


found  that  the  three  Jet*  coalesced  a short  distance  above 
the  nozzles,  and  at  low  firing  rates  combustion  proceeded 
on  either  side  of  the  central  core  thus  created.  As  the 
firing  rate  was  increased  one  of  these  two  zones  became 
unstable.  Consequent  to  this  partial  blowout  there  was  a 
readjustment  of  the  flow  pattern  in  which  the  remaining 
active  zone  expanded  in  size.  At  a higher  firing  rate  it 
was  possible  to  blow  out  the  remaining  active  zone. 


In  the  one-jet  chamber  no  symmetrical  combustion  pattern 
was  observed.  The  Jet  from  the  centrally  located  feed  bent 
over  against  one  of  the  chamber  walls,  and  combustion  proceeded 
in  a circulatory  fashion  in  the  large  zone  thus  created,  at 
low  firing  rates  combustion  appeared  principally  around  the 
edges  of  the  circulatory  zone  with  the  "eye"  of  the  zone 
containing  only  completely  burned  gases.  As  the  flow  rate 
was  increased  the  burned  center  shrank  as  the  zone  of  most 
intense  combustion  moved  inward,  until  at  blowout  only  a 
small  residual  flame  in  the  center  of  the  vortex  remained. 

Pour  models  were  considered  in  this  work.  In  the  first 
of  these,  an  interaction  between  the  pressure  effect  of 
combustion  and  the  chemical  reaction  rate  was  considered. 

The  model,  purporting  to  simulate  a chamber  with  feed  and 
discharge  along  lines  located  centrally  on  opposite  walls, 
consisted  of  two  well-stirred  reactors  of  constant  total 
volume  separated  from  each  other  by  an  elastic  barrier 


capable  of  moving  in  response  to  pressure  differences  so  that 
the  deflection  was  directly  proportional  to  the  pressure 
difference  between  the  reactors.  The  two  reactors  were 
assumed  to  be  fed  equally  with  a pre-mixed  fuel-air  mixture 
and  when  the  pressure  difference  between  the  reactors  was 
zero,  the  two  reactors  had  equal  volumes. 

When  the  firing  rate  to  this  model  was  changed,  two 
stability  limits  were  predicted,  the  first  corresponding 
to  an  initial  blowout  of  one  of  the  two  reactors  followed 
by  a shift  in  the  barrier  position  toward  the  colder  side 
and  the  second  to  blowout  of  the  remaining  active  reactor 
at  a higher  feed  rate.  These  two  blowout  limits  corresponded 
to  loading  rates  in  the  zones  equivalent  to  stirred  reactor 
blowout  rates.  This  was  qualitatively  the  behavior  observed 
experimentally  in  the  three- Jet  chamber  studied.  Although 
a quantitative  fitting  of  the  predicted  performance  to  all 
of  the  data  on  the  experimental  chamber  was  not  possible, 
the  relative  firing  rates  between  the  two  blowout  limits 
was  predicted  with  some  success. 

The  second  combustor  model  considered  allowed  for  re- 
circulation of  partially  burned  combustion  products  in  a 
pattern  such  as  that  typically  produced  by  a Jet  firing 
fuel-air  mixture  into  a combustion  zone.  This  model  is 
similar  to  one  suggested  by  Spaulding  (24)  after  the  major 
portion  of  this  work  had  been  completed;  however,  the 


7 


performance  predicted  by  the  model  suggested  here  was 
drastically  different  from  that  predicted  by  Spaulding.  In 
the  recirculation  model  two  zones  are  postulated,  represent- 
ing the  central  Jet  core  and  recirculation  zone  found  in 
actual  combustors.  Each  of  these  zones  is  postulated  to  be 
a well-stirred  reactor.  Fresh  fuel-air  mixture  is  mixed 
instantaneously  with  partially  burned  recirculated  gases 
from  the  recirculation  zone  and  the  resulting  mixture  is 
then  fed  into  the  Initial  stirred  reactor  representing  the 
jet  core,  bases  leaving  the  central  core  volume  are  divided 
into  two  portions,  one  being  recirculated  through  the 
recirculation  zone  and  the  other  leaving  the  model  as  burned 
products.  This  model  showed  two  interesting  characteristics 
at  low  recirculation  ratios,  residual  flame  phenomena  and 
potential  flickering  flame  operation.  At  higher  recircula- 
tion ratios  these  phenomena  disappeared.  In  a third  model, 
similar  to  the  second,  allowance  was  made  for  non-uniform 
distribution  of  fuel.  This  was  achieved  by  introducing 
only  a part  of  the  fuel  but  all  of  the  air  into  the  central 
core  zone  and  the  remainder  of  the  fuel  into  the  recirculation 
zone  along  with  the  recirculated  gases.  The  third  model, 
with  the  split-fuel-feed,  showed  essentially  the  same 
behavior  as  that  of  the  second  model. 

Study  of  the  performance  of  the  experimental  one- Jet 
combustor  Indicated  that  a simple  recirculation  model  was 


inadequate  to  describe  its  performance.  To  fit  the  data  a 
fourth  model  was  considered  in  which  successive  recirculation 
zones  were  added.  In  this  model  a typical  zone  received 
recirculated  gases  from  its  immediate  successor  and  supplied 
recirculated  gases  to  its  immediate  predecessor.  This 
model  predicted  a performance  similar  to  that  of  the  simple 
recirculation  model,  but  the  relative  blowout  limits  for  the 
main  flame  and  the  residual  flame  were  in  better  agreement 
with  experimental  data  on  the  one- Jet  chamber. 

The  performance  characteristics  of  combustion  chambers 
predicted  by  these  models,  - including  the  existence  of 
primary  and  secondary  blowout,  the  effects  of  air-fuel  ratio 
on  blowout,  the  effects  of  recirculation  and  of  non-uniform 
fuel-air  ratio  on  combustion  efficiency,  the  existence  of 
conditions  which  can  lead  to  pulsing  combustion,  the 
magnitude  of  the  combustion  rate,  — all  these  are  in 
encouraging  agreement  with  measured  performance  of  actual 
chambers. 

Additional  experimental  work  is  required  to  establish 
more  quantitatively  the  ability  of  simple  conceptual  models 
to  predict  actual  combustor  performance,  but  the  experimental 
work  discussed  here  provides  reasonable  assurance  that  such 
models  have  practical  utility.  It  is  recommended  that  reasons 
for  the  deviations  between  the  model  predictions  and  the 
experimental  performance  be  investigated  further. 


CHAPTER  II 
INTRODUCTION 


Soope  of  Thesis 

Mixing  has  long  been  recognized  as  a dominant  factor 
In  many  Industrial  combustion  processes.  However,  there 
has  been  a growing  realization  that  chemical  reaction  rates 
place  important  limitations  on  the  performance  of  hlgh- 
output  combustors,  especially  those  found  In  aviation 
propulsion  systems.  This  nas  led  In  recent  years  to  tne 
study  of  the  well-stirred  reactor  as  a device  which  may 
be  used  to  relate  burning  rate  to  chemical  kinetics  as 
affected  by  composition,  temperature  and  pressure  rather 
than  to  mixing  processes  alone. 

Practical  combustors  for  hlgh-output  comDustion 


systems  lie  Intermediate  with  respect  to  volumetric  heat 
release  rates  between  purely  mixing-controlled  industrial 
furnaces  and  the  chemically  rate-limited  well-stirred 
reactor.  This  suggests  that  mixing  may  be  reduced  in 
importance  aB  a rate  limiting  process  in  many  practical 
combustors.  Furthermore,  much  experimental  data  obtained 
on  high-output  combustors  can  only  be  explained  by 
postulating  some  chemical  rate  limitations.  In  principle, 
therefore,  to  predict  the  behavior  of  hlgh-output  combustors 
from  fundamental  data,  one  would  have  to  account  for  the 
simultaneous  effects  of  both  chemical  reaction  rates  and 
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mixing  rates  in  the  combustion  chamber.  Both  of  these 
phenomena  are  highly  complex,  and  their  interactions  may 
lead  to  unexpected  results.  Typical  of  these  results  are 
the  marked  effects  of  different  firing  rates  on  combustion 
efficiencies,  and  discontinuous  shifts  in  flow  patterns 
and  efficiency  in  certain  firing  rate  ranges. 

One  might  hope  to  divide  the  various  phenomena  occur- 
ring in  a combustion  chamber,  particularly  of  efficiency 
and  stability,  into  regions  of  operation  where  either 
mixing  rate  or  chemical  reaction  rate  is  the  dominant 
phenomenon.  These  could  then  be  treated  as  independent 
problems  in  terms  of  the  present  theories  of  mixing  and 
of  chemical  kinetics.  This  method  has  the  deficiency, 
however,  of  being  unable  to  predict  performance  in  the 
transition  regions  where  many  interesting  problems  occur. 

To  study  these  regions  the  approach  proposed  here  is 
to  use  simplified  models  of  both  the  mixing  and  combustion 
rate  phenomena,  and  to  combine  them  to  allow  for  inter- 
action in  the  transition  regions.  The  simplified  conceptual 
models  of  combustion  chambers  cannot  be  expected  to  Include 
all  of  the  variables  or  effects;  however,  in  principle  a 
sufficiently  complex  model  should  be  able  to  completely 
describe  the  performance  of  any  combustion  chamber. 

In  setting  up  the  models  to  be  presented  here  knowledge 
of  simple  combustion  chamber  fluid  dynamics  and  of  well- 
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stirred  reactpr  theory  Is  required.  The  introductory  material 
is  thus  divided  into  three  sections.  The  first  section  deals 
with  the  fluid  dynamics  of  combustion  chambers,  the  second 
discusses  well-stirred  reactor  theory,  and  the  third  reviews 
the  literature  on  combustion  chamber  models. 

Combustion  Chamber  Fluid  Dynamics 

The  fluid-flow  pattern  of  a modem  can-type  combustor, 
typical  of  the  combustors  used  in  nearly  all  current  gas- 
turbines,  is  highly  complex,  and,  as  yet,  has  not  received 
adequate  attention.  A sectional  view  of  a typical  combustor 
is  shown  in  Figure  II-l.  Liquid  fuel  is  sprayed  into  the 
front  end  of  the  chamber,  and  combustion  air  enters  through 
a pattern  of  holes  around  the  surface  of  the  can. 


Figure  II-l: 


Typical  can-type  combustor 
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In  a device  such  as  this  moderate  mixing  Is  achieved  by  the 
"Jet"  effect  of  the  air  and  fuel  inlets.  This  mixing  must 
serve  three  purposes:  (l)  to  contact  the  fuel  and  oxidant 
so  that  combustion  can  occur,  (2)  to  provide  sufficient 
recirculation  of  hot  combustion  products  to  the  fuel  inlet 
so  that  the  flame  will  stabilize,  and  (3)  to  create  sufficient 
turbulence  to  burn  the  fuel  at  the  desired  volumetric  heat 
release  rate. 

Little  work  has  been  done  on  the  analysis  of  Jet  mixing 
in  chambers  where  there  is  mixing  of  recirculating  fluid 
with  the  Jet.  On  the  other  hand,  considerable  work  has  been 
done  on  Jet  mixing  in  pipes  where  the  recirculation  phenomenon 
has  been  eliminated  by  the  introduction  of  a secondary  fluid 
(11,  12,  29).  The  general  conclusion  is  that  so  long  as  the 
Jet  boundary  does  not  reach  the  wall,  the  behavior  of  the 
enclosed  Jet  is  identical  to  that  of  a free  Jet.  It  is 
reasonable  to  expect  that  the  same  results  would  be  found 
in  an  enclosed  system  with  recirculation. 

In  a practical  chamber  the  problem  is  considerably 
complicated  by  the  interaction  of  two  or  more  Jets  as  well 
as  by  deviations  of  the  Jet  from  idealized  theory  as  the 
confining  walls  are  reached.  In  a multi-jet  system,  individual 
Jets  tend  to  lose  their  identity  and,  instead,  make  their 
contribution  to  the  overall  turbulence  level  in  the  chamber. 

The  arrangement  of  the  Jets  may  be  such  as  to  give  flow 
patterns  ranging  from  no  identifiable  pattern  to  a well- 
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defined  pattern. 

In  the  combustion  chambers  considered  here,  the  flow 
pattern  Is  imagined  to  be  of  this  latter  type.  The  Jets 
are  arranged  to  give  a well -developed  circulatory  pattern, 
and  flame  is  stabilized  by  the  recirculation  of  hot  burned 
gases  to  the  fresh  feed  inlets.  Typically,  fresh  feed  and 
burned  products  counter-diffuse  (turbulently ) along  a 
portion  of  the  circulatory  pattern.  This  process  provides 
ignition  for  fresh  gases  in  the  mainstream  of  the  flow 
pattern  as  well  as  maintaining  active  combustion  in  the 

i 

sheltered  recirculation  zones. 

The  approach  to  well-stirredness  in  systems  of  this 
type  is  influenced  by  two  factors:  (l)  the  portion  of 
the  chamber  devoted  to  recirculation  and  (2)  the  level 
of  turbulence  in  the  chamber.  If  only  a small  portion  of 
the  chamber  is  devoted  to  gross  recirculation,  the  chamber 

[ 

will  approach  well-stirredness  only  through  extremely 
intense  turbulence  in  the  chamber.  Likewise,  if  the 

'3 

gross  recirculation  pattern  is  of  a large  diameter,  high 
turbulence  is  required  to  couple  the  central  portion  of  the 
pattern  with  the  outer  regions.  Spaulding  (24)  suggests 
that  gross  circulation  is  a more  efficient  means  of  achieving 
mixedness  than  is  high  turbulence.  (Efficiency  is  in  terms 
of  pressure  drop  required  to  achieve  a given  level  of  mixedness.) 
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Because  of  high  reaction  rates  of  combustion,  incomplete 
mixing  within  a chamber  leads  to  zones  of  high  and  low  com- 
bustion rates.  Consequently,  completeness  of  mixing  has  been 
associated  with  high  volumetric  heat  release  rates.  It  has 
been  postulated  that  it  is  possible  to  obtain  even  higher 
volumetric  combustion  rates  by  proper  arrangement  of  combustor 
elements  (6);  however,  these  have  not  yet  been  achieved  in 
practice.  In  visualizing  models  of  combustion  chambers  such 
as  these,  the  effects  of  recirculation  which  can  be  important 
in  some  regions  must  be  recognized  in  the  model. 

Combustion  Beactlon  Kinetics:  tfell-Stlrred  Reactor  Behavior 

The  subject  of  combustion  reaction  kinetics  has  been 
treated  by  a large  number  of  authors.  General  reviews  of 
high  temperature  kinetics  have  been  given  by  Shuler  (21), 
Ubbelholde  (27),  and  vonElbe  (28).  This  subject  will  not  be 
treated  in  detail  here;  only  those  aspects  of  the  problem 
directly  applicable  to  the  present  work  will  be  discussed. 

By  far  the  largest  proportion  of  chemical  kinetics 
studies  have  been  made  in  the  region  of  low  or  moderate 
temperatures  compared  to  the  temperatures  of  normal  combustion 
reaction  (for  example,  extensive  work  on  cool  flames).  These 
studies  describe  the  combustion  reactions  in  highly  complex 
terms,  and  because  of  the  relatively  low  temperatures  Involved, 
there  is  considerable  doubt  that  these  data  can  be  extrapolated 
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to  normal  combustion  temperatures.  Furthermore,  because  of 
the  primary  Interest  here  In  the  Interaction  of  combustion 
and  fluid  flow,  It  Is  desirable  to  employ  a much  simpler 
picture  of  the  chemical  kinetics. 

Such  a simple  picture  of  combustion  kinetics  is  one  in 
which  it  is  hoped  that  a simple  nth  order  reaction  between 
fuel  and  oxygen  will  adequately  describe  the  combustion 
kinetics.  This  assumption  has  been  used  by  Longwell  (13) 
and  Avery  (2)  in  the  analysis  of  a well-stirred  reactor. 

It  is  assumed  that  the  reaction  order  is  of  a order  in 
fuel  and  n-a  order  in  oxygen  and  that  Arhennius  form  of 
activation  energy  represents  the  effect  of  temperature  on 
the  reaction  rate. 

Let  the  overall  reaction  between  fuel  and  oxygen  be: 

F + a 03  — > p P II-l 

and  let  0 times  the  stolchiometrically  required  fuel  be  used. 

If  the  identifiable  components  of  the  mixture  at  any 
stage  of  completion  of  combustion  are  products  of  complete 
combustion,  fuel,  and  oxygen,  specification  of  a single  quan- 
tity, for  example  0,  the  burned  fraction  of  the  fuel  fed, 
together  with  the  stoichiometry,  is  sufficient  to  identify 
the  mixture.  For  each  additional  intermediate  or  competing 
reaction  involving  a new  specie  (or  species),  one  more 
quantity,  7,  6 , etc,  must  be  specified  to  fix  the  composition; 
and  the  choice  of  how  best  to  define  the  new  term  will  depend 


on  the  postulated  chemistry  of  the  process.  The  mole  fraction, 
f x»  of  any  component  "X"  is  thus  representable  as  f^ 

(0>  P,  y,  8 i ...)»  or  for  a fixed  value  of  0,  f^  (0,  y, 
fe  , ...),  and  the  single  argument  0 is  sufficient  for  the 
simplest  possible  picture  of  the  combustion  process. 

The  well-stirred  reactor  is  by  definition  a chamber  in 
which  the  composition,  measured  at  the  inlet  by  0Q,  yQ,  gQ, 
changes  instantly  to  that  measured  by  0,  y,  $ ; and  the 
composition  is  everywhere  that  of  the  outlet.  In  this  system 
the  temperature,  T,  is  fixed  by  the  exit  composition  and  the 
inlet  temperature  Tq  (i.e.  an  adiabatic  reaction),  and  for  a 
fixed  0 can  be  represented  as  fT  (0,  y,  g , ...,  TQ) . For 
the  simplest  case  in  which  the  product  composition  and, 
consequently,  T is  assumed  to  be  fixed  by  the  single  variable 
0,  it  is  obvious  that  either  0 or  T may  be  used  as  a measure 
of  combustion  progress  or  burnedness.  Assumption  of  linearity 
between  T and  0,  leading  to  the  relation 


(where  I ^ is  the  adiabatic  flame  temperature)  simplifies 
numerical  calculation,  but  is  not  necessary. 

As  many  chemical  rate  equations  are  required  as  there 
are  arguments  0,  y,  £ , ...  necessary  to  fix  the  composition. 


Additional  rate  equations  will  be  required  if  alternative 
routes  between  reactants  and  products  exist.  The  present 

Lu.  . i . 
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treatment  is  restricted  to  the  simplest  description  of  the 
combustion  process,  requiring  only  one  rate  equation.  This 
is  most  conveniently  the  rate  of  consumption  of  fuel,  given 
by; 


- & - Cpa 


II-2 


on  the  assumption  that  the  reaction  is  of  overall  nth  order 
and  of  a order  in  fuel.  Substituting  fp(0)(P/B3)  for  Cp 
and  fQa  (0)  (P/RT)  for  C^,  one  has; 

dNs>  n _ a _ ( -a) 

“ Vdt  “ **  gnJ(n-l/2)  fP  (P)  f°2  (P)  11-3 

Since  T and  0 are  uniquely  related,  Equation  II-3  may  be 
written: 


dN, 


P^V  dt 


f i (P> 


11- 3a 


In  a well-stirred  reactor  fed  at  the  molal  fuel  feed 
rate  G with  the  burnedness  changing  from  an  inlet  value  0Q 
to  an  outlet  value  0,  the  fuel  consumption  rate  is  g(0-0q); 
this  replaces  - dN^/dt.  Therefore,  Equation  II-3a  becomes; 

11-4 


G _ (0) 

VP*  " £ - B 


and  for  a reactor  being  fed  an  unburned  fuel-air  mixture, 
this  may  be  written: 

G _ Vx  (0) 

VP*1  0 


$(0> 


II-5 


- 18  - 


A typical  curve  representing  Equation  II-5  appears  as  Fig.II-2. 
The  numerical  magnitude  of  G/VP11  shown  in  Fig.  II-2  is  a 
consequence  of  the  normalization  proceeding  and  is  not 
indicative  of  the  actual  values  experimentally  observed. 

The  normalization  of  this  curve  is  discussed  in  Appendix  A. 

It  will  be  noted  that  there  are  three  solutions  to 
Equation  II-5  for  a given  value  of  G/VP31  at  low  values  of 
this  quantity.  Of  these,  only  two  may  be  experimentally 
found.  The  uppermost  (solid  line)  solution  represents  the 
hot  , stably -burning  system,  and  the  lowermost  (po  axis) 
the  cold,  "non-reacting"  system.  Solutions  along  the  dotted 
middle  branch  of  Figure  II-2  are  thermally  unstable,  since 
if  a reactor  were  operating  on  this  section  of  the  curve,  a 
small  perturbation  in  burnedness  or  flow  rate  would  cause 
the  reactor  efflux  either  to  move  up  to  the  top-most  branch 
of  the  curve  or  to  drop  down  to  the  lowest  branch  of  the 
curve.  As  the  feed  rate  is  increased,  the  burnedness 
decreases,  until  a feed  rate  is  reached  where  dp/dG  is 
infinite.  This  point  represents  the  blowout  feed  rate,  since 
beyond  this  feed  rate  the  only  solution  to  Equation  11-5 
corresponds  to  the  "non-reacting"  system.  This  blowout  point 
is  a unique  function  of  the  fuel/air  ratio,  and  has  been  a 
basic  measurement  in  the  work  which  has  been  reported  in 
the  literature  on  the  well-stirred  reactor. 

In  applying  the  well-stirred  reactor  theory  to  experi- 
mental data,  Longwell  and  Weiss  (14)  and  Baker  (3)  have  used 
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FIGURE  H-2:  WELL-STIRRED  REACTOR  PERFORMANCE 
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the  simple  (one  rate  equation)  theory.  In  this  approach  no 
allowance  is  made  for  the  existence  of  intermediate  combustion 
products.  Hence,  any  products  of  partial  combustion  are 
replaced  by  stoichlometrically  equivalent  quantities  of 
unburned  fuel  and  completely  burned  products  such  that  the 
oxygen  consumption  efficiency  remains  unchanged.  The 
burnedness,  p,  is  defined  in  terms  of  this  latter  quantity 
in  these  experiments. 

Since  it  is  admittedly  unrealistic  to  assume  that  there 
are  no  Intermediate  combustion  products  in  the  well-stirred 
combustor,  both  Longwell  and  Veiss  and  Baker  have  spent 
considerable  effort  in  postulating  stoichiometric  and  re- 
action schemes  which  will  give  a more  realistic  picture  of 
the  combustion  products  actually  found.  In  the  experimental 
work  reported  by  Longwell  and  Weiss  (14)  gas  analyses  of 
samples  drawn  from  the  interior  of  the  reactor  were  reported. 
These  data  show  that  very  little  unburned  fuel  appears  in 
the  reactor,  and  that  only  when  the  combustor  is  fed  a rich 
fuel-air  mixture  do  any  appreciable  quantities  of  hydrocarbons 
appear.  In  these  samples  the  major  unburned  component 
reported  is  carbon  monoxide,  and  hydrogen  is  the  only  other 
combustible  reported  in  significant  quantities  for  lean 
mixtures.  Baker  (3)  in  sampling  data  taken  in  another  well- 
stirred  reactor  experiment  confirms  these  observations. 
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Longwell  suggested  the  following  reaction  equations  to 
explain  the  data: 

Fuel  (hydrocarbon  C H ) + aOa-*-xCO  + l/2yH30 (instantaneous) 

x y 

II-6 


CO  + 1/2  0a— *-C0a  (rate  controlling) 


11-7 


This  scheme  allowed  Longwell  to  fit  his  fuel-lean  data; 
however,  it  is  not  adaptable  to  fuel-rich  data  since  there  is 
insufficient  oxygen  to  allow  completion  of  the  first  step 
beyond  a certain  fuel-air  ratio  (depending  on  the 
carbon/hydrogen  ratio  of  the  fuel) . 

In  fitting  this  reaction  scheme  to  the  simple  well- 
stirred  reactor  equations  (Equation  II-5) , one  must  be 
careful  not  to  confuse  the  fractional  oxygen  consumption 
with  the  fractional  carbon  monoxide  consumption.  Since 
the  rate -control ling  reaction  is  equation  H-7,  the 
substitutions  into  II-5  must  be  appropriate  to  this 
stoichiometric  equation.  The  fractional  fuel  consumption, 0, 
is  now  defined  as  the  fraction  of  CO  available  from  the  fuel 
which  has  been  converted  to  C0a.  Temperatures,  accounting 
for  the  heat  liberated  by  Equation  II-6,  and  concentrations 
are  calculated  as  a function  of  this  redefined  0.  It  should 
be  noted  that  this  definition  of  burnedness  is  different 
from  oxygen  consumption  efficiency,  since  at  zero  CO 
consumption  the  oxygen  consumption  efficiency  is  that 
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required  to  react  all  of  the  fuel  to  CO  and  Ha0.  Baker  polnta 
out  that  proper  substitution  into  Equation  II-5  leads  to  the 
conclusion  that  for  near-stoichiometric  mixtures  on  the  lean 
side,  enough  heat  is  liberated  by  reaction  II-6  to  prevent 
blowout  of  the  reactor,  no  matter  how  fast  it  is  fed. 

Longwell's  correlation  of  his  lean  blowout  data,  unfortunately, 
confuses  this  distinction.  Although  his  theoretical  blowout 
curve  extends  to  stoichiometric  fuel-air  ratio,  close 
examination  of  this  derivation  reveals  that  an  improper 
substitution  was  made  for  p in  the  denominator  of  Equation 
II-5. 

Baker  suggested  two  additional  reaction  schemes,  one 
for  the  rich  side  and  one  for  the  lean  side.  Because  of  the 
anomolous  behavior  of  the  reaction  scheme  of  Equations  II-6 
and  II-7  (i.e.  no  blowout  found  near  stoichiometric),  Baker 
suggested  an  alternative  scheme  in  which  there  was  less 
heat  liberated  in  the  first  stage  of  the  reaction.  This 
scheme  is  as  follows: 

C H + aOa  — xCO  + 1/2  yHa  (instantaneous)  11-8 

X,  J 

CO  + 1/2  0a  — *■  C0a 

(rate  controlling)  II-9 

Ha  + 1/2  0a  — ♦■Ha0 


and  the  restriction  is  placed  on  11-9  that: 
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where  k la  the  water-gas  equilibrium  constant  at 
the  temperature  of  the  reactor. 

On  the  rlgh  aide  it  la  postulated  that  fuel  and  oxygen 
are  consumed  in  the  same  ratio  as  the  prevailing  fuel/alr 
ratio  to  form  the  equilibrium  mixture  of  carbon  monoxide, 
carbon  dioxide,  hydrogen  and  water.  This  scheme  is  the 
same  as  that  proposed  by  Longwell  and  Weiss  for  rich 
mixtures. 

In  Baker's  lean-side  reaction  scheme,  II-8  to  11-10, 
there  is  a much  lower  temperature  rise  in  the  first  step  than 
in  the  scheme  proposed  by  Longwell  and  Weiss,  and  Baker  was 
aole  to  calculate  blowout  rates  from  this  mechanism  for  all 
lean  fuel/alr  mixtures.  However,  it  was  pointed  out  that 
the  ratios  of  CO,  C0a,  Ha  and  Ha0  experimentally  found  in 
the  gas  samples  did  not  correspond  to  water-gas  equilibrium. 

Edgerton,  Saunders,  and  Spaulding  (8)  have  discussed 
the  well-stirred  reactor  theory  in  general,  and  they  have 
pointed  out  some  of  the  implications  of  the  simplified 
theory  dealt  with  thus  far.  The  composition  of  a partially 
burned  mixture  has  been  defined  by  a single  quantity.  This 
definition  was  originally  restricted  to  that  case  where  the 
rate-limiting  step  in  the  combustion  of  fuel  was  the  initial 
reaction  between  fuel  and  oxygen.  However,  in  any  reaction 
scheme  in  which  a single  reaction  is  the  rate-limiting 
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step,  the  composition  of  the  mixture  may  be  uniquely  defined 
by  a single  argument  such  as  0.  Thus,  a partially  burned 
mixture  whose  composition  Is  measured  by  0 may  be  obtained 
either  by  partial  combustion  to  this  composition  or  by 
mixing  unburned  fuel/air  mixture  with  completely  burned 
gases.  In  this  latter  case,  all  Instantaneous  reactions 
must  be  allowed  to  occur  (e.g.  the  conversion  of  hydrocarbon 
to  CO  and  Ha0,  Equation  11-6)  without  any  conversion  of  the 
reactants  involved  in  the  rate-limiting  step.  This  postu- 
lation that  the  composition  of  a partially  burned  mixture  is 
uniquely  defined  by  a single  parameter  Is  of  fundamental 
Importance  in  calculating  the  behavior  of  combustion 
chamber  models. 

DeZubay  (6)  has  called  attention  to  an  ambiguity  In 
the  meaning  of  "burnedness'1  which  arises  from  Baker's 
and  Longwell's  assumption  of  stoichiometry  for  rich  fuel-air 
mixtures  (e.g.  fuel  and  oxygen  disappear  In  the  same  pro- 
portion as  the  prevailing  fuel/air  ratio).  For  this  stoichi- 
ometry a burnedness  of  unity  for  a rich  fuel-air  mixture 
corresponds  to  a mixture  In  which  there  are  no  unreacted 
hydrocarbons,  in  conformity  with  the  experimental  facts. 

It  Is  evident,  however,  that  the  result  of  mixing  a 
completely  burned  rich  and  lean  mixtures  would  be  a mixture 
In  which  further  combustion  could  occur.  beZubay  has 
suggested  that  this  ambiguity  in  the  meaning  of  burnedness 
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can  be  clarified  by  the  use  of  a "fuel  availability  factor" 
in  addition  to  "bumedness"  to  specify  the  composition  of 
the  mixture.  This  composition  variable  permits  one  to 
account  for  the  additional  combustion  which  could  occur  if 
air  or  a lean  fuel-air  mixture  is  mixed  with  a burned  or 
partially  burned  rich  mixture. 

The  introduction  of  multi-step  reactions  with  more 
than  one  rate  constant  and  composition  variable  considerably 
complicates  the  stirred  reactor  problem.  This  possibility 
has  not  yet  received  adequate  attention  in  the  literature. 

Edgerton,  Saunders  and  Spaulding  (8)  In  their  paper 
discuss  this  problem  qualitatively.  They  point  out  that  so 
long  as  one  of  the  reaction  rates  is  substantially  slower 
than  the  others  the  simple  formulation  of  the  stirred- 
reactor  problem  is  still  valid.  They  further  suggest  that 
near  the  completion  of  the  combustion  process,  a single 
reaction  predominates,  thus  making  valid  the  assumption  that 
composition  is  independent  of  the  path  by  which  the  burned- 
ness is  reached.  In  the  intermediate  ranges  of  burnedness, 
however,  it  is  postulated  that  the  assumption  of  only  one 
significant  rate  constant  is  unjustified.  If  more  than  one 
rate  constant  is  required  to  describe  the  combustion  reaction, 
then,  in  the  general  case,  the  composition  and  temperature 
are  not  uniquely  defined  by  a single  variable,  but  they 
depend  on  the  path  by  which  a given  burnedness  level  is 

l—*- _ - _ Jl 
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reached.  To  allow  for  variation  in  the  burning  rate  at  a 
given  temperature,  Edgerton,  Saunders  and  Spaulding  suggest 
a "band  of  uncertainty",  which  includes  all  possible  burning 
rates  corresponding  to  a given  temperature  level. 


Behavior  of  Plug-flow  and  Staged  Stirred  Beactors 

Because  the  performance  of  two-stage  stirred  reactors, 
multi-stage  stirred  reactors,  and  plug-flow  reactors  are 
all  derivable  from  the  same  basic  rate  equation,  (Equation 
II-5) , they  are  all  derivable  from  one  another.  Consider 
first  the  plug-flow  reactor.  Such  a reactor  may  be 
pictured  as  an  Infinite  series  of  Infinitesimal  well-stirred 
reactors.  Since  the  flow  rate  through  each  of  these  is  the 
same,  they  are  best  considered  as  stirred  reactors  in  which 
volume  is  the  Independent  variable.  Equation  11-5, 
restricted  to  the  case  where  only  the  single  argument  0 is 
required  and  for  constant  pressure  operation,  is  inverted 
to  give: 


VP»  ^ 0 - 0p 

G if/  i (0) 


11-11 


A typical  slab  of  a plug-flow  reactor,  having  a volume  of 
dV,  being  fed  with  reactants  of  composition  0,  and  dis- 
charging products  of  composition  0 + d0,  is  represented  by 
Equation  11-11  with  V replaced  by  dV  and  the  accomplished 
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burnedness,  p - po  , by  dp  ; 


£ 

0 


dV 


11-12 


If  the  reaction  rate  function,  /'i,  Is  available,  as  for 
example  from  a postulated  chemical  mechanism  of  combustion, 
Equation  11-12  may  be  Integrated  directly  to  give  the 
performance  of  a plug-flow  reactor.  If,  on  the  other  hand, 
data  are  available  from  a stirred  reactor  experiment  In  the 
form  of  (fj  8 (p)  ( ipz  - (ft i/P  ),  then  the  stirred 

reactor  relationship  may  be  used  to  integrate  11-12: 

Equation  11-12  may  be  written: 

£av  = fyrm  dlne  H-13 


or,  integrating  between  the  limits  papo, 
P=P,  V=V  : 


V«VQ  and 


11-14 


If,  one  plots  the  stirred  reactor  performance  curve  In  the 
form  1/  jtx  (p)  vs  p on  a logarithmic  scale,  the  area 
under  the  curve  from  PQ  to  p (see  Pig.  II-3  with  x and  y 
coordinates  reversed  to  correspond  to  previous  plots) 
represents  the  generalized  volume  (V-Vo)Pn/G  necessary  to 
carry  out  that  much  reaction  In  plug-flow;  and  that  area 
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FIGURE  H-3: 


INTEGRATION  OF  WELL-STIRRED  REACTOR  CURVE  TO  FIND 
PLUG-FLOW  REACTOR  PEFORMANCE 


can  be  plotted  on  the  same  plot  (dotted  line) . In 
Figure  II -4  the  results  of  Integration  over  the  full  range 
of  0 Is  shown  on  arithmetic  coordinates  (dotted  line) 
together  with  the  basic  stirred  reactor  curve.  The  plug- 
flow  curve  has  been  displaced  laterally  to  be  tangent  to 
the  stlrred-reactor  curve  at  a generalized  volume  of 

the  significance  of  which  requires  comment:  On  a 
plot  of  p vs  VP^/G  on  arithmetic  coordinates,  the  point 
of  tangency  to  a line  from  the  origin  represents  a stirred 
reactor  operating  at  the  maximum  possible  rate  per  unit 
volume,  and  the  condition  for  this  point  is: 

£ « ££ 11-15 

1/  tp  3(p)  d [1/  Ip  2(p)] 

Differentiating  Equation  II-5  and  11-14  with  respect  to 
volume  (these  are  the  curves  represented  on  Figure  II-4), 
one  obtains: 

^ dV  a da/  #2  fg)3  P^dV  _ 1 fllBfi. 

G dp  dp  G dp  <lJ 3 (p)  dp 

VSR  PFR 

Noting  that  the  point  where  these  two  slopes  are  equal  is 
given  by: 

d (1/ if/  a (p) ) * d In  p / (p9  (p)  II-16 

one  sees  that  tangency  occurs  at  the  p,  V combination 
corresponding  to  the  maximum  burning  rate  In  a stirred 
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reactor.  Note  that  for  any  value  of  p above  p,  the 
Incremental  volume  In  a plug-flow  reactor  Is  less  than 
in  a well-stirred  reactor,  but  to  increase  the  burnedness 
from  any  point  lying  between  p2,  the  stability  limit  of 
the  stirred  reactor,  and  Px  to  the  value  Pi,  a greater 
volume  is  required  in  a plug  flow  reactor  than  in  a well- 
stirred  reactor.  For  points  below  p2  the  well-stirred 
reactor  is  unstable.  The  minimum  total  volume  required  to 
accomplish  any  burnedness  lying  above  Px  is,  therefore, 
the  sum  of  the  volume  of  a well-stirred  reactor  operating 
to  give  products  at  Px  followed  by  a plug-flow  reactor  to 
reach  the  burnedness  p. 

It  is  noticed  that  the  plug-flow  curve  of  Fig.  II-4 
looks  very  much  like  the  profile  of  a laminar  flame  front 
in  which  bumedness  has  been  plotted  against  L/S^,  the 

measure  of  time.  However,  in  the  calculation  for  Figure  II-3 
there  is  no  allowance  for  diffusion  of  heat  or  mass  such  as 
gives  rise  to  a laminar  flame  propagation  rate,  and  although 
the  abscissa  has  dimensions  proportional  to  a length/flame 
speed  ratio,  there  is  no  definite  distance  associated  with 
the  abscissa;  consequently,  there  is  no  clear  analogue  to 
flame  propagation. 

The  performance  of  a well-stirred  reactor  being  fed  a 


partly  burned  feed  is  readily  obtained  from  the  primary 
stirred  reactor  curve  such  as  Figure  II-5,  in  which  p is 
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plotted  against  1/  [ft  a (equal  to  VI^/G) , The  slope  of  a 
line  from  the  origin  to  a point  on  the  curve  represents  the 
burning  rate  per  unit  volume  In  any  system  In  which  the  bur- 
nedness is  0.  Since  when  the  feed  Is  partially  burned,  the 

burnedness  to  be  accomplished  is  0 - 0Q  rather  than  0,  the 
generalized  volume  necessary  Is  (VP^/G)  (0-0o)/(0). 

This  appears  In  Figure  II-5  as  MVP^/G).  Figure  II-5  can 
also  be  considered  an  example  of  operating  two  reactors  In 
series.  The  first  reactor,  fed  with  the  unburned  fuel 
(0Q  = 0),  requires  the  volume  V^P/'G  to  deliver  products 
of  composition  0i,  and  the  second  reactor,  fed  with  products 
of  burnedness  0i  and  delivering  products  of  burnedness  0a » 
requires  the  additional  volume  A (VF^/G) . Application  of 
this  principle  leads  directly  to  the  following  conclusions 
for  a pair  of  stirred  reactors:  (Note  the  inset  in 
Figure  II-5.  Two  lines  to  the  origin  have  been  drawn,  one 
tangent  to  the  nose  of  the  curve,  and  the  other  passing 
through  the  blowout  point  of  the  curve.)  For  cases  where 
the  burnedness  leaving  the  second  reactor  lies  above 
point  "C",  and  all  combinations  of  reactor  pairs  will 
require  less  volume  than  a single  stirred  reactor  giving 
the  same  final  burnedness.  If  the  burnedness  leaving 
the  second  reactor  lies  between  points  "B"  and  "C"  it  is 
always  possible  to  have  a pair  of  reactors  operating  more 
favorably  than  a single  reactor,  although  not  all  reactor 
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pair*  In  this  range  are  better.  If  the  final  burnedness 
leaving  the  second  reactor  lies  between  "A"  and  "B"  there 
Is  no  possible  combination  of  reactor  pairs  which  requires 
less  total  volume  than  a single  reactor. 


The  discussion  of  reactor  combination  has  been  limited, 
thus  far  to  the  case  where  only  one  composition  variable  Is 
required  to  describe  the  mixture.  If  two  or  more  variables 
are  required,  the  treatment  of  staged  reactors  is  considerably 
complicated.  The  problems  involved  are  best  illustrated  by 
consideration  of  a simple  case. 

Consider  a combustion  process  in  which  fuel  is  partially 
oxidized  in  the  first  step  to  form  "1"  moles  of  an  Inter- 
mediate "1"  and  "I"  then  reacts  with  additional  oxygen  to 
form  complete  combustion  products.  For  this  reaction  two 


rate  equations  are  required. 

F + a03  *• 


II  + a30s 


ll  Step  1 
pP  Step  2 


ii-i? 

11-18 


11-19 
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In  a manner  similar  to  the  previously  discussed  equation 
for  the  simple  stirred  reactor  kinetics,  these  equations 
may  be  reduced  to  functions  of  two  variables  0 and  y, 
representing  the  progress  of  reaction  of  steps  1 and  2 
respectively: 

^ (P  - PG>  * H-21 

— - (nr  - y)  = fa (0,7)  n-22 

vph  O 

Simultaneous  solution  of  these  equations  is  required  to  find 
the  performance  of  a one-stage  stirred  reactor  at  any  given 
feed  rate.  If  it  is  desired  to  operate  a two-stage  reactor, 
the  gas  composition  of  the  final  stage,  in  general,  will 
not  be  the  same  as  the  composition  of  a one-stage  reactor 
whose  exit  temperature  (or  fractional  oxygen  consumption) 
is  the  same;  and  a second  simultaneous  solution  is  required 
for  the  0Q,  70  combination  corresponding  to  the  partially 
burned  gases  leaving  the  first  stage. 

To  simplify  calculations  of  these  equations  it  is 
desirable  to  represent  graphically  the  simultaneous 
solutions  obtained  from  any  given  0Q,  *Y0  starting  combina- 
tion (obtained  in  an  unspecified  manner) . The  composition 
of  the  exit  gases  from  any  reactor  may  be  represented  on 
a diagram  of  0 vs  y for  this  two  reaction  postulated 
mechanism.  On  this  same  diagram  lines  of  0 y combinations 


35 


having  a constant  ratio  of  the  relative  disappearance  rates 
of  fuel  and  intermediate  may  also  be  mapped  (d0/d7  = C) 

(Figure  II-6  ).  Thus,  in  a stirred  reactor  whose  outlet 
composition  is  represented  by  0 7 , the  dp/ d 7 in  this  reactor 

may  be  read  directly  from  this  diagram,  and  this  ratio  must 
be  equal  to  the  ratio  of  (0  -P0)/(Y  -7g)  in  that  reactor. 
Therefore,  for  any  combination  representing  either  outlet 
or  inlet  conditions  of  a stirred  reactor,  lines  representing 
possible  inlet  or  outlet  0,7  combinations  may  be  plotted. 

Having  found  the  inlet  and  outlet  conditions  from  a diagram 
such  as  Figure  II-6,  recourse  to  either  Equation  11-21  or 
11-22  will  give  the  volume  required  in  that  reactor.  These 
equations  may  be  represented  graphically  to  speed  calcula- 
tions (Figure  II-7). 

Using  these  two  plots  the  performance  of  any  desired 
combination  of  plug-flow  and  stirred  reactors  in  which  a 
two-step  kinetic  scheme  is  postulated  may  be  calculated. 

It  is  evident,  however,  that  the  addition  of  a third 
variable  would  complicate  the  problem  still  further,  because 
of  the  unwieldy  calculation  required  for  rigorously  treating 
the  combination  of  stirred  reactors  with  multi-step 
reactions,  one  must  hope  that  a simple  one-step  reaction 
scheme  will  provide  an  accurate  enough  representation  of  the 
gross  kinetics  to  permit  calculation  of  model  performances. 

This  will  be  true,  if,  as  suggested  by  Edgerton,  Saunders 
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and  Spaulding,  one  reaction  predominates  over  all  others 
in  the  range  of  practical  interest. 

Experimental  Stirred  Seactor  Performance 

Experimentally,  the  study  of  the  stirred  reactor  in 
combustion  problems  presents  primarily  the  problem  of 
mixing  fresh  feed  gases  with  the  gases  in  the  reactor  at 
an  extremely  high  rate.  This  high  mixing  rate  is  required 
to  mix  the  gases  before  a significant  amount  of  reaction 
can  occur.  This  obviously  cannot  be  accomplished  by  simple 
mechanical  stirring.  After  some  developmental  work, 

Longwell  and  Weiss  (14)  proposed  the  now  familiar  stirred 
reactor  in  which  a combination  of  Intense  turbulence  and 
gross  recirculation,  caused  by  a large  number  of  small 
Jets  issuing  from  a small  concentric  sphere  into  the 
spherical  reactor,  mix  the  reacting  gases  moderately 
well. 

The  data  which  have  been  published  on  this  device 
have  been  primarily  blowout  data.  Three  investigators 
(3,  5,  14)  have  studied  blowout  loading  factors  as  a 
function  of  fuel-air  ratio  in  reactors  of  this 
design.  Longwell  and  Weiss  (14)  were  the  first  to 
study  this.  They  reported  a pressure  exponent  of 
1.8  as  giving  the  best  correlation 
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of  their  data,  and  they  claimed  that  there  was  no  significant 
effect  of  hole  size.  Baker,  In  a somewhat  later  Investiga- 
tion, reported  a trend  of  pressure  exponent  with  fuel-air 
ratio  and  recommended  a pressure  exponent  of  1.3.  He  also 
found  a consistent  effect  of  hole  size.  Because  of  the 
differences  In  the  results  reported  by  Longwell  and  Weiss 
and  by  Baker,  Blichner  (5)  conducted  a further  Investigation 
on  Baker's  reactor.  Blichner' s results  substantiated  those 
of  Baker.  He  found  an  exponent  lower  than  that  reported  by 
Longwell  and  Weiss,  and  showed  a significant  effect  of  hole 
size.  A comparison  of  the  blowout  data  reported  in  these 
three  Investigations  Is  shown  In  Figure  11-8.  Throughout 
this  thesis  data  of  this  type  are  reported  on  coordinates 
of  air  loading  rate,  NA/V  P^’^and  generalized  fuel  fraction 
F.  This  latter  quantity  is  given  by  0/1  + 0 , where  0 is 
the  fraction  of  stoichiometric  fuel.  F Is  0.5  for  a 
stoichiometric  mixture,  lean  mixtures  are  less  than  F = 0.5, 
and  rich  mixtures  greater.  In  this  comparison  all  three 
sets  of  data  are  shown  for  a pressure  exponent  of  1.8,  and 
no  distinction  between  hole  size  Is  made.  All  three  sets 
of  data  are  corrected  to  an  inlet  temperature  of  300  °K  using 
the  temperature  correction  factors  suggested  by  Longwell 
and  Weiss. 


Both  Longwell  and  Weiss  and  Baker  have  reported 
sampling  data.  Sampling  traverses  were  made  along  a 
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radius  of  the  reactor  to  determine  to  what  extent  composition 
varied  within  a reactor.  Longwell  and  Weiss  reported  that 
the  composition  was  substantially  constant  over  about  95^ 
of  the  volume.  Baker  points  out,  however  that  for  their 
sampling  runs  the  burnedness  level  in  the  reactor  was  high 
(over  90%),  and  consequently  even  if  there  were  a composition 
change  it  would  not  be  great.  Liongwell  and  Weiss'  data  if 
plotted  on  percent  unburned  basis  show  a much  greater 
relative  change  in  composition  within  the  reactor.  Since 
measurement  of  the  variation  in  burnedness  at  high  burned- 
ness levels  is  not  a very  sensitive  measure  of  the  mixedness 
of  the  reactor,  Baker  proposed  to  measure  the  variation  of 
burnedness  within  the  reactor  at  a loading  factor  near 
blowout.  He  found  a greater  variation  of  burnedness  with 
position  as  might  be  expected  with  this  more  stringent  test, 
but  even  so,  the  reactor  still  appears  to  be  moderately 
well-mixed. 

Despite  the  evidence  of  the  effect  of  hole  size  and  of 
sampling  traverses  that  the  reactor  is  not  completely  mixed, 
the  stirred  reactor  still  remains  a useful  device  for 
studying  the  overall  kinetic  behavior  of  burning  gases  since 
the  evidence  indicates  that  the  mixing  is  rapid  enough  so 
that  chemical  reaction  rates  are  the  predominant  phenomena 
controlling  the  reactor  performance.  Thus,  using  the  blow- 
out data  of  the  stirred  reactor,  it  is  possible  to 
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evaluate  both  a pseudo-collision  constant  and  a pseudo- 
activation  energy  for  the  reaction.  These  quantities,  though 
representing  the  gross  reaction  rather  than  the  detail  of  the 
individual  stages  of  the  combustion  process  are  useful  in 
the  evaluation  of  various  fuels  and  in  the  prediction  of 
practical  combustor  performance. 

In  addition,  Longwell  and  Weiss  have  reported  sampling 
runs  in  which  combustion  efficiency  as  a function  of 
loading  factor  for  a constant  fuel-air  ratio  was  measured. 
These  data  showed  a poor  agreement  with  the  theoretical 
curve  for  a simple  bi-molecular  collision  process  with 
reaction  rate  constants  predicted  from  blowout  data. 

However,  other  than  this  one  bit  of  data,  no  work  appears 
to  have  been  done  on  the  operating  characteristics  of  the 
stirred  reactor  within  the  stable  combustion  region.  In 
the  application  of  stirred  reactor  data  to  the  prediction 
of  practical  combustor  performance,  it  is  necessary  to 
have  such  information  since  combustion  efficiency  as  well  as 
combustor  stability  is  Important  in  chamber  operation  and 
design.  Since  the  requisite  data  of  feed  rate  vs  efficiency 
are  not  available,  a calculated  curve  will  be  used  here. 

The  curve  is  the  one  shown  in  Figure  II-2.  ThlE  curve  was 
calculated  using  combustion  stoichiometry  in  which  propane 
is  burned  to  C0a  and  Ha0  in  one  step.  A stoichiometric 
mixture  was  postulated  with  an  inlet  temperature  of  500°K, 


a flame  temperature  at  complete  combustion  of  2772°K,  and 
an  activation  energy  of  40,000  cal/mol.  The  magnitude 
obtained  of  the  quantity  G/VP11  is  a consequence  of  the 
normalizing  procedure  used  In  calculating  this  curve,  and 
these  values  must  be  multiplied  by  a constant  factor  to 
make  the  blowout  value  correspond  to  the  observed  blowout 
rate  for  a stirred  reactor. 

Combustion  Chamber  Models  In  the  Literature 

Based  on  the  performance  of  the  stirred  reactor  various 

i 

simple  combinations  of  reactors  can  be  postulated  and  their 
behavior  calculated.  Some  models  of  this  type  have  already 
been  presented  and  discussed  In  the  literature.  Postulation 
of  such  models  may  have  as  objectives:  (l)  In  theoretical 
studies  one  may  attempt  to  find  the  ideal  or  optimum  per- 
I formance  which  can  be  expected  from  combustion  chambers,  and 

(2)  In  the  study  of  real  combustor  performance  one  may 
seek  to  find  a simple  combination  of  combustor  elements 

1 

which  will  describe  the  performance  of  the  chamber  voider 
study  with  respect  to  changes  in  the  various  operating 
or  design  variables.  Several  examples  of  the  former  type 
have  appeared  in  the  recent  literature. 

Avery  and  Hart  (2)  considered  a model  of  a can-type 
combustor  in  which  a pre-mixed  fuel-air  mixture  was  fed 
Into  the  can.  The  combustor  was  piloted  in  an  unspecified 
manner  by  a partly  burned  stream  of  the  same  fuel-air  ratio 
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as  the  fresh  feed  to  the  can  and  of  the  same  burnedness  as 
the  gases  leaving  the  can.  The  rate  of  feed  of  the  fresh 
fuel-air  mixture  along  the  length  of  the  combustor  was  so 
adjusted  that  the  composition  within  the  can  remained 
constant. 

Since  under  these  restrictions  the  composition  of  the 
gases  in  the  combustor  is  the  same  throughout,  the  burning 
rate  within  the  chamber  was  calculated  directly  from  well- 
stirred  reactor  principles.  In  addition  the  pattern  of 
fuel-air  mixture  injection  required  to  maintain  this  optimum 
condition  was  calculated.  Although  no  experimental  work 
has  been  reported  on  this  model  it  is  of  interest  since  it 
is  the  first  such  proposed  used  of  basic  stirred  reactor 
information. 

In  this  original  model  of  Avery  and  Hart  no  allowance 
was  made  in  the  energy  balance  for  the  effects  of  directed- 
veloclty  kinetic  energy.  In  a later  paper  Bosen  and  Hart 
(19)  proposed  a modification  of  this  model  to  account  for 
kinetic  energy  effects.  These  authors  considered  a 
straight,  tubular  combustor  piloted  by  a partially  burned 
gas  in  the  same  manner  as  in  the  Avery  and  Hart  model  and 
fed  with  fresh  fuel-air  mixture  through  the  combustor 
walls  in  such  a manner  as  to  maintain  a constant  static 
temperature  within  the  chamber.  These  authors  then  wrote 
differential  equations  of  continuity,  and  momentum  and 
energy  balances.  Prom  these  it  is  possible  to  derive: 


- 44 


(1)  the  maximum  throughput  of  unhurried  gases  permitted 
by  the  assumption  of  constant  static  temperature  (this 
Is  not  blowout,  but  analogous  to  "choking"  In  compressible 
flow);  (2)  the  relation  required  between  kinetic  energy, 
temperature  and  combustion  efficiency,  and  (3)  the  minimum 
temperature  (or  maximum  amount  of  enthalpy  which  may  be  put 
Into  kinetic  energy)  allowing  stable  combustion  for  a con- 
ventional nth  order  reaction. 

The  restriction  of  this  model  to  a constant  temperature 
case  and  the  restriction  of  the  Avery  and  Hart  model  to  a 
constant  efficiency  case  makes  these  models  of  relatively 
small  value  In  analyzing  the  behavior  of  real  combustors. 

beZubay  (6)  suggested  a different  approach  to  the 
modeling  problem.  He  considered  a combustor  to  which  a 
rich  fuel-air  mixture  was  fed,  and  after  combustion  was 
partially  completed,  additional  secondary  air  was  added. 

This  problem  was  approached  by  considering  the  chamber 
to  be  divided  into  two  portions,  each  a well-stirred  reactor. 
The  first  portion  was  fed  the  rich  fuel-air  mixture,  and  the 
second  received  the  partly  burned  products  from  the  first 
chamber  plus  secondary  air.  In  this  model  beZubay  demon- 
strated that  it  was  possible  for  two  stirred  reactors  to  give 
a higher  overall  heat  release  rate  than  a single  reactor 
fed  the  same  overall  fuel-air  mixture  and  operating  under 
the  same  conditions  as  the  final  reactor  of  the  pair. 
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However,  the  force  of  this  conclusion  is  somewhat  clouded 
by  the  effects  of  the  addition  of  secondary  air. 

The  results  which  DeZubay  presented  also  showed  the 
quenching  effect  of  adding  air  to  the  partly  burned  gases. 

It  was  shown  that  the  addition  of  a large  enough  quantity 
of  air  would  quench  the  reaction  in  the  second  reactor. 

Curves  were  given  of  heat  release  rate  vs  combustion 
efficiency  leaving  second  reactor  for  various  values  of 
fuel-air  ratio  in  the  first  of  the  two  reactors.  In  all 
cases  calculated  by  DeZubay  the  fuel-air  ratio  in  the 
second  stage  and  the  relative  volumes  of  the  two  stages  were 
the  same.  For  those  cases  where  a large  amount  of  secondary 
air  was  added  more  than  one  blowout  point  was  found;  i.e. 
as  feed  rate  increased  the  second  stage  became  unstable,  and 
further  increases  in  feed  rate  caused  the  first  stage  to 
blow  out.  It  was  also  found,  as  might  be  expected,  that  the 
best  operating  conditions  — maximum  heat  release  rates  — 
were  found  when  the  first  stage  had  a fuel-air  ratio  of 
stoichiometric . 

The  models  discussed  thus  far  have  treated  the  combus- 
tion chamber  problem  without  cognizance  of  the  effect  of 
flow  recirculation  on  combustion  chamber  performance.  Since 
recirculation  is  generally  postulated  to  have  an  important 
role,  especially  in  combustor  stability,  any  chamber  model 
which  is  to  be  representative  of  combustor  performance  must 

id 
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make  allowance  for  the  effect  of  recirculation.  In  a recent 
paper  Spaulding  (24)  has  presented  an  analysis  of  combustion 
chamber  pressure  drop  requirements  which  does  account  for 
the  effect  of  recirculation. 

Spaulding  considered  a combustion  chamber  in  which  the 
gases  flowed  in  a circulatory  manner.  Part  of  the  chamber 
was  devoted  to  a mixing  section  in  which  fresh,  unburned 
gases  mixed  with  recirculated  gases;  a portion  of  the 
mixture  was  recycled  and  a portion  left  the  chamber  as 
burned  products.  No  combustion  was  postulated  to  occur  in 
the  mixing  section.  The  remainder  of  the  cnamber  was 
called  the  recirculation  zone,  and  combustion  occurred  in 
this  region  as  the  gases  recirculated  in  a plug- flow  fashion. 

Using  an  assumed  stirred  reactor  performance  curve, 

Spaulding  could  then  calculate  the  performance  of  the  model 
in  response  to  changes  in  either  firing  rate  or  recirculation 
ratio.  This  model  is  very  similar  to  the  one  proposed  in 
this  work  for  modeling  combustion  chamber  performance. 

To  facilitate  comparison,  Spaulding's  model  has  been 
recalculated  using  the  stirred  reactor  curve  which  was  used 
through  this  work.  The  performance  of  Spaulding's  model 
is  shown  in  Figure  II-9. 

Although  no  experimental  work  was  done  directly  on 
this  model,  Spaulding  made  some  comparison  of  the  model 
performance  with  the  performance  of  some  combustor  - 
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FIGURE  11-9:  RECALCULATED  PERFORMANCE  OF 
SPAULDING'S  SIMPLE  MODEL  (£4) 
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configurations  in  his  electric  analogue  to  combustion  (23). 
Because  of  differences  between  the  model  and  combustor 
configurations,  it  was  not  expected  that  a good  correlation 
would  be  obtained,  however,  the  results  of  this  comparison 
were  encouraging. 

Spaulding  also  considered  an  application  of  this 
model  technique  to  the  well-stirred  reactor.  From  an 
assumed  flow  pattern  of  the  jet,  a model  was  established  in 
which  the  outwardly-flowing  cold  gases  received  recirculated 
gases  from  the  reverse-flowing  recirculating  gases  in  a 
continuous  manner.  The  recirculation  ratio  was  related 
in  this  model  to  the  down-stream  distance  from  the  Jet 
nozzle.  Spaulding  concluded  from  this  analysis  that  the 
well-stirred  reactor  design  suggested  by  Longwell  and  Weiss 
was,  in  fact,  capable  of  producing  moderately  well-mixed 
zones,  although  there  was  some  departure  from  the  Ideally 
well-mixed  reactor  performance. 

One  additional  model  of  combustor  performance  has  been 
proposed  by  Berl,  Rice  and  Rosen  (4) . They  analyzed  the 
volume  requirements  for  relatively  low  intensity  combustion 
in  terms  of  an  extended  laminar  flame  front  theory.  This 
model  is  not  closely  related  to  the  models  to  be  discussed 
here  and  will  not  be  considered  further. 
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Summary 

In  recent  years  It  has  become  evident  that  the  analysis 
of  hlgh-output  combustor  performance  must  take  Into  account 
the  rate-limiting  effect  of  chemical  reaction  as  well  as 
the  effect  of  mixing.  However,  It  was  not  until  the  concept 
of  grossly  simplified  kinetics  was  applied  to  combustion 
problems  that  the  required  analysis  became  practical.  The 
development  of  an  experimental  stirred  reactor  has  permitted 
the  gross  kinetic  constants  In  such  a simple  analysis  to  be 
estimated.  The  availability  of  stirred  reactor  performance 
data  has  tempted  some  investigators  to  devise  conceptual 
models  of  practical  combustor  performance  in  terms  of  stirred 
reactor  combinations.  With  one  exception,  however,  no 
allowance  has  been  made  for  the  effect  of  recirculation,  nor 
has  much  experimental  work  been  done  to  test  the  validity 
of  these  models.  In  this  thesis  some  new  models  of 
combustor  performance  are  proposed  and  compared  with 
experimental  combustor  performance. 
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CHAPTER  III 

MODELS  OP  COMBUSTION  CHAMBERS 

The  classical  analysis  of  combustion  chambers  has 
been  in  terms  of  flame  speeds.  Prom  the  spreading  angle 
of  a stabilized  flame  in  a duct  and  a knowledge  of  the 
stream  velocity  the  velocity  of  the  turbulent  flame  in 
its  normal  direction  is  calculated.  The  classical  analyst 
then  attempts  to  relate  this  turbulent  flame  speed  to  the 
laminar  flame  speed  of  the  combustible  mixture  and  to  the 
turbulence  level  in  the  chamber  which  is  assumed  to  cause 
extension  of  the  area  available  for  flame  propagation. 

In  more  complex  chambers  it  is  difficult  to  extend 
this  analysis,  since  as  the  heat  output  rate  of  a chamber 
increases  it  becomes  more  difficult  to  define  a meaningful 
boundary  to  represent  a time-average  flame  front.  As  one 
approaches  the  chemically  limiting  rates  in  a device  such 
as  a well-stirred  reactor  it  becomes  difficult  to  Justify 
the  use  of  an  extended  flame-front  model  since  substantial 
regions  of  the  combustor  are  given  over  to  homogeneous 
combustion. 

The  analysis  of  this  latter  class  of  chambers  can  be 


stated  in  terms  of  well-stirred  reactor  performance.  Using 
stirred  reactor  data  to  give  the  burning  rate,  or  fuel 
consumption  rate  per  unit  of  volume,  as  a function  of 


temperature  and  a variety  of  composition  variables,  the 
continuity  equation  for  a differential  segment  of  chamber 
volume  can  be  written.  The  time-rate  of  accumulation  of 
each  chemical  specie  plus  the  rate  of  creation  or  destruction 
of  that  specie  by  chemical  reaction  is  equated  to  the 
diffusive  and  convective  fluxes  of  the  specie  into  the 
chosen  control  volume.  An  equation  of  this  sort  must  be 
written  for  each  of  the  species  present  and  the  resulting 
set  must  be  solved  either  with  an  apriorl  knowledge  of  the 
fluid-flow  pattern  or  with  the  momentum  balance  differential 
equations  for  the  same  control  volume.  Considerable 
simplifications  of  these  equations  must  be  made  for  a 
practical  solution  even  in  the  simplest  of  chambers. 

As  an  alternative  analysis,  one  may  imagine  that 
relatively  large  regions  of  the  combustor  may  be  assigned 
average  compositions  which  will  give  the  burning  rate  in 
each  region,  and  that  these  compositions  may  be  related  to 
the  feed  rate  and  composition  of  gases  entering  each  zone. 
Prom  a knowledge  of  the  fluid-flow  pattern  the  interchange 
of  materiel  between  zones  may  be  evaluated  and  the  per- 
formance of  the  chamber  synthesized. 

A trivial  model  might  consist  in  identifying  the  modest 


performance  of  a practical  combustor  with  that  of  a well- 
stirred  reactor  operating  at  the  same  efficiency  and  firing 
rate  but  of  much  smaller  volume,  and  in  making  the 
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assumption  that  the  departure  of  the  practical  chamber  from 
well-stirred  reactor  performance  is  due  solely  to  Incomplete 
use  of  the  volume  available.  One  would  not  expect  so  simple 
a model  to  fit  over  any  extreme  range  of  operating  variables. 

If,  however,  one  were  to  postulate  a combustor  composed 
not  of  Just  a single  well-stirred  volume  and  a single 
Inactive  volume  but  of  a small  number  of  volumes  character- 
ized as  well-stirred  or  plug-flow  elements,  representing 
various  regions  of  the  real  chamber  connected  in  a 
combination  of  series,  parallel,  and  recirculating  flow 
patterns  which  is  representative  of  the  chamber  flow  pattern, 
one  could  expect  such  a conceptual  model  to  predict  the 
performance  of  a real  chamber  over  a much  wider  range  of 
operating  variables.  Simplicity  in  these  models  is  para- 
mount, for  only  if  the  model  consisted  of  a small  number 
of  interacting  elements  would  the  results  be  impressive. 

In  any  real  combustor  several  types  of  zones  are 
expected.  With  a Jet  of  fuel-air  mixture  feeding  a chamber, 
one  expects  the  central  core  to  be  a relatively  unburned 
zone.  Around  the  sides  of  the  Jet  the  more  nearly  burned 
gases  in  the  recirculation  zones  are  in  a well-mixed  region 
because  of  the  Intense  turbulence  created  by  the  central 
Jet.  In  the  corners  of  the  chambers,  regions  of  inactivity 
are  expected  since  these  regions  are  relatively  stagnant. 

If  the  fuel  and  air  are  fed  separately  into  a moderately 
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well-mixed  chamber,  similar  zones  are  expected.  In  this 
case  one  would  also  expect  to  find  variations  In  the  fuel/ 
air  ratio  in  the  various  portions  of  the  chamber. 

In  these  models  values  must  be  assigned  to  the 
relative  rates  of  Interchange  between  the  various  zones  and 
to  their  relative  sizes.  The  magnitude  of  the  variation 
In  the  fuel/air  ratio,  If  any,  must  also  be  chosen.  These 
values  are  associated  with  chamber  design  variables,  and 
to  the  extent  that  reasonable  values  may  be  associated  with 
the  theoretical  models,  the  models  may  be  held  to  be 
successful. 

Three  conceptual  models  have  been  developed  which 
predict  some  of  the  phenomena  of  practical  combustors. 

One  of  these  phenomena  Is  the  abrupt  change  In  flow  pattern 
observed  In  some  combustors  when  the  feed  rate  Is  increased. 
When  a certain  feed  rate  is  reached  the  flame  in  one  portion 
of  the  chamber  disappears,  and  there  is  an  accompanying 
Increase  in  the  volume  of  the  remaining  active  zones  at 
the  expense  of  the  zone  which  was  extinguished.  This 
behavior  has  been  observed  in  a two-dimensional  chamber 
fed  by  a mult-jet  system  of  a premixed  fuel/alr  mixture. 

A simplified  conceptual  model  of  such  a Jet-chamber 
combination  will  now  be  presented. 

Figure  III-l  shows  such  a chamber  with  the  fuel/air 
mixture  entering  a single  Jet  and  products  leaving  at  the 
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FIGURE  HL-3  WELL-STIRRED  REACTOR  CURVE  FOR 
TWO-STABILITY -LIMIT  MODEL 
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top.  Let  the  dotted  lines  bound  the  edges  of  those  areas, 
taken  normal  to  the  axis  of  flow,  through  which  the  total 
upward  gas  flow  is  constant  and  equal  to  the  feed  rate. 

The  volume  thus  bounded  is  called  the  "central  core",  and 
the  region  outside  of  this  volume  is  called  the 
"recirculation  zone".  It  should  be  noted  that  the 
recirculation  zone  defined  in  this  manner  includes  a 
complete  vortex  with  gases  flowing  both  upward  and  down- 
ward in  contrast  to  a definition  of  the  recirculation 
zone  as  that  region  where  there  is  only  downward  gas  flow. 
The  dotted  lines  defining  the  boundary  of  the  central  core 
represent  surfaces  through  any  element  of  which  there  is 
no  net  mass  flow,  but  there  is  equal  and  opposite  turbulent 
diffusion  through  them.  Because  the  feed  to  the  zones 
external  to  the  central  core  is  on  the  average  less  burned 
than  the  gas  which  returns  to  the  central  core,  the  mean 
gas  temperature  in  the  outer  zones  Is  higher  than  that  of 
the  central  core.  Consequently  the  linear  velocity  of  the 
outwardly  moving  gas  is  less  than  that  of  the  inwardly 
moving,  more  completely  burned  gas;  and  the  force  producing 
this  net  inward  momentum  flux  is  a higher  static  pressure 
existing  in  the  outer  zones. 

If  both  recirculation  zones  are  equally  active,  the 
pressure  effects  on  the  two  sides  counter-balance.  But 
if  one  zone  is  active  and  the  other  inactive,  a pressure 
Imbalance  occurs,  causing  a deflection  of  the  Jet  toward 


the  inactive  side.  This  increases  the  volume  of  the  active 
zone,  thereby  improving  the  opportunity  for  completion  of 
combustion  in  it.  The  end  effect  is  a bending  of  the  Jet, 
in  response  to  the  pressure  exerted  by  the  larger  recircu- 
lation zone,  to  a position  of  equilibrium  between  the 
Jet-momentum  force  and  the  pressure-area  force  of  the  hotter 
recirculation  zone. 

Since  the  system  is  two-dimensional  the  difference  in 
the  volumes  of  the  recirculation  zones  may  be  assumed  to  be 
proportional  to  the  pressure  difference  between  them,  and 
since  the  pressure  effect  is  postulated  to  be  the  result 
of  a tranverse  momentum  flux,  the  pressure  difference  be- 
tween the  two  zones  should  be  roughly  proportional  to  their 
temperature  difference.  To  simplify  this  treatment,  the 
central  core  will  be  assumed  to  be  a high-speed,  non- 
reacting Jet  stream  in  which  the  gases  reside  for  a 
relatively  short  time,  serving  to  supply  feed  gases  to  the 
recirculation  zones,  and  the  latter  zones  will  be  considered 
well-stirred  reactors. 

The  block  diagram  representing  this  model,  Figure  111-2, 
consists  of  two  well-stirred  reactors  of  constant  total 
volume,  separated  from  each  other  by  an  elastic  barrier 
zone,  the  Jet,  which  has  the  property  of  bending  in 
response  to  pressure  differences  so  that  the  change  in 
volume  of  either  of  the  reactors  is  directly  proportional 
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to  the  pressure  difference  between  the  two.  In  the  absence 
of  a pressure  difference  the  system  is  symmetrical.  The 
fuel-air  mixture  will  be  assumed  to  be  fed  equally  to  the 
two  reactors. 

The  deflection  of  the  barrier  is  expressed  in  terms 
of  the  difference  in  the  volumes  of  the  two  reactors; 


V - V. 


b _ 


- k>  - V 
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and  the  pressure  difference  created  by  a difference  in 
burnedness  in  the  two  reactors  is: 


P - P, 

a b 


<0.  - fyj 
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Assume  a well-stirred  reactor  performance  curve  is  available 
(e.g.  Figure  III-3) , giving  the  relation  between  burnedness 
and  normalized  volume. 

In  this  model  a series  of  stirred  reactor  pairs  of 
various  volume  differences  may  be  considered  at  each  firing 
rate.  For  each  volume  difference  considered  the  correspond- 
ing pressure  difference  between  the  two  reactors  may  be 
found  by  combining  Figure  II-3  and  Equation  III-2.  These 
pressure  differences  found  may  then  be  compared  with  the 
pressure  difference  required  to  maintain  the  postulated 
volume  difference  from  Equation  III-l.  For  each  new  firing 
rate  considered  a new  series  of  volume  differences  and 
pressure  differences  may  be  found. 


- 58 


This  calculation  may  be  represented  graphically,  bet 
the  total  volume/feed  ratio,  equal  to  the  Initial  ratio  for 
each  reactor,  be  represented  by  point  "0"  of  Figure  III-3, 
and  let  the  barrier  be  deflected  by  an  amount  (V  - V^i/2, 
causing  zones  * and  t to  Increase  and  decrease  respectively 
by  this  amount.  These  volume  changes  are  accompanied  by 
corresponding  changes  In  burnedness  which  may  be  found 
from  Figure  1II-3  where  a.  and  b represent  the  new  volume/ 
feed  ratios  of  sides  a and  b.  Solutions  for  the  burnedness 
In  the  new  volumes  exist  on  the  branches  at  a»  fi." » and 
b and  b"  (a*  and  b*  are  thermally  unstable  reactor  conditions). 

These  burnedness  differences,  4 of  them  (plus  5 involving 
unstable  reactor  conditions) , can  be  evaluated  and  from 
Equation  III-2  the  corresponding  (Pft  - P^) 's  can  be  found. 

These  points  for  a typical  case  are  plotted  on  Figure  III-4. 

(The  designations  ak»  ak'  etc.  represent  reactor  conditions 
of  Figure  III-3.)  depitition  of  this  calculation  for  a 
series  of  volume  differences  leads  to  branches  I and  11 
of  Figure  III-4.  Those  portions  of  the  curves  which  Involve 
thermally  unstable  reactor  conditions,  i.e.  the  dotted 
portions  of  Figure  III-3,  are  included  for  completeness  but 
are  dotted.  Figure  1II-3  shows  rather  inadequately  that 
Branch  II  Is  symmetrical  around  the  origin;  no  symmetry 
around  either  axis  Intended. 

The  solution  occurs  when  both  reactors  are  hot, 

i . iJ 


FIGURE  m - 4:  PERFORMANCE  OF  A TWO-STABILITY- LI  MIT 

COMBUSTION  CHAMBER  MODEL.  GROSS  VOLUME- 
FEED  RATIO  AND  FEED  COMPOSITION  FIXED. 
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and  a low  P#  - Is  developed,  lying  on  Branch  11  in  the 
first  quadrant.  If  side  is  blown  out,  a high  P^  - P^ 
is  developed  at  the  point  ab"  on  Branch  I,  also  in  the  first 
quadrant.  Similar  reasoning  leads  to  the  remainder  of  the 
points. 

For  each  new  value  firing  rate  chosen  a new  pair  of 
branches  I and  II  will  be  found.  It  will  be  noted  that 
at  a sufficiently  high  firing  rate  the  reactor  pair  will 
Just  be  able  to  operate  with  both  sides  active.  This  occurs 
when  point  "0"  of  Figure  IH-3  is  Just  at  the  blowout  point 
of  the  well-stirred  reactor  curve.  If  the  feed  rate  to  the 
same  reactor  paid  is  increased  further,  it  may  still  be 
possible  for  one  side  of  the  model  to  operate  on  the  stable 
"hot"  branch  of  the  reactor  curve  if  the  volume  difference 
is  large  enough;  however,  since  the  maximum  volume  difference 
possible  occurs  when  V is  twice  its  initial  value,  the 
maximum  feed  rate  which  permits  a possible  "hot"  reactor  is 
twice  the  feed  rate  corresponding  to  the  above  blowout 
point. 

To  find  the  stable  operating  points  on  Figure  1II-4, 
the  line  representing  Equation  III-l  is  superimposed  on 
branches  I and  II  for  the  specific  "0"  chosen.  (A  non- 
linear relation  between  V - V.  and  P_  - PK  could  be  used 

a b a d 

equally  well).  The  equilibrium  values  of  are  found 

at  the  intersections  of  the  line  of  Equation  III-l  and  either 
branch  I or  branch  II.  Branches  I and  II  and  the  curve  of 


- 61  - 


equation  III-l  are  symmetrical  around  the  origin,  i.e.  it 
makes  no  difference  which  side  becomes  colder.  Three 
stable  solutions  are  found,  the  origin  and  the  outermost 
intersections,  A,  D,  and  D* . The  remaining  four  inter- 
sections are  all  unstable;  C and  C'  because  one  side  of 
the  reactor  pair  is  operating  on  the  unstable  branch  of 
the  well-stirred  reactor  relation,  and  B and  B'  because  the 
barrier  is  in  a metastable  position.  A stable  barrier 
position  is  one  in  which  a momentary  change  in 
causes  a change  in  Pft  - Pfc  such  that  a net  restoring 
force  results.  This  means  that  at  a stable  intersection 
of  the  branches  I or  II  with  the  line  of  Equation  III-l 
both  curves  must  have  a positive  slope  and  the  slope  of 
the  former  must  be  greater  than  the  slope  of  the  latter. 

The  physical  significance  of  the  constants  in 

Equations  III-l  and  III-2  is  of  importance  in  relating 
the  theory  to  experimental  and  design  problems.  The 
interpretation  of  these  constants  may  be  seen  by  reducing 
Equations  III-l  and  III-2  to  dimensionless  relations. 

In  any  chamber  which  has  a pressure  difference  trans- 
verse to  the  flow  direction  of  a Jet  a deflection  of  the 
Jet  will  be  found.  This  deflection  would  be  expected  to  be 
a function  of  the  ratio  of  pressure-area  forces  to  the  Jet- 
momentum  force,  and  assuming  dimensional  similitude,  this 
function  may  be  expressed: 
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Ayr  * ^(APri  xr) 


II 1-3 


where  Ayr  * Ay/yN,  xr  = x/xN,  and  APr  * AP/P^*  • 
yjj  - nozzle  width,  nozzle  velocity, 


PN  = nozzle  density. 

AP  = the  difference  in  pressure  on  the  two  sides 
of  the  Jet  averaged  from  the  nozzle  to  the 
distance  x from  the  nozzle. 

Ay  = the  deflection  of  the  Jet  centerline  at  the 
distance  x from  the  nozzle. 

The  change  in  volume,  V - VQ,  for  either  side  may  be  computed 
by  integration  of  the  deflection  of  the  Jet  over  the  length 
L of  the  chamber: 


V - V = 
o 


■/ 


Ay  dx  = wyjj  | Ayr  dx  r 


wyjj 


I 

/ 


(AP  ,X  ) dx^ 
vi  r r 
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This  reduces  to: 


- «yj  $UP*r  , L.)  WI-6 

c* 

where  APr  is  the  reduced  pressure  difference  averaged  over 
the  entire  chamber  length  and  w is  the  width  of  the  two- 
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dimensional  Jet.  This  may  be  compared  with  Equation  III-l 
directly,  and  to  be  in  conformity  with  it,  i.e.  to  make  the 
displacement  proportional  to  AP^  , Equation  III-5  must  be 
of  the  form: 


V - V. 


wyN  APV  F (LrJ 
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This  equation  may  be  rewritten  in  dimensionless  form: 


V.  - 


V + V.  V + Vv 
a b a b 


AP*  P (L  ) 
r r 
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The  "jet*  will  now  be  defined  more  precisely  as  that 
part  of  the  flow  pattern  of  constant  axial  mass  flow  rate 
(i.e.  the  same  portion  of  the  flow  referred  to  in  Figure 
III-l) ; and  that  part  of  the  system  external  to  the 
boundary  surfaces  of  this  special  Jet  will  be  called  the 
"recirculation  zone".  (Since  that  part  of  the  zone  in 
contact  with  the  Jet  core  has  a large  positive  axial 
velocity  component,  this  is  not  the  conventional  Jet 
referred  to  in  the  literature  on  Jet  structure.) 

Because  of  the  definition  of  the  recirculation  zones, 
these  zones  must  be  fed  with  fresh  combustible  gases  by  an 
eddy  diffuslonal  process,  and  burned  gases  must  leave  them 
similarly.  These  diffusion  terms  may  be  assigned  average 
velocities  Uq  outward  and  inward  which  are  applicable 
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to  the  entire  interface  between  the  Jet  mainstream  and 
the  recirculation  zone,  and  by  the  similarity  principle: 

uo  * aUN  1II-8 

Using  this  outward  velocity,  the  outward  momentum  flux  mQ 

associated  with  it  can  be  evaluated: 

m = Prt  U ■ = a8P  U*  I I 1-9 

O 0 0 o « 

By  applying  a material  balance,  the  gases  leaving  the 
recirculation  zone,  P^  U^,  must  be  equal  to  the  gases 
entering  that  zone: 


PoUo  - P1U1 
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The  Inward  momentum  flux,  m^ , is  then 

“l  “ piu!  - <poUS>  pc/pl  - **  poUN  po/pl  II1-11 


Hence  the  net  momentum  flux  at  the  boundary  is: 

si  = «»P0ujS  - 1) 
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and  for  two  zones,  the  pressure  difference  between  them  may 
be  calculated,  assuming  there  is  negligible  temperature 
rise  in  the  Jet  mainstream,  by  the  equation; 
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Now  by  an  energy  balance: 


tm  ■ " ( 


+ 1 


111-14 


- 65  - 


Hence: 


ap.  . ,» 


f--1) 

iN  • 
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This  equation,  when  compared  with  Equation  III-2  is  of 
identical  form,  i.e.  the  pressure  difference  is  directly 
proportional  to  the  burnedness  difference. 

Figure  III -4  may  be  replotted  using  coordinates  of 
(Vft  - V^)/ + Vfc)  and  , yielding  a generalized 

expression  of  the  relative  performance  of  the  two  sides  of 
the  model.  In  this  diagram  all  possible  combinations  of 
and  must  lie  within  the  square  bounded  by 

the  limits  of  (V^  - + * ± 1 and  Pa  - (3^  = + 1. 

On  this  same  diagram  the  response  of  the  Jet  to  a burnedness 
difference  may  also  be  plotted,  this  relation  being  obtained 
by  a combination  of  Equations  III-7  and  III-15: 

-» ~ Vfe  - — — ?<v]  <?«-  pb> 

v.  +vt  L \ +vb  J 

II 1-16 

This  plot  of  the  relations  embodied  in  Equations 
III-l  and  III-3  has  the  advantage  that  all  of  the  geometrical 
design  variables  are  lumped  into  one  coefficient,  namely  the 
slope  of  the  line  representing  III-16.  Inspection  will 
verify  that  for  a given  value  of  the  overall  loading  factor, 
G/fV^  + V^JP11  » the  relative  performance  of  the  two  reactors 
of  the  model  (branches  I and  II  of  Figure  III-4)  when 


plotted  on  these  coordinates  will  yield  the  same  curve, 
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independent  of  the  shape  of  the  chamber.  Different  curves 
of  branches  I and  II  will  be  obtained  for  each  value  of 
the  loading  factor  chosen.  However,  from  Equation  III-16, 
the  response  of  the  Jet  to  the  burnedness  differences  is 
independent  of  the  velocity  (or  loading  factor).  Hence, 
in  comparing  the  response  of  a given  chamber,  fed  a given 
fuel/air  ratio  and  fuel  type,  to  changes  in  the  feed  rate 
a family  of  curves  of  branches  I and  II  will  be  obtained 
which  intersect  the  one  line  representing  Equation  III-16. 

Design  variables  such  as  chamber  and  nozzle  shape,  on 
the  other  hand,  are  reflected  only  in  the  slope  of  the 
line  representing  Equation  III-16.  These  variables,  together 
with  the  volumetric  expansion  of  the  burned  fuel/alr  mixture, 
(T^f/Tjj)  - 1,  determine  the  fractional  volume  change  produced 
by  Jet  movement  in  response  to  a unit  of  burnedness  difference 
in  the  two  sides.  In  Equation  III-16  the  variables  affected 
by  chamber  design  are  the  constant  aa  which  is  related  to 
the  Jet  spreading  angle  and  the  quantity  (fewy^)/  (V^  + V^Pd^) 

which  is  related  to  the  chamber  geometry.  Changes  in  the 
feed  composition  and  inlet  temperature  will  affect  both 
branches  I and  II  (Bince  this  will  change  the  basic  stirred 
reactor  relation  shown  in  Figure  III-3)  and  the  Jet 
sensitivity  through  <T^ - 1. 

Using  a diagram  such  as  Figure  III-4  in  the  generalized 


form  discussed  above,  both  the  response  of  the  chamber  to 
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changing  feed  rates  and  the  stability  characteristics  of 
the  chamber  may  be  Investigated.  Consider  the  chamber  to 
be  burning  Initially  In  a stable,  symmetrical  pattern 
(operation  corresponding  to  the  origin  of  Figure  III-4) . 

If  the  total  loading  factor  of  the  chamber  is  increased, 
point  "0"  of  Figure  III-3  will  move  toward  the  p axis. 

As  it  approaches  the  nose  of  the  curve  of  Figure  III-3, 
branch  II  of  Figure  III-4  becomes  smaller  and  thinner 
(i.e.  its  slope  at  the  origin  becomes  smaller  and  its 
Pa  - Pb  intercept  approaches  the  origin).  This  change  in 
operating  conditions  will  not  affect  the  slope  of  the 
curve  representing  jet  position.  When  the  loading  factor 
has  increased  sufficiently  to  make  branch  II  become  tangent 
to  the  Jet  position  curve  at  the  origin  (Figure  III-5), 
the  origin  becomes  an  unstable  intersection  of  the  jet 
position  line  and  the  reactor  operating  lines.  Hence 
the  only  remaining  stable  solution  is  the  intersection  of 
the  Jet  position  line  with  branch  I of  Figure  III-5, 
corresponding  to  one  side  burning  stably  and  the  other 
blown  out.  This  calls  for  a sudden  deflection  in  the  jet 
to  the  new  stable  position  of  - V^.  Further  increase 
in  the  feed  rate  will  cause  the  disappearance  of  branch  II, 
and  eventually  branch  I will  become  tangent  to  the  line  of 
Equation  III-16  (Figure  III-6) . This  represents  the  condi- 
tion for  total  blowout  of  the  chamber,  since  at  this  point 


FIGURE  HI-  5t  ONE-SIDED  BLOWOUT  CONDITION  IN  THE 
TWO-STABILITY -LIMIT  MODEL 


FIGURE  JE-  6:  FINAL  STABILITY  LIMIT  FOR  THE  TWO- 
STABILITY  - LIMIT  MODEL 
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and  higher  feed  rates  there  are  no  stable  intersections  of 
branch  I of  Figure  III-6  with  the  line  of  Equation  III-16. 

This  model  also  has  interesting  implications  con- 
cerning behavior  at  other  chamber  feed  rates.  Figures  III-7a 
to  III-7f  represent  six  feed  rates  with  six  corresponding 
pairs  of  branches  I and  II,  along  with  a common  curve  showing 
the  Jet  stiffness  independent  of  feed  rate.  Here  only  the 
stable  portions  in  the  first  quadrant  of  Figure  III-4  are 
shown.  On  each  of  these  diagrams  a horizontal  dotted  line 
is  drawn  tangent  to  the  nose  of  branch  II.  This  represents 
the  blowout  point  of  branch  II,  i.e.  the  maximum  volume 
difference  allowable  for  a given  volume/  feed  ratio  if  both 
zones  are  to  be  simultaneously  stable.  If  exceeds 

this  level  the  smaller  side  will  be  blown  out,  and  if 
Va  - is  below  this  limit,  either  the  smaller  side  will 
burn  stably,  or  if  it  is  blown  out  it  may  be  relighted  by 
a sufficiently  strong  source. 

Consider  a chamber  operating  at  a volume/feed  ratio 
represented  by  Figure  III-7b  with  both  sides  active,  i.e. 
with  operation  at  the  origin  of  Figure  III-7b.  If  a dis- 
turbance, , is  imposed  on  the  Jet,  the  Jet  will 

• nd  to  return  to  its  central  position  unless  the 

• is  greater  than  the  blowout  limit  represented  by 

■>  • • •■■  . ; n» . If  tnls  limit  is  exceeded  then  one  side 

. **  . . «x.<i  a high  P#  - Pfa  will  be  developed. 


The  Jet  will  consequently  move  to  the  one-sided  blowout 
position  (the  branch  I Intersection)  Instead  of  returning 
to  its  central  position.  If  the  chamber  is  operating  at 
the  one-sided  blowout  position  when  the  disturbance  is 
Imposed,  then  the  reverse  argument  is  applied;  however, 
in  this  case  the  problem  of  re-light  must  be  considered. 
Even  if  the  disturbance  is  large  enough  to  bring  the  small 
side  to  a re-lightable  position,  the  small  side  will  not 
relight  if  the  ignition  source  is  too  weak;  and  the  Jet 
will  return  to  the  one-sided  blowout  position.  Since  in 
this  case  the  disturbance  required  to  reach  the  critical 
level  1b  much  larger  than  in  the  first  case,  it  would  be 
expected  that  the  one-sided  blowout  position  would  be 
the  more  stable  one. 

At  a lower  feed  rate  than  in  Figure  III-7b,  such  as 
represented  by  Figure  III-7c,  the  blowout  point  of  branch 
II  is  about  centrally  located  with  respect  to  the  two 
stable  Jet  positions.  In  a highly  turbulent  combustor  in 
which  there  are  high  transfer  rates  from  active  zones  to 
extinguished  zones  (conditions  typical  of  many  high-output 
combustors)  it  is  easy  to  imagine  that  if  the  volume/feed 
ratio  for  any  cold  zone  is  above  the  minimum  required  for 
stability  by  the  well-stirred  reactor  kinetics,  ignition 
will  occur.  If  simultaneously  there  is  an  equal  liklihood 

that  a disturbance  in  V - V,  from  either  stable  Jet 

a b 
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poBition  will  cross  the  critical  blowout  level,  conditions 
for  a flickering  flame  exist. 

At  even  lower  feed  rates,  represented  in  Figure  III-7d, 
the  situation  is  reversed  from  that  represented  by 
Figure  III-7b.  In  Figure  III-7d  the  two-act lve-zone  solution 
is  the  relatively  more  stable  since  a smaller  disturbance 
is  required  to  go  from  the  one-sided  blowout  position  to 
the  maximum  relight  level  than  in  the  opposite  direction. 

In  the  sixth  diagram,  Figure  III-7f,  the  feed  rate 
is  so  low  that  the  one-sided  blowout  solution  is  below  the 
critical  blowout  level  of  branch  II.  Thus,  assuming  that 
relight  will  always  occur  if  the  opportunity  exists,  the 
only  stable  solution  at  this  rate  is  the  origin  (two  active 
zones) . 


This  type  of  argument  also  leads  to  the  conclusion 
that  a hysteresis  effect  exists  in  this  combustion  chamber 
model.  In  reaching  the  transition  from  two-active  zones 
to  the  one-sided  blowout  from  low  feed  rates,  chamber 
operation  moves  from  a low  feed  rate  condition  represented 
by  Figure  III-7f  to  that  high  flow  rate  condition  repre- 
sented by  Figure  III-7a,  where  the  origin  becomes  an 
unstable  solution.  Once  the  transition  occurs,  relight 
can  be  caused  only  by  reducing  the  feed  rate  to  a point 
where  the  one-sided  blowout  solution  is  below  the  blowout 
point  of  branch  II.  The  feed  rate  where  this  critical 
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level  is  reached  is  represented  in  Figure  III-7 
Thus,  it  may  be  concluded  that  this  simple  model  of  a com- 
bustion chamber  is  useful  not  only  for  studying  pressure 
interactions  of  two  or  more  zones  but  also  for  suggesting 
mechanisms  for  other  important  chamber  problems. 

Model  Showing  Effect  of  Heclrculatlon 

Flow  recirculation  is  another  phenomenon  of  combustion 
chambers  which  is  characteristic  of  many  systems.  This 
has  been  recognized  by  Spaulding  (24)  in  a paper  which 
appeared  after  the  present  model  had  been  developed. 
Spaulding  proposed  a model  of  a recirculating  system 
similar  to  the  model  developed  here,  although  the  treat- 
ment and  conclusions  are  markedly  different.  In  a chamber 
fed  by  a Jet  of  fuel/air  mixture  the  region  of  the  Jet 
near  the  base  of  the  chamber  aspirates  burned  gases  from 
the  surrounding  region.  These  gases  are  replaced  by 
reverse-flow  from  near  the  exit  of  the  chamber. 

In  constructing  a model  of  this  behavior,  Spaulding 
suggested  that  the  gases  flow  around  this  loop  in  plug- 
flow.  The  central  core  serves  only  to  mix  the  recirculated 
gases  with  the  fresh  feed.  A portion  of  this  mixture  is 
discharged  from  the  model  as  burned  products,  and  the 
remainder  recirculates  in  a plug-flow  combustion  zone. 
Spaulding  has  calculated  two  cases  for  this  model;  one 
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where  there  1b  no  interchange  between  the  Jet  mainstream 
and  the  recirculation  zone  (all  recirculated  material 
enters  at  the  base  of  the  Jet)  and  a second  case  where 
allowance  is  made  for  continuous  aspiration  of  burned 
gases  from  the  recirculation  zone  along  the  length  of 
the  Jet.  The  first  of  these  cases  is  directly  comparable 
to  the  results  obtained  here,  and  the  behavior  of  this 
model  is  shown  in  Figure  11-9. 

No  allowance  has  been  made  in  Spaulding's  model  for 
combustion  occurring  in  the  central  core  of  the  Jet.  A 
model  in  which  such  an  allowance  is  made  is  probably  in 
better  agreement  with  the  actual  situation.  The  con- 
sequences of  allowing  for  combustion  in  the  central  core 
are  illustrated  by  the  model  considered  here.  In  this 
model  both  the  central  core  and  the  recirculation  zone 
have  been  treated  as  well-stirred  reactors.  Figure  III-8 
shows  this  model  dlagrammatically . Fresh  feed  enters 
and  is  mixed  instantaneously  with  recycled  burned  gases. 

The  mixed  gases  then  enter  zone  Vi,  representing  the 
central  core  of  the  Jet.  A portion  of  the  partially  burned 
products  leaving  Vi  is  recycled  through  the  recirculation 
zone  Va,  where  further  combustion  occurs,  and  then  is 
returned  to  Vi  again. 

To  find  the  feed  rate  to  each  of  the  zoneB,  a 
recirculation  ratio  is  defined  as: 


- 76 


B B Mplal  recirculation.  rate_to_zoBe_ya 
Wolal  fresh  feed  to  model 

This  definition  Is  the  same  as  that  used  by  Spaulding. 

From  a performance  curve  showing  the  relation  between 
feed  rate  and  burnedness  (i.e.  Figure  III-3),  the  response 
of  this  system  to  feed  rate  changes  may  now  be  calculated 
for  a given  ratio  of  Va/Vx  and  recirculation  ratio  B. 

This  calculation  has  been  carried  out  for  a fixed 
value  of  Va/Vx  * 2.0  and  for  various  recirculation  ratios. 
Burnedness,  as  a function  of  feed/total  volume  is  shown 
in  Figure  II 1-8. 

Several  performance  characteristics  of  this  system 
should  be  noted.  If  the  chamber  is  operating  at  a recir- 
culation ratio  of  zero,  the  volume  Va  allowed  for  recir- 
culation represents  "dead"  space,  and  the  central  core  Vx 
has  only  fresh,  unburned  gases  entering  it.  The  central 
core,  therefore,  operates  as  a simple  well-stirred  reactor 
without  any  additional  assistance  from  the  recirculated 
gases.  The  curve  representing  the  behavior  of  zero 
recirculation  may  be  found  by  multiplying  the  curve  for 
the  simple  well-stirred  reactor  by  the  fraction  of  the 
volume  in  the  model  which  represents  the  active  central 
core,  (i.e.  by  Vi/(VX  + Va);  curve  E * 0 of  Figure  III-8) . 

If,  on  the  other  hand,  the  recirculation  ratio  is 
infinite,  the  composition  of  both  zones  is  identical,  and 
the  performance  becomes  that  of  the  well-stirred  reactor 
performance  used  to  generate  all  of  the  other  curves. 


77 


(Curve  R = <co,  Figure  III-8) . 

For  recirculation  ratios  neither  zero  nor  Infinite, 
significant  interaction  occurs  between  the  central  core 
and  the  recycle  volume.  This  is  Illustrated  by  considering 
the  relative  contributions  of  each  of  the  reactors  to  the 
total  amount  of  combustion.  Let  the  burn  Iness  of  the 
streams  leaving  reactors  Vx  and  V3  be  represented  by  px 
and  p3  (the  former  also  represents  the  burnedness  of  the 
gases  entering  V3  and  leaving  the  system) ; and  let  p0 
represent  the  burnedness  of  the  gases  fed  to  Vx,  obtained 
by  mixing  gas  of  composition  Pa  with  fresh  feed.  The 
contribution  to  burnedness  of  each  reactor  may  be  expressed 
quantitatively  as  the  flow  rate  through  the  reactor  times 
the  increase  in  burnedness  in  that  reactor: 
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By  dividing  through  by  the  total  burning  rate  of  the  system, 
G*Px  , two  terms  are  obtained  representing  the  fractional 

contributions,  J , of  each  reactor: 

Pi  - P„  . 


-2-  (1  + R) 


- Pi 


78  - 


/r 


Choose  an  operating  point  corresponding  to  a moderately 

low  feed  rate,  and  consider  the  response  of  a fixed- 

volume  reactor  pair  to  an  increase  in  feed  rate. 

Figure  III-9  shows  the  behavior  of  various  points  in 

the  model  for  a recirculation  ratio  of  1/2  and  a 

Va/V i of  2.  Curve  pi  is  the  burnedness  in  the  reactor  Vx 

and  also  the  burnedness  leaving  the  system  (this  curve  is 

also  K = 1/2  of  Figure  III-8) . Curve  Pa  is  the  burnedness 

in  the  recirculation  zone:  and  curve  B is  the  burnedness 

o 

of  the  feed  to  Vi  obtained  by  mixing  fresh,  unburned  feed 
with  recirculated  products.  At  a low  feed  rate 
(represented  by  dotted  line  "A")  the  feed  material  to  Va, 
of  composition  Pi,  has  a high  burnedness  and  there  is  very 
little  burning  remaining  to  be  accomplished  in  reactor  Va. 
Under  these  conditions,  the  fractional  utilization,  Ja, 
of  the  recirculation  zone  is  low,  i.e.  most  of  the  burning 
occurs  in  Vx.  At  a higher  feed  rate,  represented  by  line 
"C"»  the  composition,  when  approached  from  low  flow  rates, 
corresponds  to  points  1,  3,  and  5;  and  most  of  the  burning 
still  occurs  in  Vi.  However,  since  the  burnedness  in  Vx 
has  decreased  slightly,  Ja  increases  correspondingly.  As 
the  feed  rate  increases  further  past  line  *D",  the  central 
core  becomes  chilled,  and  blowout  occurs,  after  which  only 
the  recycle  volume  is  active  (i.e.  a residual  flame).  At 
feed  rates  beyond  line  "D"  combustion  occurs  only  along 
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FIGURE  m-10: 


EFFECT  OF  SPLIT  FUEL  FEED  ON 
RECIRCULATION  MODEL 
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the  lowest  branch  of 
line  "C").  At  these 
combustion  occurs  In 
unity . 


the  curve  (corresponding  to  dotted 
feed  conditions  practically  no 
the  central  core,  and  Ja  Is  nearly 


If  the  flow  rate  is  now  decreased  from  this  high 
level,  the  central  core  stays  cold  past  line  "D"  and  at 
line  " C * the  compositions  2,  4,  and  6 are  reached. 

When  line  "B"  Is  reached  there  is  a sudden  relighting  of 
the  central  core.  This  hysteresis  loop  is  the  most 
Interesting  part  of  these  models,  since  it  is  a potential 
source  of  a flickering  flame.  In  this  region  the  central 
core  may  light  or  blowout  in  response  to  small  disturbances 
in  operating  conditions  with  very  little  change  in  the 
burnedness  of  Va. 

At  higher  recirculation  ratios  the  behavior  is  similar 
but  because  of  the  "stirring"  effect  of  higher  recycle 
rates,  the  ratio  of  the  fractional  contributions  of  the 
two  volumes,  Ji/Jz,  approaches  unity  for  a wider  range  of 
feed  rates.  As  Figure  III-8  shows,  the  undercut  of  the 
curve  Pi  disappears  at  higher  recirculation  ratios,  and 
as  the  recycle  is  continually  increased,  the  inflection 
in  the  curves  becomes  less  pronounced  until  it  also 
disappears. 
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Model  Showing  Effect  of  Non-Uniform  Fuel/Air  Batlo 

Another  problem  of  interest  is  that  of  allowing  for 
the  effect  of  unequal  fuel/air  ratios  in  the  various 
chamber  zones.  Such  a condition  would  arise  in  a comoustor 
where  the  fuel  and  air  were  injected  separately  into  a 
moderately  well-mixed  system.  A possible  model  of  such  a 
combustor  is  a variation  of  the  above  recirculation  model, 
in  which  only  a portion  of  the  fuel,  but  all  of  the  air, 
is  fed  into  the  central  core  and  the  remainder  of  the 
fuel  is  fed  into  the  recirculation  zone.  This  produces 
a richer  mixture  in  the  recycle  volume  than  in  the  central 
core.  The  one  additional  variable  required  may  be  thought 
of  as  either  the  fraction  of  the  original  fuel  introduced 
into  the  central  core  or  the  fuel/alr  ratio  of  the 
recirculation  zone.  The  former  has  the  advantage  of 
showing  the  limitations  on  the  variation  of  the  fuel/air 
ratio  of  the  recycle  volume.  The  latter  shows  more 
directly  the  effect  of  the  split  fuel  feed.  Since  the 
fuel/air  ratio  in  the  central  volume  is  the  same  as 
the  overall  fuel/air  ratio,  it  will  not  be  affected  by 
this  split. 

The  results  of  calculations  made  on  this  recircula- 
tion model  are  shown  in  Figure  III-10  for  an  overall 
fuel/air  ratio  which  is  stoichiometric  for  consistency 
with  the  previous  models.  Because  of  this  condition, 
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other  fuel/air  ratios  in  the  recirculation  zone  can  only 
be  a detriment  to  the  overall  heat  release  rate  since  the 
recycle  zone  must  operate  on  less  than  the  optimum  effi- 
ciency curve.  At  low  feed  rates  this  effect  is  of  little 
consequence  since  the  gases  in  the  recirculation  zone  are 
essentially  completely  burned  in  either  case.  However, 
if  the  feed  rate  is  increased  sufficiently  the  burning 
rate  in  the  recirculation  zone  becomes  significant,  and 
the  consequence  of  operating  this  volume  at  other  than 
stoichiometric  fuel/air  ratio  is  apparent. 

If  the  chamber  were  operating  with  a lean  overall 
fuel/air  ratio  and  the  fuel/alr  ratio  in  the  recirculation 
zone  were  stoichiometric  the  reverse  situation  would  be 
expected.  In  this  case  the  recirculation  zone  would  be 
operating  on  a more  favorable  efficiency  curve  than  the 
central  core.  In  those  feed  rate  ranges  where  the  per- 
formance of  the  model  for  an  overall  lean  fuel/air  ratio 
would  be  higher  than  for  the  case  of  feeding  all  the 
fuel  into  the  central  core. 

Models  of  More  Complex  Situations 

Thus  far,  the  models  discussed  have  been  restricted 
to  specific  phenomena.  In  modeling  practical  situations, 
however,  one  may  wish  to  consider  situations  where  two  or 
more  phenomena  are  simultaneously  significant.  If  models 
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for  the  separate  phenomena  of  interest  are  available, 
models  of  such  situations  may  be  built  from  simple 
combinations  of  the  individual  models.  Thus,  for 
instance,  in  a situation  where  both  residual  flame 
and  two-stability -limit  phenomena  are  observed,  it  may 
be  postulated  that  the  chamber  consists  of  two  or  more 
recirculation  zones  as  represented  by  the  two-stability 
limit  model  and  that  the  structure  of  these  zones  is 
represented  by  models  of  the  recirculation  type  rather 
than  a simple  well-stirred  reactor. 

Summary 

It  has  been  hypothesized  that  stirred  reactors  may 
be  used  in  simple  combinations  to  predict  the  behavior  of 
combustion  chambers.  Examples  considered  in  this  chapter 
have  shown  that  behavior  characteristics  found  in  practical 
combustion  chambers  may  be  duplicated  by  such  models  of 
interacting  elements.  Models  illustrating  pressure  inter- 
actions, recirculation,  and  non-uniform  fuel/alr  ratio 
effects  have  been  considered.  Although  the  calculated 
performance  curves  given  here  have  been  based  on  a 
theoretical  well-stirred  reactor  curve,  one  could  in 
principle  use  an  experimentally  determined  curve  of 
efficiency  vs  feed  rate  for  a well-stirred  reactor. 
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The  criterion  of  success  applied  to  these  models 
is  that  only  a relatively  few  variables  be  required  to 
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specify  a model  performance  in  conformity  with  reality  and 
that  these  variables  have  physical  significance.  These 
models  have  shown  that  it  is  possible  under  these  restric- 
tions to  simulate  a variety  of  known  combustion  chamber 
phenomena.  It  is  now  proposed  to  compare  quantitatively 
the  model  performance  with  the  performance  of  a simple 
experimental  combustion  chamber. 


CHAPTER  IV 

EXPERIMENTAL  EQUIPMENT 
AND  PROCEDURES 

The  experimental  equipment  used  in  this  work  con- 
sisted primarily  of  a specially  constructed  combustion 

chamber  of  a simple  design  to  facilitate  the  comparison 
with  the  theory.  Ancillary  equipment  Included  metering 
devices,  a gas  sampling  train  and  gas  analysis  equipment. 

Combustion  Chamber 

In  the  Initial  stages  of  this  investigation  the 
objective  was  to  determine  the  influence  of  recirculation 
on  combustion  chamber  performance.  Consequently,  the 
chamber  desired  was  one  in  which  recirculation  played  an 
important  role,  and  in  which  some  control  could  be 
exercised  over  the  recirculation.  At  the  same  time,  care 
was  necessary  to  insure  that  the  combustor  chosen  be 
representative  of  conventional  chamber  design,  since  it 
was  felt  desirable  that  the  experimental  program  should 
avoid  the  necessity  of  studying  phenomena  which  were  the 
results  of  a peculiarity  of  a particular  chamber  rather 
than  principles  of  general  applicability  to  chamber  design. 

Several  alternative  chambers  were  considered.  Since 
recirculation  was  the  dominant  characteristic  desired,  a 
chamber  in  which  the  flow  was  circulatory  was  immediately 
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suggested.  Three  such  chamber  designs  have  appeared  in 
the  recent  literature.  Schmidt  and  Schoppe  (19)  presented 
two  chambers,  both  for  non-premixed  systems  in  which  gases 
flowed  in  a spiral  manner.  In  these  chambers  stabiliza- 
tion was  accomplished  by  the  recirculation  of  burning 
gases  to  the  cold  zones  near  the  fuel  and  air  inlet 
points.  Very  little  data  was  given  on  the  behavior  of 
these  chambers  except  that  the  authors  claimed  wide 
stability  limits  and  high  heat  liberation  rates,  roughly 
i comparable  to  can-type  combustors  (l  to  10  x 10® 

Btu/hr  ft3  at  1 atm  pressure).  Another  combustor  of  this 
general  type  was  suggested  by  Hottel  and  Person  (9).  They 
studied  a vortex  type  system  similar  to  the  one  suggested 
by  Schoppe,  but  with  a different  fuel  injection  arrange- 
ment; however,  as  did  Schoppe's,  this  chamber  used  a 
non-premixed  fuel  and  air  feed. 

None  of  these  chamber  designs  seemed  particularly 
adapted  to  the  present  study.  Person  gave  moderately 
detailed  information  concerning  the  flow  pattern,  but  he 
found  that  there  was  little  recirculation  in  his  chamber. 

Schmidt  and  Schoppe  gave  no  Information  at  all  about  flow 
patterns  in  their  chambers  nor  were  any  details  about  the 
performance  of  their  chambers  presented.  Other  designs 
Involving  forced  cyclonic  flow  were  also  considered. 

These  various  alternatives,  however,  were  all  considered 
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to  be  so  non-typical  of  a conventional  chamber  as  not  to 
merit  detailed  study. 

A simple,  rectangular  two-dimensional  chamber  fed  by 
a Jet  of  a pre-mixed  fuel/air  mixture  was  then  considered. 
This  had  the  advantage  of  being  somewhat  more  analogous  to 
a practical  chamber,  yet  simple  enough  in  design  to  permit 
easy  study  of  its  behavior.  In  such  a chamber,  recircu- 
lation is  achieved  by  the  aspirating  effect  of  the  Jet. 
hear  the  chamber  exit,  material  equivalent  to  the  material 
aspirated  must  reverse  direction  and  form  recirculation 
zones  around  the  central  core  of  the  Jet.  Control  over 
the  amount  of  recirculation  can  be  achieved  by  varying 
the  ratio  of  chamber  width  to  nozzle  with,  ana  wltn  quite 
different  results,  by  varying  the  chamber  dimension  ratios. 

In  the  initial  stages  of  this  problem  there  was 
concern  over  the  possible  effects  of  heat  loss  on  the 
experimental  measurements.  To  minimize  wall  heat  losses 
the  chamber  was  built  with  three  parallel  feed  slots  in 
its  floor.  It  was  hoped  that  the  experiments  could  be 
carried  out  with  no  significant  Interaction  of  the  Jets; 
however,  Jet  interactions  proved  to  be  so  large  that  it 
was  necessary  to  return  to  a single- Jet  design  to  study 
recirculation  in  this  system. 

It  is  interesting  to  note  the  similarity  of  this 
chamber  to  the  well-stirred  reactor.  If  the  volume 
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associated  with  a single  jet  of  the  stirred  reactor  is 
compared  with  the  jet  of  a chamber  into  which  a single 
Jet  flows,  similarity  of  the  two  systems  becomes  apparent. 
One  may  then  hope  that  a slngle-Jet  system  will  provide 
an  insight  into  the  operation  of  the  stirred  reactor. 

The  design  of  the  chamber  was  based  on  this 
similarity  to  the  stirred  reactor.  Three  parallel  slots, 
each  having  an  opening  of  0.030"  x 2"  were  placed  along 
the  floor  of  a rectangular  combustion  chamber  which  was 
3 1/2"  x 3"  x 2"  in  Internal  dimensions.  (The  floor  of 
the  chamber  is  represented  by  the  3"  dimension.)  Along 
the  roof  of  the  chamber  were  three  opposing  exit  flues 
1/2"  x 2"  in  size.  The  combustion  chamber  was  formed  by 
placing  insulating  firebrick  inside  of  the  frame  which 
held  the  nozzles.  The  inside  chamber  dimensions  would 
be  readily  altered  by  varying  the  size  of  the  firebrick 
used  as  lining,  and  the  nozzle  spacing  was  controlled  by 
an  assortment  of  Interchangeable  spacer  blocks  separating 
the  nozzle  bodies.  The  front  and  back  of  the  chamber  were 
formed  by  Vycor  windows  held  by  steel  frames  and  a sheet- 
metal  retaining  piece.  The  entire  chamoer  was  placed  Inside 
a vacuum  tank  to  control  the  pressure  of  operation.  A 
sketch  of  the  combustion  chamber  assembly  is  shown  in 
Figure  IV-1,  and  photographs  of  the  apparatus  are  shown 
in  Figures  IV -2  to  IV -4 . 


Pre-mixed  fuel/air  mixtures  were  fed  to  this  chamber. 
Fuel  and  air  were  metered  separately,  mixed  in  a length  of 
pipe  containing  baffles,  and  fed  to  a manifold.  To  Insure 
that  each  Jet  was  operating  at  the  same  feed  rate,  the 
manifold  mixture  was  metered  to  each  Jet.  Standard 
square-edges  orifices  were  used  for  all  metering.  A flow 
diagram  is  shown  in  Figure  IV-5,  including  the  metering 
and  temperature  measurement  points.  Photographs  of  the 
control  panel  are  shown  in  Figures  IV-6  and  IV-7. 

After  the  experimental  program  was  started  it  became 
evident  that  a one-jet  chamber  should  be  studied  as  well 
as  the  three- Jet  chamber.  This  was  done  in  the  apparatus 
described  above  by  replacing  the  side  bricks  so  that  a 
narrower  interior  width  was  obtained  and  shutting  off  the 
feed  to  the  outer  two  Jets.  In  the  original  one- jet 
chamber  the  Inside  width  was  made  Just  1/3  of  the  three- 
Jet  chamber  width,  i.e.  1".  This  proved  to  be  too  narrow 
to  achieve  flame  stabilization  at  reasonable  feed  rates. 

The  inside  chamber  dimension  was  therefore  enlarged  to  a 
2"  width  (by  2"  deep  and  3 1/2"  high)  for  the  study  of  the 
one-jet  chamber.  This  chamber  is  shown  in  the  photograph 
IV-8  (also  shown  is  a mounted  sampling  probe  used  to 
sample  the  Interior  gases) . 

This  chamber  was  placed  inside  a small  vacuum  tank 
so  that  it  could  be  operated  at  pressures  ranging  from 


QUENCH  AIR 


FIGURE  EZ-5:  FLOW  DIAGRAM  OF  APPARATUS 


FIGURE  UC-7:  PHOTOGRAPH  OF  CONTROL  PANEL  SHOWING  SAMPLING 
BOARO  ASSEMBLY. 
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about  0.1  atm  to  0.5  atm.  However,  except  for  one  series 
of  runs  made  to  determine  the  effect  of  pressure  on  the 
blowout  feed  rate,  the  data  were  all  taken  at  a pressure 
of  about  0.34  atm,  or  10"  Hg  absolute. 

Sampling  Equipment 

Some  additional  equipment  was  required  for  sampling 
the  combustion  gases.  This  consisted  of  a water-cooled 
sampling  probe,  sampling  rate  measurement  and  control 
equipment,  and  analytical  equipment.  For  sampling  gases 
Inside  the  chamber  probes  of  the  same  general  design  used 
by  Baker  (3)  were  used.  This  design  consisted  of  three 
concentric  hypodermic  tubes  assembled  as  shown  in  Figure 
IV-9,  forming  a central  tube  and  two  concentric  annuli. 

High  pressure  water  passed  at  a rate  of  10  to  50  cc/min 
down  the  inner  annulus  and  returned  to  drain  in  the  outer. 

About  200  to  400  psi  were  required  to  send  this  amount  of 
water  through  the  probe  assembly.  Sample  was  withdrawn 
through  the  Innermost  tube. 

A somewhat  different  design  was  used  for  the  probe 
which  sampled  the  exit  gases.  For  this  a single,  long 
tube  entered  through  one  side  of  the  supporting  framework, 
crossed  the  flue  in  the  plane  of  the  chamber  roof,  and 
left  through  the  opposite  side  of  the  framework.  A sample 
inlet  was  drilled  into  the  side  of  this  long  cooling  water 
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FIGURE  U-9:  REVERSE-FLOW  COOLING  WATER  SAMPLING  PROBE 
(10  ft  ACTUAL  SIZE) 


COOLING  WATER  IN 


FIGURES- 10:  STRAIGHT-THROUGH  FLOW  COOLING  WATER 
SAMPLE  PROBE  (10  X ACTUAL  SIZE) 
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tube  and  connected  to  a concentric  inner  tube.  A sketch 
of  this  probe  Is  shown  In  Figure  IV-10.  Photographs  of 
the  assembled  probes  used  are  shown  In  Figure  IV-11. 

The  analysis  of  the  samples  in  this  work  was 
accomplished  In  one  of  two  ways:  (1)  For  those  samples 
where  oxygen  consumption  efficiency  was  desired,  the  analysis 
of  the  gases  was  made  with  a Beckman  Oxygen  Meter,  operating 
on  the  para-magnetic  property  of  oxygen.  (2)  For  those 
samples  where  a more  complete  analysis  was  desired,  the 
analysis  was  performed  by  standard  Orsat  techniques. 

To  sample  properly  gases  flowing  at  high  velocity  and 
varying  In  composition  it  Is  desirable  to  match  the  stream 
velocities  Just  inside  and  outside  of  the  tube.  To 
accomplish  this  the  impact  pressure  at  the  tip  of  the 
probe  must  first  be  measured,  and  the  sample  then  with- 
drawn through  the  probe  at  that  rate  calculated  to  match 
the  velocities.  This  Is  a trial  and  error  procedure  since 
the  gas  temperature  must  be  calculated  from  the  composition 
In  order  to  find  the  desired  sampling  rate.  The  sampling 
train  required  must  include:  (l)  a manometer  to  measure 
the  Impact  pressure,  ( 2 ) flow  measurement  and  control 
equipment,  and  (3)  either  an  oxygen  meter  or  a sample 
collection  bottle.  A diagram  of  the  sampling  train  used 
Is  shown  In  Figure  IV-12. 
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FIGURE  CC  - 12:  FLOW  DIAGRAM  OF 
SAMPLING  TRAIN 
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Experimental  Procedures 

Three  types  of  experiments  were  made  on  this  apparatus: 

(1)  stability  limits  at  constant  pressure,  giving  fuel* 
air  ratio  as  a function  of  the  blowout  air  feed  rate; 

(2)  efficiency  measurements,  giving  oxygen  efficiency 
vs  feed  rate  for  a constant  fuel/air  ratio  and  pressure; 
and  (3)  pressure  exponent  measurement,  giving  pressure  vs 
feed  rate  at  blowout  for  a constant  fuel/air  ratio.  The 
procedures  used  In  these  experiments  are  described  in 
greater  detail  below. 

Blowout  Measurements  at  Constant  Pressure 

The  constant-pressure  stability  limits  were  determined 
for  both  the  one-jet  and  the  three-jet  combustion  chambers. 
The  flame  was  lit  by  a high-voltage  spark  above  the  exit 
flue  (or  flues)  at  a very  low  feed  rate  and  high  fuel/air 
ratio.  If  the  fuel/air  ratio  was  then  carefully  reduced  the 
flame  would  flash  back  into  the  combustor.  The  fuel  and 
air  rates  were  then  slowly  raised  until  the  desired 
operating  air  rate  was  reached.  Simultaneously  the 
brickwork  became  a bright  cherry  red.  In  the  case  of  the 
three-jet  chamber,  flash-back  was  more  easily  accomplished 
by  feeding  one  of  the  outside  Jets,  and  then,  after  that 
Jet  had  been  lit,  gradually  introducing  the  fuel/air 
mixture  feed  to  the  other  two  Jets.  After  the  chamber 
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had  reached  a steady-state  operation  at  the  desired 
feed  rate,  the  desired  blow-out  measurements  could 

be  made.  With  pressure  and  air  feed  rate  constant, 

( 

the  fuel-air  ratio  was  slowly  changed  until  blowout 
was  reached.  As  the  stability  limit  was  approached, 
a few  minutes  were  allowed  after  each  small  change 
in  fuel  rate  so  that  the  chamber  operation  could 
reach  thermal  equilibrium.  (The  use  of  kaolin  brick 
minimized  the  necessary  waiting  time.) 

In  the  three-jet  chamber  the  three  Jets  coalesced 
a short  distance  downstream  of  the  nozzles  so  that  only 
two  major  flame  areas  were  observed.  At  sufficiently 
low  flow  rates,  two  small  zones  trapped  between  pairs 
of  the  three  Jets  near  the  base  of  the  chamber  would 
also  light,  but  in  the  regions  studied  they  did  not 
contribute  to  the  chamber  stability,  and  in  general, 
gases  did  not  burn  stably  in  them.  In  this  chamber, 
as  the  fuel/air  ratio  was  changed  two  blowout  limits 
were  observed,  one  in  which  one  of  the  two  active 
(large  recirculation)  zones  blew  out,  with  a consequent 
Increase  in  the  volume  of  remaining  active  zone,  and 
the  second  when  the  flame  in  the  remaining  active  zone 
was  no  longer  stable. 

In  the  one- Jet  chamber  the  phenomenon  of  the 
two  stability  limits  was  not  observed  because  the 
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flame  would  not  burn  stably  in  a symmetrical  position. 
The  only  stable  combustion  pattern  was  the  one  in  which 
the  jet  was  bent  over  against  one  wall  with  combustion 
visible  only  in  the  enlarged  recirculation  zone.  This 
undoubtedly  resulted  from  the  narrowness  of  the 
chambers,  as  in  a 3B  wide  chamber  (as  compared  with 
the  2"  chamber)  a symmetrical  combustion  pattern  could 
be  obtained.  A second  interesting  characteristic  of 
the  blowout  in  the  2",  one-jet  chamber  was  that  a 
residual  flame  was  consistently  observed.  The  combus- 
tion pattern  which  developed  in  this  chamber  had  a 
circulatory  pattern  in  which  the  mainstream  of  the  Jet 
went  up  along  one  wall  and  recirculating  gases  returned 
along  the  other.  A small  recirculation  zone  developed 
in  the  center  of  the  vortex  created  by  this  motion. 

At  low  feed  rates  combustion  proceeded  primarily  at 
the  edge  of  the  mainstream  of  the  Jet  and  the  edge 
of  the  recirculating  gases.  The  central  portion  of 
the  recirculation  zone  was  nearly  Inactive  since  the 
gases  reaching  it  were  almost  completely  burned  (see 
Figure  V-7).  As  higher  rates  were  reached,  approaching 
the  blowout  limit,  there  was  a gradual  transition  in 
which  the  active  combustion  region  shifted  to  the 
central  regions,  forming  a residual  flame  near  blowout 
(Figure  V-8  shows  a stage  in  this  transition).  The 
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data  reported  for  this  chamber  are  the  gross  chamber 
loading  rates  necessary  to  blow  out  this  residual 
flame.  Constant-pressure  blowout  data  were  taken  in 
the  same  manner  as  that  described  for  the  three- Jet 
chamber. 

For  these  chambers  blowout  loading  rates  were 
measured  as  a function  of  the  fuel-air  ratio.  data 
were  taken  at  an  approximately  constant  absolute 
pressure  of  10"  mercury.  Several  air  rates  were  used 
for  each  chamber,  and  at  each  air  rate  the  fuel  rate 
was  varied  to  obtain  both  lean  and  rich  blowout  limits. 
An  ambient  inlet  temperature  of  approximately  540  °H 
was  used,  deviations  in  both  pressure  and  temperature 
from  the  desired  values  were  corrected  for  by  well- 
stirred  reactor  theory.  The  temperature  correction 
charts  of  Longwell  and  Weiss  (14)  were  used  to  correct 
the  experimental  data  to  540 °H,  and  a pressure  exponent 
of  1.8  was  used  to  correct  for  deviations  in  pressure. 
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The  observed  change  in  visible  flame  pattern 
with  changes  in  feed  rate  to  the  2",  one-jet  chamber 
suggested  that  its  variation  of  efficiency  vs  loading 
rates  might  conform  approximately  to  the  calculated 
behavior  of  the  recirculation  model.  Hence  this 
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chamber  was  chosen  for  a more  detailed  study  of  Its 
performance.  A series  of  runs  were  made  to  determine 
the  relation  between  combustion  efficiency  and  feed 
rate  and  another  group  was  made  to  determine  the  effect 
of  pressure  on  the  blowout  feed  rate  at  several  fuel- 
air  ratios. 

The  sampling  of  this  chamber  was  accomplished  by 
means  of  the  gas  sampling  probes  described  above.  In 
these  runs  a fixed  fuelralr  ratio  of  1.35  times  stoi- 
chiometric was  used.  Various  air  feed  rates  were  used 
at  a constant  pressure  to  determine  a curve  of  loading 
rate  vs  the  oxygen  consumption  efficiency.  Analysis 
was  done  by  means  of  a paramagnetic  oxygen  analyzer. 

A traverse  of  oxygen  efficiency  vs  distance  Into  the 
chamber  at  one  level  was  made  for  various  air  rates 
to  locate  the  position  and  size  of  the  combustion 
zones  as  well  as  to  determine  the  combustion  efficiency 
of  the  gases  within  these  zones.  The  position  of  the 
probe  tip  In  these  traversing  measurements  was  deter- 
mined by  use  of  a cathetometer . 

Buns  were  also  made  to  determine  the  gross 
efficiency  of  the  combustion  at  the  chamber  outlet. 

In  these  runs  a probe  of  the  straight-through  design 
described  In  Figure  IV-10  was  used,  burled  In  the  roof 
of  the  chamber  so  that  the  sampling  opening  waB  In  the 
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plane  of  the  chamber  roof  and  pointed  toward  the  on- 
coming gases.  In  sampling  this  gas  stream,  a case 
where  the  gases  Impinged  on  the  sample  opening,  an 
effort  was  made  to  match  velocities.  The  stream 
velocity  was  estimated  by  making  an  Impact  pressure 
measurement  through  the  sample  probe  before  the  sample 
was  withdrawn.  Since  the  gas  temperature  was  assumed 
to  be  dependent  only  the  gas  composition  as  read  from 
the  oxygen  meter,  the  required  trial-and-error  procedure 
was  relatively  simple. 

The  measurement  of  a satisfactory  static  pressure 
against  which  to  compare  the  Impact  pressure  presented 
a problem.  Purge  air  was  introduced  Into  the  vacuum 
tank  to  sweep  the  tank  clean  of  combustion  products, 
and  to  keep  the  tank  Internals  cool.  The  flow  of  this 
purge  air  created  a pressure  difference  of  several 
tenths  of  an  Inch  of  water  between  the  static  pressure 
measured  at  the  tank  wall  and  the  Impact  pressure  open- 
ing at  the  exit  flue.  An  Improvement  In  this  measure- 
ment was  achieved  by  an  extension  tube  Inside  of  the 
tank  wall  to  a point  adjacent  to  the  combustion  chamber 
frame.  Although  this  poor  pressure  measurement  may 
have  caused  some  error  In  matching  the  velocities,  the 
effect  of  such  a mismatch  on  the  measured  oxygen  effi- 
ciency was  found  experimentally  to  be  small.  Comparison 
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of  several  sampling  velocities  for  a given  chamber 
condition  showed  that  mismatching  ranging  from  50^ 
to  120^  of  the  stream  velocity  caused  a variation  of 
less  than  in  the  oxygen  efficiency  measured. 

In  measuring  the  gross  efficiency,  attention  had 
to  be  given  to  the  variation  of  combustion  efficiency 
across  the  flue  opening  even  though  it  was  only  1/2" 
wide.  Since  the  combustion  pattern  in  the  chamber  was 
highly  assymetric,  the  combustion  efficiency  on  one 
side  of  the  flue  reflected  the  behavior  of  the  cold 
Jet  mainstream,  while  that  on  the  other  side  of  the 
flue  reflected  behavior  in  the  more  completely  burned 
regions  of  the  chamber.  To  determine  the  gross  com- 
bustion efficiency,  therefore,  gas  samples  were  taken 
at  each  edge  of  the  flue  and  at  the  center.  The 
position  of  the  probe  opening  was  determined  visually 
in  these  runs;  this  was  considered  adequate  because 
the  variation  found  in  gas  composition  was  small. 
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Effect  of  Pressure  on  Blowout 

The  effect  of  pressure  on  the  blowout  loading 
was  also  studied  at  four  different  fuel/alr  ratios. 

In  these  runs  several  feed  rates  were  established  at 
constant  fuel/alr  ratio,  and  after  each  air  rate  was 
set,  the  pressure  in  the  vacuum  tank  was  slowly  reduced 
until  blowout  was  reached.  After  each  blowout  this 
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procedure  was  repeated  for  the  next  air  rate.  The  fuelr 
air  ratios  studied  corresponded  to  fuel  fractions,  F, 
of  0.600  and  0.575  on  the  rich  side  and  0.450  and  0.475 
on  the  lean  side. 

Combustion  Efficiency  Measurements  - Three-Jet  Chamber 

In  addition  to  the  above  data  on  the  one-jet 
chamber  a few  sampling  data  were  obtained  on  the  ori- 
ginal three- Jet  chamber.  These  were  taken  near  the 
rich  blowout  limit  at  a constant  air  rate  and  varying 
fuel-air  ratios.  A probe  of  the  reversed-flow  cooling 
water  design  was  built  so  that  it  could  be  dropped  down 
from  the  top  through  each  of  the  three  exit  flues. 

Probe  movement  was  accomplished  by  means  of  an  air 
solenoid  and  an  electric  actuator  motor.  The  analysis 
was  by  means  of  the  oxygen  meter.  These  data  were  taken 
to  establish  that  the  three-jet  chamber  was  operating 
in  a reasonable  range  of  fuel/air  ratio  rather  than  to 
establish  the  efficiency-loading  rate  curve  for  the 
chamber.  Consequently  the  data  are  very  limited. 
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CHAPTER  V 

EXPERIMENTAL  RESULTS 

Three-Jet  Chamber  Operation 

The  three-jet  combustion  chamber  was  the  first  of 
the  chambers  to  be  investigated.  Lata  taken  on  this 
chamber  were  primarily  blowout  data.  Some  photographic 
measurements  of  the  combustion  pattern  and  a few 
efficiency  measurements  were  also  made. 

The  flow  pattern  in  this  chamber  showed  a strong 
interaction  of  the  jets.  Convergence  of  the  two  out- 
side jets  was  observed  both  with  and  without  combustion 
occurring  in  the  chamber.  The  three  jets  formed  four 
recirculation  zones  whose  pattern  was  symmetrical  around 
the  center  plane  of  the  chamber.  The  two  inside  recir- 
culation zones  trapped  by  pairs  of  the  jet  mainstreams 
were  relatively  small.  The  outer  two  occupied  the  space 
made  available  by  the  convergence  of  the  three  Jets. 

The  outer  zones  behaved  essentially  as  recirculation 
zones  to  the  large  central  Jet  created  by  the  convergence 
of  the  three  Jets. 

When  combustion  in  the  chamber  was  first  established 


at  very  low  flow  rates,  the  flame  occupied  four  zones, 
corresponding  to  the  four  recirculation  zones  of  the 
cold  flow  pattern.  There  was  no  visible  evidence  in 
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this  chamber  of  combustion  occurring  in  the  mainstreams 
of  the  jets.  A photograph  of  the  chamber  operating  under 
these  conditions  is  shown  in  Figure  V-l.  It  will  be 
noted  that  the  combustion  pattern  is  slightly  assymetrlc, 
indicating  a possible  mis-matching  of  nozzle  sizes  or  a 
mis-alignment  of  the  nozzle  axes.  This  assymetry  is 
discussed  in  detail  in  the  next  chapter  (Errors  in  the 
Experimental  Work). 

As  the  feed  rate  was  increased  the  two  smaller 
zones  became  unstable,  and  combustion  continued  only 
in  the  two  larger  outer  zones,  Figure  V-2.  The  final 
stable  combustion  pattern  occurred  when  one  of  these  two 
larger  recirculation  zones  "blew  out",  and  the  remaining 
active  zone  expanded  to  a more  stable  position,  Figure  V-3. 
Continued  increase  in  the  feed  rate  caused  blow-out  of 
this  final  flame. 

Blowout  measurements  on  this  chamber  were  made  to 
determine  the  loading  rate  at  which  the  shift  from  the 
two-actlve-zones  to  the  one-active-zone  pattern  occurred, 
and  to  determine  the  gross  loading  factor  at  the 
stability  limit  of  the  one-active-zone  pattern.  Wo 
measurements  were  made  on  the  stability  of  the  two 
smaller  zones  since  these  did  not  contribute  materially 
to  the  operation  of  the  chamber.  Figure  V-4  shows 
these  results.  To  facilitate  comparison  with  the  data 
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FIGURE  V-4:  BLOWOUT  DATA  OF  THREE-JET  CHAMBER 

COMPARED  TO  STIRRED  REACTOR  (14) 


on  the  well-stirred  reactor,  uhe  data  are  plotted  on 
coordinates  of  N /VP1,8  vs  generalized  fuel  fraction, 

fi 

P,  along  with  the  corresponding  well-stirred  reactor 
data.  These  runs  were  all  made  at  an  absolute  pressure 
of  approximately  ten  inches  of  mercury,  and  an  inlet 
temperature  of  540 °H  (ambient) . No  temperature  control 
was  used,  but  the  data  were  corrected  for  the  small 
variations  that  occurred  in  inlet  temperature  by  using 
the  temperature-correction  charts  of  Longwell  and  Weiss 
(14) . 

Another  series  of  runs  using  this  three-jet 
chamber  was  made  to  estimate  the  size  of  the  various 
zones.  Measurements  of  zone  size  were  made  from 
photographs.  This  technique  for  the  estimation  of  the 
volume  of  combustion  zones  has  been  used  by  Simon  and 
Wagner  (22).  Since  these  photographic  measurements 
were  to  be  used  in  conjunction  with  the  blowout  data, 
the  photographs  were  taken  with  the  chamber  feed  rate 
near  blowout  for  the  range  of  fuel/air  ratios  of 
interest.  Figure  V-5  gives  the  data  on  the  estimated 
flame  volumes.  A typical  photograph  showing  the  outlines 
of  the  flame  zone  which  were  planlmetered  to  determine 
volumes  is  also  shown  in  Figure  V-5. 

Some  measurements  of  combustion  efficiency  in  this 
chamber  were  also  made.  A gas-sampling  probe  was 
inserted  between  the  flue  bricks  to  the  level  of  the 
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Inside  chamber  roof.  Samples  were  withdrawn  through  this 
probe  and  analyzed  for  oxygen  content  using  a paramagnetic 
oxygen  analyzer.  The  few  data  which  were  taken  are 
shown  In  Figure  V-6. 

Single-Jet  Chamber  Operation 

In  the  2"  wide,  single-jet  chamber  no  symmetrical 
combustion  pattern  was  observed.  The  axis  of  the  Jet 
curved  over  against  one  side  of  the  chamber,  with  com- 
bustion occurring  in  the  large  recirculation  zone  thus 
created.  A photograph  of  a typical  flame  pattern  at 
moderate  firing  rates  is  shown  in  Figure  V-7.  Flame  in 
this  chamber  appeared  in  the  boundaries  of  the  recircu- 
lation zone,  where  relatively  unburned  gases  were  avail- 
able. But  in  the  "eye"  of  the  zone  no  unburned  gases 
could  penetrate,  and  consequently  no  flame  appeared  in 
this  volume.  There  was  no  visible  evidence  of  combustion 
occurring  in  the  mainstream  of  the  Jet,  presumably  due 
to  the  high  velocity  and  consequent  short  residence  time 
of  this  region. 

As  the  feed  rate  to  this  chamber  increased,  the 
zone  of  active  combustion  moved  toward  the  center, 
occupying  more  of  the  "eye"  and  receding  from  the  wall, 
until  at  rates  near  the  ultimate  blowout  limit,  all 
that  remained  of  the  active  zone  was  a small  ball  of 
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FIGURE  2-10:  OVERALL  COMBUSTION  EFFICIENCY  OF  ONE- JET  CHAMBER 
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flame  in  the  place  where  the  dark  "eye"  of  the  low-flow- 
rate pattern  had  been.  A photograph  of  the  residual 
flame  partially  formed  (but  still  not  at  the  blowout 
limit)  is  shown  in  Figure  V-8. 

This  one-jet  chamber  was  the  most  extensively 
studied.  Blowout  limits  were  the  first  characteristic 
Investigated.  Since  the  only  unambiguous  blowout  was 
the  disappearance  of  the  residual  flame,  only  this 
stability  limit  is  reported  here.  These  blowout  data 
are  shown  in  Figure  V-9.  As  in  the  case  of  the  stability 
data  for  the  three- Jet  chamber  these  data  were  taken 
at  a pressure  of  10"  of  mercury  absolute  and  a tempera- 
ture of  540°H.  They  are  plotted  on  coordinates  of 
N^/VP1,8  vs  F together  with  the  Longwell  and  Weiss  data. 

It  was  observed  that  the  behavior  of  this  chamber 
qualitatively  followed  that  which  would  be  expected  from 
the  recirculation  model  discussed  in  Chapter  III.  To 
investigate  this  further,  combustion  efficiency  was 
determined  as  a function  of  loading  rate  at  a constant 
fuel-air  ratio  and  pressure.  These  measurements  were 
made  by  means  of  a horizontal  probe,  with  a side-entrance 
sampling  port,  burled  in  the  roof  of  the  chamber  and 
extending  across  the  flue.  Because  of  the  assymetry 
of  the  combustion  pattern,  three  samples  were  taken  at 
each  air  rate,  one  from  each  flue  edge  and  one  from 
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the  center  of  the  flue.  The  results  are  plotted  In 
Figure  V-10.  Since  this  plot  shows  that  there  Is  only 
a small  difference  among  these  three  samples,  the 
simple  arithmetic  mean  will  be  used  In  the  analysis  of 
the  data. 

To  compare  the  chamber  performance  with  the  theo- 
retical model,  quantitative  data  on  the  composition 
and  sizes  of  the  various  zones  are  also  desirable.  These 
were  obtained  by  a combination  of  sampling  In  the  chamber 
Interior  and  photographic  measurement.  Samples  of  the 
chamber  contents  were  obtained  by  a sampling  traverse 
across  the  chamber.  This  was  done  for  several  feed 
rates,  but  at  only  one  level  in  the  chamber,  about  2" 
from  the  floor,  and  one  fuel/alr  ratio.  One  unfortunate 
result  of  this  procedure  was  that  the  flow  pattern  in 
the  chamber  was  upset  by  the  Introduction  of  the  probe. 
This  was  indicated  by  a reduction  in  the  stability 
limits  for  the  chamber  when  sampling  measurements  were 
being  made.  This  reduction  was  a function  of  how  far 
the  probe  had  been  Inserted  Into  the  chamber.  These 
data  are  shown  In  Figure  V-ll. 

The  existence  of  the  central  eye  and  outer  ring 
zones  Is  Indicated,  and  an  estimate  of  zone  volumes  can 
be  made  by  considering  the  midpoint  In  efficiency  between 
the  maximum  and  wall  values  to  be  the  zone  boundary.. 
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However,  the  zone  concept  Is  better  Illustrated  by  cross- 
plotting Figure  V-ll  to  give  efficiency  versus  loading 
rate  for  each  of  several  positions  In  the  chamber 
(Figure  V-12).  Examination  of  that  figure  shows  that 

t 

the  2 and  3 cm.  positions  are  characterized  by  a quite 
different  response  to  loading  from  that  of  the  0,  1, 
and  4 cm.  positions.  The  5 cm.  position  (the  cold 
wall  of  the  chamber) , where  the  mainstream  of  the  non- 

’ 

burning  Jet  passes  shows  some  combustion  because  of 

i 

mixing  with  the  adjacent  zone. 

Photographs  of  the  flame  for  two  loading  rates 
with  the  probe  withdrawn  are  shown  in  Figure  V-13.  The 
probe  positions  in  the  sampling  are  Indicated  on  these 
photographs . 

i Finally  the  effect  of  pressure  on  the  blowout 

limit  was  studied.  Buns  were  made  by  establishing  a 
stable  flame  at  a predetermined  value  of  feed  rate  and 
fuel/alr  ratio.  The  pressure  on  the  chamber  was  then 
gradually  reduced  until  blowout  occurred.  During  this 
series  of  runs  the  pattern  of  the  approach  to  blowout 
was  not  noticeably  affected  by  the  pressure  level. 

Data  showing  the  effect  of  pressure  on  blowout  are  given 
in  Figure  V-14.  Four  fuel-air  ratios  were  studied.  A 
plot  of  pressure  exponent  vs  the  generalized  fuel  frac- 
tion is  shown  In  Figure  V-15.  Also  shown  In  Figure  V-14 
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FIGURE  31-15:  EFFECT  OF  FEED 

COMPOSITION  ON  PRESSURE  EXPONENT 


CHAPTER  VI 


DISCUSSION  OP  RESULTS 


The  experimental  program  was  undertaken  to  provide 
data  for  examining  the  question  of  whether  stirred 
reactor  performance  data  could  be  combined  with  fluid 
flow  patterns  to  predict  the  performance  of  practical 
combustors.  Using  stirred  reactor  blowout  data,  model 
performance  for  the  two-stability-limit  model  may  be 
predicted.  Since  sufficient  experimental  data  on  the 
efficiency-loading  relation  for  the  stirred  reactor 
are  not  available,  performance  of  the  recirculation 
model  must  be  predicted  from  theoretical  stirred-reactor 
performance.  The  analysis  of  the  data  consists  princi- 
pally of  direct  comparison  of  these  theoretical  models 
with  the  experimental  reactor  results. 

Experimental  Errors 

The  discussion  of  the  experimental  work  must  be 
prefaced  by  commentary  on  the  errors  involved  in  these 
experiments.  Two  significant  sources  of  experimental 
error  were  present;  air  leakage  into  the  combustion 
zone  and  non-ideality  of  the  stirred  reactor.  A third 
lay  in  the  construction  of  the  nozzle  header  assembly 
which  ideally  should  have  perfectly  aligned  and  matched 
nozzles.  This  last  was  not  as  much  a problem  as  the 
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first  two  because  the  effect  of  mis -matched  or  mis- 
aligned nozzles  was  easily  detected.  Each  of  these 
problems  will  be  discussed  In  turn. 

Air  Leakage.  The  combustion  chamber,  consisting 
as  described  previously  of  an  arrangement  of  firebrick 
and  nozzles  between  two  frames  carrying  the  Vycor  windows, 
Is  fed  by  high  velocity  Jets  along  Its  floor.  In  any 
system  of  this  sort  the  Jets  aspirate  gases  from  the 
surrounding  fluid,  and  near  the  top  of  the  chamber  a 
reverse  flow  of  the  chamber  a reverse  flow  occurs  to 
replace  the  gases  aspirated  from  the  recirculation  zone. 

To  create  this  reverse-flow  pattern  requires  (from  a 
momentum  balance)  a low-pressure  region  In  the  lower 
portion  of  the  chamber.  If  there  Is  an  Improper  seal 
between  the  chamber  and  its  side  walls,  substantial 
quantities  of  air  can  be  aspirated  from  the  ambient  air 
Into  this  low  pressure  region  of  the  combustlbn  chamber. 

In  the  construction  of  the  apparatus,  precautions  were 
taken  to  seal  the  chamber  tightly.  However,  from  time 
to  time  thermal  warpage  of  the  side  frames  or  Improperly 
fitted  bricks  made  a satisfactory  seal  impossible.  This 
problem  may  be  alleviated  to  some  extent  by  careful 
sizing  of  the  exit  flue  or  flues.  If  these  are  made 
narrow  enough  the  pressure  drop  required  to  force  the 
exhausted  gases  through  them  may  be  sufficient  to  raise 
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the  pressure  at  the  top  of  the  chamber  so  that  the  low- 
pressure  region  in  the  recirculation  zone  is  at  or 
above  ambient  pressure. 

The  presence  of  an  air  leak  in  the  chamber 
assembly  could  be  most  easily  detected  by  a shifting 
of  the  stability  limits  of  the  chamber  toward  the  rich 
side  (l.e.  stability  data  which  lacked  symmetry  around 
a stoichiometric  fuel/air  ratio) . This  criterion, 
however,  was  not  highly  sensitive,  as  a substantial 
shift  was  required  to  differentiate  assymetry  due  to 
air  leakage  from  the  small  assymetry  inherent  in  the 
well-stirred  reactor  stability  curves. 

To  test  whether  the  chamber  was  properly  sealed 
when  the  blowout  curves  were  symmetrical,  one  run  was 
made  in  which  a sample  of  gas  was  withdrawn  from  the 
Interior  of  the  chamber.  An  Orsat  analysis  of  this 
sample  for  the  carbon/ nitrogen  ratio  was  compared  to 
the  ratio  caTculated  from  the  metered  rates  of  fuel  and 
air.  The  results  are  tabulated  in  Table  VI-1,  Included 
in  this  table  is  the  hydrogen/carbon  ratio  which  serves 
as  a check  by  comparison  with  the  hydrogen/carbon  ratio 
for  propane.  Since  only  enough  data  were  obtained  to 
calculate  the  H/C  ratio  from  the  known  ratio  of 
nitrogen/oxygen  in  air  (or  visa-versa) , only  one  check 
on  the  analysis  could  be  calculated.  Indicated  in 


Table  VI-1  Is  the  H/C  ratio  of  propane,  2.67.  Analyses 
of  the  fuel  have  shown  the  propane  to  be  about  96^ 
pure  (10) , the  principal  impurities  being  ethane, 
propylene,  and  butane.  Assuming  a typical  analysis  as 

- 

given  by  Woo  (30)  of  0.8^  ethane,  1.5^  propylene,  1.5 % 

a 

butane,  and  the  balance  propane,  the  calculated  atoms 
of  carbon  per  mole  of  propane  are  3.01  and  the  H/C  ratio 
is  2.66. 

TABLE  VI-1 

[ i 

Comparison  of  Carbon/Hi trogen  ratio  Determined  by 
Sampling  with  the  Metered  Hatio 

Value  from  Sample  Calculated  by  Metering 
C/N  0.150  0.156 

1H/C  2.68  2.56 

When  the  sampling  data  reported  in  this  table  were 
taken  the  stability  limits  of  the  combustor  were  ap- 
proximately symmetrical.  In  this  table  it  is  seen  from 
the  comparison  of  H/C  that  the  analysis  of  the  sample 
is  reliable.  The  agreement  between  the  C/N  ratios, 
though  it  indicates  a small  leakage,  is  close  enough 
to  constitute  a good  check  on  the  absence  of  serious 
leaks  at  the  time  of  this  test.  On  the  basis  of  this 
analysis  it  is  concluded  that  symmetrical  blowout  limits 
are  a sufficient  test  of  the  absence  of  air  leaks. 

L._i  . A 
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Departure  from  Ideal  reactor  conditions  A second 
source  of  error  In  these  experiments  is  the  departure 
of  the  chamber  from  Ideal  reactor  operation.  In  the 
theoretical  development  of  Chapter  III,  the  models  were 
based  on  the  assumption  of  an  Ideal  stirred  reactor  In 
which  operation  was  adiabatic  and  the  reaction  rate  was 
first  order  In  each  of  the  reactants.  In  principle, 
the  models  discussed  are  not  necessarily  thus  restricted. 

Any  relationship  of  combustion  efficiency  to  feed  rate 
or  blowout  loading  to  fuel/air  ratio  determined  experi- 
mentally on  a stirred  reactor  may  be  used  In  place  of 
the  theoretical  curves  used  In  the  previous  chapter. 

As  pointed  out,  however,  this  generalization  Is  restricted 
to  reaction  schemes  in  which  a single  reaction  controlls 
the  rate.  Where  successive  or  competitive  reactions  are 
possible,  one  cannot  use  a single  experimental  stirred 
reactor  performance  curve  to  compute  model  performance. 

In  comparing  combustor  performance  with  the  theo- 
retical model  prediction,  one  need  consider  only  the 
first  of  these  problems  (l.e.  heat  loss)  in  the  two- 
stabllity-limit  model,  since  the  reactors  of  the  model 
are  operating  in  parallel.  In  the  recirculation  model, 
both  heat  loss  and  multi-step  reaction  departures  are 
significant. 

In  the  present  work,  the  comparison  of  the  experl- 
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mental  results  with  the  theory  Is  based  on  experimental 
stlrred-reactor  performance  for  the  two-stability-llmlt 
model.  For  work  with  the  recirculation  model  the  Ideal- 
ized stirred  reactor  theory  Is  used.  Although  It  would 
be  clearly  preferable  to  use  an  experimental  stirred 
reactor  performance  curve,  such  data  have  not  yet  been 
reported  In  the  literature  In  sufficient  detail  to  be 
of  use  here.  To  the  author's  knowledge,  only  one  curve 
of  efficiency  vs  feed  rate  has  been  reported  (14),  and 
this  curve  was  at  a different  fuel/alr  ratio  than  the 
one  used  for  the  present  Investigation. 

There  are  ample  blowout  data  reported  in  the 
literature.  The  available  data  indicate  that  for 
hydrocarbons,  there  Is  no  substantial  difference  In 
the  performance  of  various  members  of  the  same  family 
(l.e.  saturates,  olefins,  etc.),  and  in  general,  most 
hydrocarbons  behave  approximately  the  same.  Acetylene 

and  hydrogen  are  different.  This  permits  analysis  of 
the  blowout  data  to  be  based  on  empirical  curves  of 
blowout  rates  vs  fuel/alr  ratio. 

Using  either  the  theoretical  performance  curves 
or  the  experimental  blowout  curves,  account  must  be 
taken  of  the  fact  that  the  heat  losses  from  this  reactor 
are  significantly  greater  than  those  for  the  stirred 
reactor  experimental  conditions.  In  the  present 
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experimental  work,  no  direct  measurement  of  the  heat 
losses  was  made.  An  estimate  of  their  magnitude  Involved 
accounting  for  losses  by  radiation  through  and  from  the 
Vycor  windows,  and  for  conduction-convection  losses 
directly  through  the  windows,  ho  correction  for  losses 
through  the  firebrick  was  made  since  it  was  felt  that 
it  was  sufficiently  thick  to  make  losses  In  this  direc- 
tion negligible  compared  with  the  other  heat  losses 
In  this  chamber. 

In  computing  the  heat  losses  from  this  chamber 
both  radiative  and  convective  heat  transfer  had  to  be 
considered.  Since  the  chamber  was  of  small  dimensions 
and  operated  at  low  pressures,  gas  radiation  could  be 
neglected.  Heat  transfer  from  the  gases  to  the  chamber 
walls  and  Vycor  windows  was  by  convection.  Heat  trans- 
ferred to  the  walls  and  roof  of  the  chamber  was  then 
radiated  to  the  Vycor  windows  and  the  outside  surround- 
ings seen  through  the  windows.  Heat  transferred  to 
the  windows  was  transferred  by  radiation  and  convection 
simultaneously  to  the  outside. 

In  computing  the  heat  transferred  by  radiation, 
the  chamber  was  divided  into  three  zones;  the  walls 
and  roof,  the  floor  and  the  windows.  The  walls  and 
roof  were  considered  to  be  a heat  source  at  a temperature 
T^ (different  from  the  gas  temperature) . The  windows 
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were  considered  to  be  a heat  sink.  The  floor  of  the 
chamber  represented  a refractory  surface,  returning  as 
much  heat  by  radiation  to  the  chamber  as  It  received 
from  the  walls  and  windows.  Vycor  windows,  made  of 
nearly  pure  quartz,  transmit  radiation  whose  wave  length 
Is  shorter  than  3.7^,  and  are  nearly  opaque  to  radiation 
of  longer  wavelengths.  Therefore,  the  absorptivitles 
and  emissivitles  of  Vycor  to  radiation  from  various 
temperature  sources  (the  walls  and  the  floor)  had  to 
be  accounted  for  in  the  calculations. 

Calculations  were  carried  out  for  the  three-jet 
chamber  using  a heat  transfer  coefficient  from  the  gases 
to  the  brick  of  10  Btu/hr  ft2  °F  at  the  blowout  velocity, 
and  it  was  assumed  that  the  heat  transfer  coefficient 
varied  linearly  with  velocity.  A convection  coefficient 
from  the  Vycor  windows  to  the  ambient  air  was  taken  to 
be  2 Btu/hr  ft2  °P.  The  brickwork  and  the  floor  of  the 

chamber  were  assumed  to  have  an  emlssivlty  of  1.0,  The 
Vycor  was  assumed  to  be  black  at  wave  lengths  above 
3.7m.,  and  transparent  at  wavelengths  below  this.  The 
reflectivity  of  the  Vycor  was  ignored  in  this  calculation. 
Figure  VI-l  shows  the  resultant  heat  loss  rate  from  the 
chamber  as  a function  of  temperature  and  relative  firing 
rate. 

Allowance  for  the  effect  of  heat  loss  in  the 
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calculation  of  chamber  performance  has  been  made  In  the 
past  (3,  14)  by  assuming  the  heat  losses  to  correspond 
to  an  effectively  lower  Inlet  temperature  to  the  reactor. 
This  technique  has  been  applied  principally  to  the  problem 
of  the  effect  of  temperature  on  blowout  limits.  Allowance 
for  the  effect  of  heat  losses  by  substitution  of  an 
artificial  inlet  temperature  requires  a previous  knowledge 
of,  or  an  Independent  calculation  of,  the  heat  loss. 

If  the  heat  loss  is  known  as  a function  of  temperature 
and/or  Beynolds  number,  a modified  approach  is  suggested. 
This  is  outlined  below. 

Consider  a reactor  for  which  the  heat  loss  rate 

may  be  expressed  by  the  relation: 

q = UA(t  - t ) VI-1 

o 

where  b is  the  heat  transfer  coefficient  based  on 
the  heat  transfer  area,  A,  the  temperature,  t,  of  the 
burning  gas  and  the  ambient  temperature,  tQ,  (also  the 
inlet  gas  temperature). 

The  simplified  heat  balance  for  the  adiabatic  reactor 
is  given  by: 

t - t + A . p VI-2 

a o af  r 

and  by  a process  analogous  to  that  used  to  obtain  Vl-2, 
the  heat  balance  for  a non-adiabatlc  reactor  may  be 
written: 

0 cp  (t  - tQ)  ♦ UA  (t  - tQ)  - Aaf  P 

k .. _ 
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where  la  the  adiabatic  flame  temperature  rise  for 

complete  combustion 

i? 

P is  the  fractional  bumedness  of  the  gases 

i 

G is  the  molal  feed  rate  to  the  chamber 

and  0^  is  the  molal  heat  capacity  of  the  gases  fed. 
fiewritlng  VI-3  leads  to  an  expression  parallel  to  VI-2; 

B 

A r 

* - *0  + » VI-4 
G°P 

; 

The  heat  balance  represented  by  Equation  VI-4  may  now 
be  combined  in  the  usual  manner  with  the  equation 
describing  the  stirred  reactor  behavior. 

Two  limiting  cases  for  VI-4  are  of  interest.  If 
the  heat  transfer  rate  is  controlled  by  forced  con- 
| vection  or  turbulent  diffusion,  then  to  a first  approxi- 

mation the  heat  transfer  coefficient,  U,  may  be  assumed 
to  be  directly  proportional  to  the  mass  feed  rate,  G. 

Then  the  denominator  of  the  fraction  in  IV-4  becomes  a 
constant,  and  the  heat  loss  rate  is  a constant  per- 
centage  of  the  total  heat  released,  l.e.  t “ V u*  p- 
This  is  nearly  equivalent  to  an  adiabatic  reactor  operating 
on  a fuel-air  mixture  diluted  with  additional  inert  gas 
such  as  nitrogen. 

If,  on  the  other  hand,  the  heat  transfer  coefficient 

L i 


is  independent  of  the  fleynolds  number,  such  as  in 
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radiation  or  conduction,  and  if  b may  be  treated  as  a 
constant  over  substantial  ranges  of  temperature,  then 
Equation  VI-4  assumes  the  form: 

t - t + p vx-s 

In  radiant  heat  transfer  the  assumption  of  linear  varia- 
tion of  heat  loss  with  temperature  over  substantial 
temperature  ranges  is  not  good,  but  it  is  Justified 
here  by  the  simplicity  thus  introduced.  Efficiency 
vs  burnedness  curves  for  a stirred  reactor  operating 
with  heat  losses  corresponding  to  the  two  limiting 
heat  loss  cases  above  are  shown  in  Figure  VI-2.  In 
this  figure  it  is  seen  that  the  shape  of  the  reactor 
performance  curves  near  blowout  is  not  greatly  affected. 

At  lower  feed  rates,  however,  the  heat  loss  in  the  con- 
stant U case  correspond  to  marked  cooling  of  the  gases 
and  results  in  quenching  the  reaction.  Thus  significant 
departure  from  the  adiabatic  stirred  reactor  curve  is 
predicted. 

A second  error  introduced  into  the  comparison  of 
the  experimental  reactor  performance  with  the  theoretical 
stirred  reactor  lies  in  the  departure  of  the  stirred 
reactor  from  the  simple  bimolecular  kinetic  mechanism. 

There  seems  to  be  general  agreement  in  the  literature 
that  the  apparent  overall  order  of  high-temperature 
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reaction  In  combustion  is  less  than  two.  Values  from 
about  1.3  to  2.C  have  been  reported,  with  no  general 
agreement  on  the  exact  value.  Values  between  1.6  and 
1.8  are  the  most  generally  accepted.  In  the  present  work 
the  pressure  exponent  of  the  reactor  at  blowout  was 
determined  for  four  different  fuel/air  ratios,  and  values 
of  "n"  ranging  from  1.3  to  2.0  were  found.  It  should  be 
pointed  out,  however,  that  this  variation  may  Involve 
the  simultaneous  effect  of  the  variation  of  reaction 
rate  with  composition  and  variation  of  mixing  rate  with 
Reynolds  number,  since  the  experiments  were  not  run  at 
a constant  Reynolds  number.  Even  when  experiments  are 
restricted  to  a single  feed  composition,  the  variation  in 
results  illustrated  due  to  pressure  may  De  to  an  unknown 
degree  due  as  well  to  Reynolds  number  variation  with 
its  associated  effect  on  mixing  and  to  variation  in  heat 
loss  with  velocity  change  in  feed  rate. 

A consequence  of  using  a second  order  equation 
when  the  actual  equation  is  of  lower  order  but  with  the 
activation  energy  unchanged  is  a change  in  the  shape  of 
the  curve  as  illustrated  in  Figure  VI-3.  In  this  figure 
the  theoretical  curve  for  a reaction  rate  first  order 
in  each  of  fuel  and  air  is  compared  with  the  curve  for 
a rate  which  is  first  order  in  air  and  0.8  order  in 
fuel.  These  latter  values  for  exponents  were  recommended 
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by  Longwell  and  Weiss  on  the  basis  of  their  studies. 

(The  lateral  displacement  of  the  curve  in  Figure  Vi-3 
is  irrelevant,  and  is  simply  a consequence  of  not 
changing  the  numerical  value  of  the  "collision  constant".) 

For  the  present  work  it  has  been  concluded  that 
although  accurate  experimental  data  on  the  loading  rate- 
efficiency  behavior  of  the  stirred  reactor  would  be 
desirable,  the  use  of  theoretical  curves  for  a compari- 
son of  the  model  behavior  with  the  experimental  behavior 
is  necessary.  Since  in  the  firing  rate  range  of  Interest 
the  difference  in  shape  between  idealized  stirred  reactor 
curves  of  p vs  (/  and  reactor  performance  corrected  for 
heat  loss  and  order  of  reaction  are  not  great,  it  does 
not  appear  Justifiable  to  Include  these  corrections. 

The  corrections  are  not  known  sufficiently  well  to  attach 
greater  certainty  to  a corrected  curve  than  to  an 
uncorrected  one. 

Errors  in  heasurement  A third  general  category 
of  error  to  which  this  experimental  program  was  suoject 
is  that  of  errors  in  measurement.  The  measurements  of 
significance  in  this  experimental  work  are: 

(1)  measurement  of  flow  rates 

(2)  measurement  of  chamber  dimension.  This  in- 
cluded the  height,  width  and  depth,  location 
and  size  of  the  nozzles  and  exit  flue 
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(3)  the  parallelism  of  the  Jets  when  the  three 

Jet  chamber  Is  being  used,  and  parallelism  of 
the  nozzle  block  assembly  and  the  sidewalls 
of  the  chamber. 

By  far  the  most  critical  measurements  were  those 
on  the  size,  spacing  and  alignment  of  the  feed- Jet 
assembly.  Even  slight  errors  here  were  sufficient  to 
cause  substantial  reductions  in  the  chamber  stability. 

Every  effort  was  made  to  obtain  uniformity  In  the  Jet 
assembly.  The  slot  openings  were  measured  with  a 
feeler  gauge,  and  the  Jet  alignment  was  checked  by  means 
of  a scale  with  1/100"  graduations  read  with  a magnifying 
glass.  In  the  final  analysis,  however  blowout  measure- 
ments provided  the  best  test  of  accurate  assembly. 

Although  it  is  believed  that  errors  of  this  type  do  not 
significantly  effect  the  blowout  and  efficiency  studies, 

[ I 

it  is  difficult  to  demonstrate  this  conclusively. 

Errors  in  the  measurement  of  flow  rate  and  of 
brickwork  dimensions  are  relatively  unimportant  in  this 
work.  Experience  with  calibrating  and  metering  with 
standard  square-edge  orifices  indicates  that  the  error 
in  metering  was  approximately  1#  in  mass  flow.  Errors 
of  the  dimension  were  kept  as  small  as  possible  with  the 
loose,  friable  type  of  refractory  used.  However,  brick- 
work dimension  errors  which  were  known  to  be  substantial 
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had  no  apparent  effect  on  the  chamber  stability;  so  It 
was  concluded  that  errors  resulting  from  this  source 
had  negligible  effect  on  the  overall  result. 

Analysis  of  the  Three-Jet  Chamber  Performance 

In  the  study  of  the  three- Jet  chamber,  two  stability 
limits  were  observed,  one  corresponding  to  a sudden  shift 
in  flow  pattern  with  blowout  of  the  flame  in  part  of  the 
chamber,  and  the  other  corresponding  to  total  blowout. 

A theory  explaining  this  behavior  qualitatively  has  been 
discussed  in  Chapter  III.  It  is  now  proposed  to  fit 
the  theory  quantitatively  to  the  data. 

Consider,  first,  approximations  which  can  be  made 
to  the  theory  of  Chapter  III.  As  the  blowout  condition 
is  approached,  (the  figure  representing  this  condition 
is  reproduced  for  conveniency  in  Figure  VI-4) , the 
stable  intersections,  corresponding  to  stable  combustion, 
approaches  the  noses  of  Branches  1 or  II,  depending  on 
the  mode  of  chamber  operation.  In  doing  this,  the 
corresponding  operating  point  on  the  well-stirred  reactor 
curve  for  the  zone  in  question  also  approaches  the  blow- 
out point  of  the  stirred  reactor  curve.  A reasonable 
approximation  which  can  be  made  to  this  situation  is  to 
say  that  the  blowout  of  a stable  zone  will  occur  when 
that  zone  has  reached  the  blowout  loading  rate  calculated 
on  the  basis  of  the  stirred  reactor.  Since  it  was 
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FIGURE  SI-4:  TWO- STABILITY-LIMIT  MODEL 

BLOWOUT  CONDITIONS 
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theoretically  predicted  that  blowout  occurred  In  this 
model  at  a somewhat  higher  burnedness  than  that 
corresponding  to  stirred  reactor  blowout,  this  assump- 
tion Is  not  rigorously  valid.  However,  In  this  range 
of  loading  rates,  burnedness  changes  very  rapidly  with 
loading,  so  that  only  a small  error  In  blowout  flow  rate 
Is  Introduced  by  this  assumption.  If  knowledge  of  the 
volumes  and  feed  rates  of  the  two  zones  In  this  chamber 
Is  postulated,  it  Is  possible  to  predict  whether  the 
chamber  will  burn  stably  in  the  two-zone  or  one-zone 
configuration,  or  be  completely  blowout  at  a given 
feed  rate. 

The  measurement  of  combustion  zone  volume  In  a 
turbulent  bunsen  cone  has  been  accomplished  by  Simon 
and  Wagner  (22)  using  flame  photographs  of  the  blue 
luminosity  of  a pre-mlxed  fuel/alr  mixture.  The  theo- 
retical validity  of  this  method  Is  uncertain,  since  some 
combustion  reactions,  such  as  carbon  monoxide  and 
hydrogen  oxidation  produce  nearly  non-lumlnous  flames, 
whereas  reactions  Involved  in  the  Initial  stages  of 
hydrocarbon  combustion,  generally  Involving  C-H  or  C-C 
bond  fracture,  give  rise  to  the  characteristic  color  of 
normal  flames.  If  the  hypothesis  is  that  the  initial 
stages  of  combustion  are  rate-controlling,  then  the 
measurement  of  the  luminous  zone  will  give  a reasonable 
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approximation  to  the  volume  of  combustion.  If  carbon 
monoxide  oxidation  is  postulated  to  be  the  rate- 
controlling  step,  then  this  measurement  may  be  grossly 
misleading. 

The  feed  rate  to  each  half  of  the  chamber  is 
postulated  to  be  a constant  fraction  of  the  total 
feed  and  in  addition,  a constant  fraction  of  the  fresh 
feed  is  assumed  to  grossly  bypass  the  active  volumes  in 
the  chamber.  In  the  case  of  a symmetrical  combustion 
pattern,  this  is  a reasonable  assumption.  For  assymetrical 
combustion  patterns  it  is  assumed  that  this  is  still 
valid  (i.e.  when  the  combustion  pattern  changes  due  to 
a blowout  of  one  side  of  the  chamber,  no  change  in  the 
fraction  of  the  fresh  feed  reaching  the  remaining  active 
zone  occurs) . 

Combining  these  assumptions  with  the  data  on  the 
photographically  measured  flame  volumes  leads  to  the 
conclusion  that  both  blowout  curves  for  the  three-jet 
chamber  should  be  displaced  from  the  well-stirred  reactor 
curve  by  a constant  multiplying  factor,  and  that  the  two 
blowout  curves  for  the  three-jet  chamber  should  be 
displaced  from  each  other  by  a constant  predicted  from 
the  difference  in  zone  volumes  only.  Since  Figure  V-6 
shows  the  zone  volumes  to  be  independent  of  fuel/alr 
ratio,  if  the  three-jet  chamber  were  composed  of  stirred 
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FIGURE  H-8:  COMPARISON  OF  STIRRED  REACTOR 

BLOWOUT  DATA  AND  THREE-JET  CHAMBER  DATA. 

ONE-SIDED  BLOWOUT  DATA  HAVE  BEEN  CORRECTED 
BY  THE  RATIO  OF  ZONE  VOLUME  AT  ONE-SIDED 
BLOWOUT  TO  ZONE  VOLUME  AT  TOTAL  BLOWOUT. 


I 


- 146 


reactor  elements,  one  would  expect  that  the  performance 
curves  of  the  well-stirred  and  the  three- Jet  combustors 
should  be  displaced  from  each  other  by  a constant  which 
is  independent  of  fuel/air  ratio.  This  is  not  observed, 
as  may  be  seen  in  Figure  VI-5.  This  deviation  from  the 
expected  behavior  may  result  from  the  inadequacy  of  the 
assumptions  above  or  it  may  reflect  the  inability  of 
the  stirred  reactor  to  describe  operation  in  the  active 
zones. 

The  model  also  predicts  that  the  two  curve  for  the 
three- Jet  chamber  should  be  separated  from  each  other 
by  a constant  which  may  be  derived  from  the  zone-volume 
data.  This  may  be  argued  formally: 

(1)  If  Vx  represents  the  volume  of  the  smaller 
zone  (of  the  two  active  zones), 

(2)  V8  represents  the  volume  of  the  large  zone 
when  there  is  only  one  active  zone,  and 

(3)  If  Y represents  the  fraction  of  the  feed 
reaching  each  zone  which  is  the  same  for 
each  active  zone  and  independent  of  flow 
rate  and  zone  volume, 

then,  Vx/GY  = Lx,  the  loading  factor  for  active  zone  1 
and  Va/GY  » La,  the  loading  factor  for  active  zone  2. 

When  an  active  zone  blows  out,  it  is  postulated 
that  it  is  at  the  critical  loading  factor,  of  the 
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stirred  reactor  at  blowout.  Hence  • l^^  ’ or 
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BOj  = Vj 
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In  Figure  VI-5,  the  stirred  reactor  curve  has  been 
compared  to  the  three-jet  chamber  stability  curves  by 
plotting  an  ordinate  of  (N^/VP1,^J3j0/  ” 

vs  generalized  fuel  fraction.  It  is  immediately  apparent 
that  the  three-jet  chamber  does  not  follow  the  expected 
performance  based  on  that  of  the  well  stirred  reactor. 
However,  the  two  curves  representing  the  blowout  limits 
of  the  three-jet  chamber  are  approximately  displaced  from 
each  other  by  the  expected  amount.  The  significance  of 
the  unexpected  relative  performance  is  not  clear.  It 
might  be  argued  that  this  resulted  from  air  leakage  into 
the  chamber.  This  argument  is  countered,  however,  with 
the  evidence  discussed  at  the  outset  of  this  chapter 
that  the  data  are  essentially  free  from  air  leakage 
errors.  A second  possible  reason  for  thiB  behavior  1b 
that  the  difference  in  the  combustion  characteristics 
of  propane  (used  here)  and  octane  (used  by  Longwell) 
is  sufficient  to  cause  this  behavior.  This  argument  is 
not  supported  by  any  evidence.  Hecent  work  done  at 
M.  I.  T.  (5)  as  yet  unpublished)  using  pure  propane  as 
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the  fuel  In  a well-stirred  reactor  gave  the  results 
Indicated  In  Figure  II-8,  showing  no  difference  In 
stirred  reactor  performance  which  could  account  for 
the  observed  difference  In  the  three-jet  chamber  and 
the  stirred  reactor  performance  characteristics. 
Therefore,  one  Is  led  to  the  conclusion  that  the  ob- 
served difference  represents  a characteristic  of  the 
operation  of  the  active  zones  not  accounted  for  by 
the  stirred  reactor  performance. 

The  efficiency  measurements  made  on  this  chamber, 
though  not  extensive  enough  to  fit  to  any  model,  are 
useful  for  showing  that  the  chamber  operates  at  a 
reasonable  efficiency  level.  From  these  measurements 
It  Is  also  possible  to  estimate  the  quantity  of  gross 
bypassing  of  the  combustion  zones.  If  one  takes  a 
value  of  0 of  90#  for  the  burnedness  of  each  of  the 
zones  (an  average  figure  representative  of  the 
efficiency  data  reported  by  Longwell  (14)),  then  for 
the  gross  efficiency  to  be  about  70%  leaving  the  three- 
Jet  chamber  requires  approximately  22#  of  the  total 

feed  to  the  chamber  to  bypass  completely  the  combustion 
zones.  The  balance  (78#  of  the  total  feed)  enters  the 
combustion  zones  and  burns  to  a 90#  efficiency.  Were 
the  ratio  of  gross  loading  factor  for  the  three-jet 
chamber  at  blowout  to  the  well-stirred  reactor  at 
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blowout  known  with  certainty,  It  would  be  possible  to 
calculate  the  fraction  of  volume  In  each  recirculation 
zone  which  was  active  and  that  fraction  which  was 
Inactive.  However,  because  of  the  oehavlor  of  the 
ratio  (Figure  VI-5)  no  constant  value  for  volume  ratio 
1 s consistent  with  the  rest  of  the  theory. 

In  addition  to  the  results  reported  here,  an 
attempt  was  made  to  use  Schlieren  photography  to  gain 
an  insight  into  the  mechanism  of  chamber  behavior. 

These  efforts,  unfortunately,  were  unsuccessful,  due 
particularly  to  extraneous  patterns  caused  by  the  hot 
glass  and  by  the  warm  turbulent  air  ambient  to  the 

chamber  in  the  vacuum  vessel  which  were  In  the  optical 
path. 

Analysis  of  the  One-Jet  Chamber  Performance 

The  theory  developed  to  explain  the  two-stability 
limits  observed  in  the  three-jet  chamber  implied  that  a 
similar  behavior  should  be  observed  In  the  one-jet 
chamber.  To  test  this  hypothesis  a one-jet  chamber  was 
built  and  tested.  The  first  design  tried  was  a chamber 
of  width  l/3rd  that  of  the  original  three-jet  chamber. 
As  mentioned  previously,  this  proved  to  be  too  narrow 
since  its  stability  llmt  was  unreasonably  low.  Stable 
operation  was  achieved  by  making  the  chamber  2W  wide 
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rather  than  the  original  1*  width,  still  narrower, 
however,  than  the  original  3"  wide  three- Jet  chamber. 
Stable  combustion  was  easily  established  in  the  one- 
jet  chamber,  though  the  combustion  pattern  was  not 
symmetrical  (see  Figure  V-7  and  8).  It  is  postulated 
that  this  asymmetrical  combustion  pattern  is  a con- 
sequence of  a basic  instability  of  the  cold  flow  pattern 
in  which  there  is  no  stable  symmetrical  flow  pattern 
regardless  of  how  slowly  the  chamber  is  fed.  Because 
no  symmetrical  combustion  pattern  existed,  there  was  no 

i 

possibility  of  the  dual  stability  limit  found  in  the 
three-jet  chamber.  From  the  behavior  of  the  combustion 
pattern  as  blowout  was  approached,  a model  analogous 
to  the  recirculation  model  discussed  in  Chapter  111 
was  suggested. 

To  compare  the  narrow  one- Jet  chamber  with  the 
theoretical  recirculation  model  of  Chapter  III,  efficiency 
vs  loading  rate  data  were  required.  Blowout  data  were 
also  taken  to  compare  the  performance  of  this  chamber 
with  that  of  the  stirred  reactor. 

The  comparison  of  the  stirred  reactor  data  and  the 
one- Jet  chamber  blowout  data  is  shown  in  Figure  VI-6. 

These  data  show  a deviation  from  the  stirred  reactor 
which  is  not  accounted  for  by  the  model.  As  in  the 
case  of  the  two-stability-llmlt  model,  the  recirculation 
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model  also  predicts  that  there  should  be  a constant 
difference  between  the  stirred  reactor  blowout  and 
the  experimental  chamber  blowout.  The  trend  observed 
here  may  reflect  the  departure  of  the  stirred  reactor 
from  the  assumed  simple  kinetic  picture. 

The  recirculation  models  and  the  actual  performance 
of  a combustor  are  compared  by  the  efficiency  vs  loading 
rate  curve.  Given  this  information  about  a combustor's 
performance,  a value  of  R,  the  recirculation  ratio,  and 
Va/V j may  be  chosen  so  that  the  theoretical  behavior 
will  match  the  model  behavior.  Roth  the  position  of  the 
curve  (blowout  loading  rate)  and  the  chamber  outlet 
efficiency  are  adjusted  in  these  models  by  the  use  of 
gross  bypassing  and  of  Inactive  volumes.  These  two 
variables  do  not  affect  the  basic  shape  of  the  efficiency 
loading  rate  curve.  In  addition,  data  on  the  size  and 
composition  of  the  individual  zones  will  make  correlation 
between  theory  and  experiment  more  meaningful  since 
fewer  constants  are  then  to  be  determined  from  the  data. 
Such  data  were  obtained  by  the  use  of  sampling  traversed 
across  the  chamber  and  by  photographic  measurements.  The 
behavior  of  the  combustor  pattern  as  the  flow  approached 
the  blowout  limit  suggests  that  the  bright  zone  of 
Figure  V-7  be  taken  as  the  central  core  (of  volume  Vi), 
and  the  dark  eye  of  the  pattern  be  taken  as  the 


model  Incorporating  this  assumption  might  fall  to 
provide  a basis  for  correlating  the  data.  As  may  be 
seen  from  the  traverses  of  composition  within  the  chamber, 
the  change  in  burnedness  in  moving  between  the  various 
zones  in  the  chamber  is  gradual,  and  a model  allowing 
for  a more  nearly  continuous  change  would  probably  be 
in  better  agreement  with  the  observed  data. 

If  one  considers  a circulatory  flow  pattern  in 
which  gases  may  enter  and  leave  at  the  circumference, 
the  inner  zone  (or  zones),  further  removed  from  the 
circumference  of  the  combustion  pattern,  receive  a 
smaller  fraction  of  the  fresh,  combustible  gases,  and 
consequently,  are  less  "heavily  loaded"  than  those 
outer  zones  nearer  to  the  circumference.  In  such  a 
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combustor  as  the  feed  rate  Is  Increased,  the  outer  zones 
will  blow  out  first  because  of  their  heavier  loading 
rate,  and  the  inner  zones  will  consequently  receive  a 
less  well-burned  feed  since  the  outer,  shielding  zones 
now  contain  more  unburned  gases.  In  a conceptual  model 
such  as  this,  as  the  feed  rate  is  increased  the  region 
of  most  intense  burning  rate  will  gradually  recede 
from  the  circumf erence  of  the  combustion  pattern  to 
the  center  until  the  inner-most  zone  of  the  model  blows 
out,  causing  total  extinction  of  the  flame. 

Although  a two-zone  approximation  to  this  behavior 
has  already  been  seen  to  be  unsatisfactory,  treatment 
of  this  model  by  integration  of  the  performance  of 
differential  elements  would  be  a formidable  task.  A 
satisfactory  approximation  considers  recirculation  models 
involving  a small  number  of  identical  elements,  each 
successive  recirculation  zone  being  fed  by  its 
immediate  predecessor  , and  receiving  recirculated 
gases  from  its  immediate  successor,  burned  gases  from 
a typical  zone  are  split  in  two,  one  portion  feeding 
the  recirculation  zone  which  is  its  Immediate  predecessor 
and  the  other  feeding  the  zone  which  is  Ita  Immediate 
successor.  A diagram  of  such  a model  is  shown  in 
Figure  VI-7.  Since  this  extended  recirculation  model 
is  an  approximation  to  the  circumferentially-fed 
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recirculation  model  described  above,  the  recirculation 
ratios  and  relative  zone  volumes  at  each  stage  of  the 
model  are  independent  of  previous  model  stations. 

The  volumes  (comparable  to  finite-difference  intervals 
in  a conventional  finite-difference  solution)  are 
taken  equal  to  each  other  (or  in  the  case  of  cylindrical 
coordinates,  for  example,  may  be  proportional  to  the 
distance  from  the  origin)  and  the  recirculation  ratios 
representing  the  diffuslonal  transport  of  matter  in 
the  chamber,  are  equal  (or  again  may  be  proportional  to 
the  distance  from  the  inner  zone) . 

For  prediction  of  combustor  performance  using 
this  model  only  variable,  the  recirculation,  is  in 
principle  required  to  find  the  behavior  of  the 
circumferentially-fed  recirculation  model.  A character- 
istic recirculation  ratio  having  been  chosen,  a model 
with  a sufficient  number  of  volumes  is  calculated  to 
make  its  behavior  practically  the  same  as  the  infinite- 
number-of -zones  case.  In  practice,  the  calculation  of 
a model  containing  more  than  four  or  five  zones  is  too 
laborious  to  consider  by  hand  calculation  techniques, 
and  for  models  containing  this  small  number  of  zones, 
the  behavior  of  the  model  is  strongly  dependent  on  the 
number  of  zones  in  the  model.  This  is  shown  in  Figure 
VI-8,  where  three  cases  are  calculated  for  three 
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different  models,  {for  this  calculation,  h In  each 
case  was  so  chosen  to  bring  the  blowouts  to  about  the 
same  burnedness.)  In  these  calculations  both  the  volume 
sizes  and  the  recirculation  ratios  have  been  taken  to 
be  independent  of  position  in  the  model. 

The  effect  of  the  number  of  volume  units,  p,  in 
a model  is  better  shown  by  comparison  of  the  models  on 
a basis  of  (total  feed)/ (first  volume)  ratio.  Since 
at  high  burnedness  the  recirculation  zones  contribute 
little  to  the  performance  of  a chamber,  performance 
curves  for  models  of  different  p's  will  coincide  when 

plotted  on  this  basis.  Using  this  basis  the  behavior 
of  a model  of  constant  recirculation  ratio  and  initial 
volume  but  of  varying  number  of  zones  is  shown  in 
Figure  VI-9. 

When  the  calculated  performance  curves  are  plotted 
in  this  manner,  two  Interesting  observations  emerge; 

(l)  additional  zones  do  not  cause  corresponding  addi- 
tional humps  in  the  performance  curve  (as  did  tne  f irst 
recirculation  zone) , and  (2)  increasing  the  numDer  of 
zones  of  the  model  increases  the  stability  of  the 
residual  flame,  but  does  not  affect  the  stability  of 
the  first  zone  significantly.  This  suggests  that  if  a 
residual  flame  zone  exists  in  a combustor,  there  will 
be  a sudden  change  In  combustion  efficiency  as  the 
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FIGURE  SI-9:  EXTENDED  RECIRCULATION  MODEL,  FIRST  ZONE  VOLUME  AS  BASIS,  R'l.O 
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main  combustion  zone  blows  out,  leaving  the  residual 
flame,  but  further  Increases  in  flow  rate  will  only 
cause  the  residual  flame  to  recede  gradually  from  the 
main  combustion  zone.  No  additional  partial  blowout 
of  zones  within  the  residual  flame  are  to  be  expected. 

The  possibility  of  designing  combustors  of  much 
Improved  stability  is  also  suggested,  by  deliberately 
allowing  for  a zone  which  is  sufficiently  sheltered 
and  which  1b  deep  enough  to  allow  the  residual  flame 
to  recede  into  it,  a substantial  distance,  a chamber 
equivalent  to  a high-p  model  is  obtained,  having  a 
high  ratio  of  ultimate  blowout  loading  to  initial 
blowout  loading. 

The  use  of  this  extended  recirculation  model  will 
now  permit  a fitting  of  the  experimental  data  obtained 
in  the  present  experiments.  There  are  two  variables 
in  the  model,  namely,  the  recirculation  ratio  and  the 
number  of  recirculation  zones.  In  addition,  allowance 
must  be  made  for  "dead”  or  inactive  zones  in  the  chamber 
and  for  gross  bypassing  of  the  combustion  zones. 

Only  the  variables  of  recirculation  ratio  and 

number  of  zones  will  affect  the  shape  of  the  performance 
curve.  Increasing  recirculation  ratio  will  cause  a 
change  in  the  sharpness  of  the  first  blowout  limit 
without  causing  much  change  in  the  ratio  of  the  two 


blowout  limits  (Figure  Vl-10) . An  Increase  In  the  number 
of  zones  will  increase  the  relative  stability  of  the 
recirculation  zone  (Figure  VI-9). 

A change  in  gross  bypassing  of  gases  around  the 
combustion  zones  affects  both  the  outlet  burnedness 
(by  dilution)  and  the  relative  stability  with  respect 
to  the  stirred  reactor.  It  does  not,  however,  affect 
the  shape  of  the  curve  at  all.  The  addition  of  "aead" 
space  to  the  combustion  chamber  model  permits  the 
apparent  volumetric  heat  liberation  rate  to  be  adjusted 
so  that  the  model  prediction  is  in  conformity  with  the 
stirred  reactor  data.  The  effectiveness  of  this  type 
of  model  in  correlating  the  data  on  chamber  efficiency 
is  demonstrated  in  Figure  VI-ll,  where  the  experimental 
data  are  compared  with  a model  having  a recirculation 
ratio  of  one,  two  or  three  recirculation  zones  plus 
bypassing  and  dead  space.  By  plotting  these  data  on 
log-log  coordinates,  movements  corresponding  to 
addition  of  dead  space  or  bypassing  correspond  to 
either  a simple  lateral  and/or  vertical  displacement 
of  the  two  curves.  (The  raw  data  have  been  shifted 
by  a constant  factor  of  3 x 10  3 to  facilitate  plotting 
on  this  graph.) 


LOADING  RATE,  NA/VP'® 

FIGURE  SI- II:  COMPARISON  OF  MODEL  THEORY  WITH  EXPERIMENTAL  COMBUSTOR 
PERFORMANCE. 
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Summary 

The  demonstration  of  the  ability  of  simple  theo- 
retical combustion  chamber  models  to  reproduce  the 
behavior  of  an  experimental  chamber  in  a semi- 
quantitative  manner  is  of  great  significance.  Previous 
models  of  combustion  chamber  performance  have  not,  as 
yet,  been  used  in  conjunction  with  an  experimental 
program.  In  this  investigation  several  theoretical 
models  were  considered,  and  the  calculated  performance 
of  the  models  was  compared  with  the  experimentally 
measured  performance  of  a simple  chamber  which  demon- 
strated characteristics  common  in  many  practical  chambers. 

Ultimately  one  would  hope  to  use  models  of  this 
nature  to  predict  the  performance  of  any  combustion 
chamber.  Using  a relatively  small  number  of  variables 
it  has  been  possible  to  generate  a wide  variety  of 
performance  curves  for  various  combustor  models.  Thus, 
for  instance,  in  moderately  well-stirred  systems,  typical 
of  can-type  combustors,  the  extended  recirculation  model 
suggested  seems  to  be  applicable  to  the  prediction  of 
residual  flame  phenomena,  and  the  use  of  the  two- 
stablllty-limlt  model  is  suggested  for  the  analysis  of 
pressure  Interaction  phenomena.  Hopefully  extension 
of  th*»se  modeling  techniques  will  lead  to  a general 
xer.hoa  of  combustion  chamber  analysis,  and  possibly 
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CHAPTEH  VII 

CONCLUSIONS  AND  BECOMMENLATIONS 


Conclusions 

1.  The  theoretical  models  of  combustion  chambers  have  been 
Investigated.  These  models  have  predicted  qualitatively 
phenomena  observed  in  actual  combustors. 

a)  A two-zone  model  of  stirred  reactors  in  parallel 
was  used  to  predict,  the  effects  of  pressure  interaction 
between  two  combustion  zones. 

b)  A two-zone  recirculation  model  in  which  one  zone 
was  fed  only  partially  burned  products  produced  in  the 
first  and  the  first  zone  received  a mixture  of  fresh  gases 
and  burned  products  from  the  second  zone  was  used  to  in- 
vestigate recirculation  phenomena.  This  model  showed  the 
conditions  for  residual  flame,  hysteresis,  and  flickering 
flame  phenomena. 

c)  An  extended  recirculation  model  was  investigated 
which  allowed  for  a more  continuous  variation  in  composi- 
tion in  the  combustor  than  could  be  obtained  in  a two-zone 
model.  This  model  showed  phenomena  similar  to  those  observed 
in  the  two-zone  model. 

2.  An  experimental  combustion  chamber,  fed  by  a two- 
dimensional  Jet  of  a pre-mlxed  fuel/air  mixture  was  built 
and  operated  to  test  the  models  of  combustion  chambers 
quantitatively . 
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a)  All  the  models  predicted  that  stirred  reactor  blow- 


out and  combustion  chamber  loading  rates  should  be  displaced 
from  each  other  by  a constant  factor.  This  was  not  observed. 
This  deviation  was  postulated  to  result  from  the  breakdown 
of  Certain  assumptions  and  from  departure  of  the  reaction 
kinetics  from  a mechanism  controlled  by  a single  reaction. 

b)  The  two-stability-llmit  model  predicted  the  dis- 
placement of  the  two  stability  limits  of  the  experimental 
reactor  to  a reasonable  approximation. 

c)  The  overall  combustion  efficiency  of  the  chamber 
was  measured  by  gas  sampling.  These  data  were  correlatable 
by  the  extended  recirculation  model. 

d)  The  combustion  efficiency  of  the  chamber  contents 
was  determined  by  gas  sampling.  These  data  showed  a 
qualitative  behavior  in  agreement  with  the  models.  No 
attempt  was  made  to  fit  the  data  quantitatively. 

e)  The  effect  of  pressure  on  the  blowout  loading, 
G/VP1*®  was  determined  experimentally.  An  average  pressure 
exponent,  "n"  of  1.6  was  found  which  compares  well  with 
reported  literature  values  of  1.6  to  1.8.  A trend  of  "n" 
was  observed  with  variation  in  fuel/air  ratio.  However, 

no  significance  was  attached  to  this  trend  because  of 
unequal  Reynolds  numbers  between  experiments. 


Recommendations 

(l)  The  prediction  of  combustor  performance  by  the  use 


of  these  models  requires  experimental  data  on  combustion 
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efficiency  vs  loading  rate  for  the  well-stirred  reactor. 
There  is  not  yet  available  an  adequate  amount  of  these 
performance  data  for  the  stirred  reactor.  It  is  recommended 
that  an  experimental  program  be  undertaken  to  determine 
these  data  and  to  determine  whether  they  are  in  agreement 
with  the  performance  predicted  from  the  stirred  reactor 
theory . 

(2)  Although  the  experimental  program  undertaken  gives 
some  support  to  the  theories,  unexplained  phenomena  still 
exist.  It  is  recommended  that  further  experimental  work  be 
undertaken  to  establish  methods  for  fitting  the  combustion 
chamber  phenomena  to  the  theory.  Investigation  of  certain 
special  phenomena,  is  also  recommended: 

a)  To  fit  the  present  data  to  existing  stirred  reactor 
data  it  is  necessary  to  postulate  substantial  bypassing  of 
the  combustion  zones.  An  experimental  verification  of  the 
extent  of  this  bypassing  is  recommended. 

b)  The  theory  has  only  been  fitted  at  one  fuel/air 
ratio  to  experimental  performance  data.  It  is  recommended 
that  other  fuel/alr  ratios  be  Investigated  and  that  par- 
ticularly, the  prediction  of  the  effects  of  fuel/air  ratio 
changes  by  the  theoretical  models  be  studied. 

3.  hata  on  the  stirred  reactor  nas  indicated  some  lack  of 
complete  mixing.  An  investigation  of  the  stirred  reactor  is 
recommended  to  determine  what  the  consequences  of  this  lack 
of  mixing  are,  and  to  determine  which  zones  in  a combustion 
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APPENDIX  A 

DERIVATION  OF  THE  WELL-STIRRED  REACTOR  EQUATION 

The  derivation  of  a generalized  well-stirred  reactor 
equation  has  already  been  discussed.  In  the  literature 
several  more  detailed  discussions  of  the  exact  algebraic 
form  of  the  equations  have  been  given  for  specific  postu- 
lated stoichiometric  schemes  (3,  14) . Since  it  was  felt 
that  the  use  of  complex  stoichiometric  schemes  did  not 
yield  a justifiable  increase  in  accuracy  of  the  model 
predictions,  the  derivation  of  the  algebraic  relation 
between  burnedness  and  loading  rate  used  in  this  thesis 
was  restricted  to  the  simplest  possible  case  in  which  no 
intermediate  products  are  postulated. 

The  stoichiometric  equation  of  the  reaction  represent- 
ing the  overall  combustion  reaction  was  written  for 
propane , 

C3Hq  + 5 0a  — ► 3 C0a  +4  Ha0  A-l 

The  heat  of  reaction  was  assumed  to  be  sufficient  to  give  a 

2270°K  temperature  rise  for  a stoichiometric  mixture,  and 

the  heat  capacities  of  the  reactants  and  products  were 

assumed  to  be  equal  and  independent  of  temperature  so  that 

the  fraction  of  fuel  consumed,  0,  was  related  to  temperature 

for  a stoichiometric  mixture  by: 

T - T 

0 = 8-  A-2 

2270 

^ -iJ 
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In  the  simplest  description  of  the  stirred  reactor  the 
identifiable  components  of  the  mixture  at  any  stage  of 
completion  of  combustion  are  burned  products  plus  fuel  and 
oxygen;  therefore  only  a single  variable,  p,  is  necessary 
to  specify  the  composition  of  the  mixture.  If  the  fraction  (p 
of  the  stolchiometrlcally  required  fuel  is  used,  then  for 
each  5 moles  of  oxygen  the  moles  of  each  component  are  given 
as  follows: 

c3h8  <p  (1  - p) 

0a  5(1  -<pp) 

Na  18.8 

Ha0  and  C0a  7 cp  P 

Total  mols  23.8  + (ft  (l  + p) 

If  a rich  fuel-air  mixture  is  used,  there  should  be  no 
unburned  fuel  appearing  at  completion  of  combustion  (as 
would  be  called  for  by  the  above  stoichiometry  for  <p>l). 
Thus  for  rich  mixtures  the  stoichiometric  equation  A-l  is 
altered  to  give: 

<P  C 3Hg  + 5 Qj  ->  <p[3(GO  + C0a)  + 4(Ha  + Ha0)] 

A-3 

In  this  equation  the  summation  of  oxygen  atoms  on  the 
right-hand  side  of  the  equation  must  be  equal  to  10. 

Longwell  and  Weiss  (24)  and  Baker  (3)  have  both  used  this 
assumption  to  fit  data  on  the  rich  side,  and  in  both  cases 
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the  authors  assumed  that  the  water-gas  equilibrium  could  be 
used  to  calculate  the  composition  of  the  products.  The 
present  calculations  were  simplified  by  assuming  that  the 
heat  of  reaction  contributed  by  the  water-gas  reaction  was 
negligible  compared  to  the  heat  released  by  ti.e  initial 
combustion  step.  With  this  assumption  it  is  not  necessary 
to  specify  the  composition  of  the  mixture  to  calculate  the 
temperature.  The  stoichiometry  of  reaction  A-3  yields  the 
following  species  (basis  - 5 moles  oxygen) : 

c3h8  <£  (1  - 0) 

08  5 (l  <=  p) 


products 

total 


18.8 
7<P  p 

23.8  + 4>  + (6ep-5)  p 


For  fuel-air  mixtures  other  than  stoichiometric  the 
heat  balance  equation,  A-2,  relating  temperature  and 
burnedness,  must  be  modified.  In  principle  an  enthalpy 
balance  should  be  made  to  account  for  dissociation  products 
and  water-gas  equilibrium,  especially  for  rich  mixtures. 
However,  calculations  may  be  greatly  speeded  by  the 
simplifying  assumptions  of  (l)  a constant  heat  capacity  of 
the  gases  and  (2)  ignoring  the  energy  of  the  water  gas 
reaction  and  dissociation  products.  For  the  present 
calculations  the  following  temperature-burnedness  relations 
were  used: 
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500. 000  <P 
ZOcp  + 190 


p * T - TQ  (lean) 


A-4.1 


500.000 
30  <p  + 190 


P * T - Tq  (rich) 


A-4.2 


It  should  be  noted  that  for  a stoichiometric  fuel-air 
mixture  both  rich  and  lean  schemes  are  identical. 

The  composition  of  the  gases  in  a stirred  reactor  in 
which  this  combustion  occurs  is  a function  of  the  firing 
rate  to  the  reactor.  This  function  is  derived  from  the 
reaction  rate  equation  applicable  to  this  reaction.  For 
this  case  here,  only  a single  equation  is  required,  and 
this  is  assumed  to  be  an  equation  of  first  order  in  each 
of  the  reactants.  The  temperature  dependence  of  the  rate 
is  indicated  by  the  conventional  exponential  term: 


— = /T  k ' e"E/HT . c , Cp  A-5 

Vdt  °8  ¥ 

From  the  stoichiometry  discussed  the  concentrations  of 
the  fuel  and  oxygen,  and  Cp,  may  be  found  as  the 
respective  mole  fraction  terms  times  the  molar  density  of 
the  gas  determined  from  the  perfect  gas  law; 


'03 


. 5 (1  -<f)B)  p_ 

23.8  +(f(p  + 1)  HT 

= Cf>(l  - 6)  P_ 

23.8  + cp  (p  + 1)  HT 


lean 


A-6.1 


23.8  + <P  + (6<P-5)P 


A U ---£l 


;j 


) rich 


A-6.2 


23.8  + cp  +(6^P  -5)p  ^ 


Since  the  rate  of  fuel  consumption  is  given  by  G«p  ; 

0_  . lC|~yaT  SfU-WU-P)  lean  A-7.1 
VP8  t3/2  p [23.8'  + cP  (l+p)J  2 


G 


- *>-E/gT  5 


(i  - e): 


rich  A-7.2 


VP2  T 3^2  23.8  + <p  + (6<p-5)p 

For  calculatlonal  purposes  it  is  desirable  to  normalize 

equations  A-7  so  that  those  changes  in  operating  variables 

which  counteract  each  other  will  be  evident.  At  the  same 

time  it  is  desirable  that  each  of  the  groups  chosen  be 

independent  l.e.  so  that  a change  in  a factor  in  one  group 

will  not  affect  the  values  of  the  other  groups. 

In  an  experiment  to  determine  the  reaction  rate  constants 

-E/HT 

k1  and  E,  the  rate  constant  of  the  reaction,  k = k 'e  , 
will  be  found  at  some  high  temperature,  and  by  varying  the 
temperature  k1  and  E may  be  found  separately.  Since  this 
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experiment  must  be  carried  out  In  a temperature  region  where 
the  reaction  proceeds  with  a reasonable  velocity,  knowledge 
of  k*  at  room  temperature  will  be  very  uncertain.  This 
suggests  that  an  appropriate  normalization  factor  would  be 

the  rate  constant,  k,  determined  at  some  high  temperature 
designated  as  the  key  temperature,  TK.  Temperatures  would 
also  be  normalized  by  TR.  This  normalization  procedure 
yields  the  reduced  equation; 

G_  B»TKVg  _ ,-B/Mk  [(l/Tjt-0  sdd.ecPHl-B) 

vpa  k «e-E/HTK  tJ5/2  [23.8  + (l-0)Ja  p 

A-8 

Por  this  Investigation  a basic  stirred  reactor  curve 
calculated  from  A-8  was  used  to  compare  the  performance 
curves  of  the  various  models.  This  curve  is  based  on  the 
following  constants: 

<P  * 1.0 

Tk  ■ 2000  °K 

T = 500 °K 

o 

Taf  = 2772°K 

e/htk  = 10 

A plot  of  (G/VPa) (B*TK/k,e"E/aTK)  vs  ^ is  shown  in 


Figure  A-l  for  these  constants.  This  curve  (plotted 
elsewhere  in  this  thesis  as  G/VPa  or  G/VP11  vs  (3  is  used  to 
generate  all  of  the  other  model  performance  curves  calculated. 
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In  the  model  Investigating  the  effect  of  non-uniform  fuel- 
air  ratio  on  combustor  performance  curves  for  stirred 
reactors  being  fed  a rich  fuel-air  mixture  were  required. 
For  this  purpose  Equation  A-7.2  was  normalized  and  curves 
calculated  for  a 4?  of  1.2  and  1.5  using  all  other  constants 
shown  above  the  same. 
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MULTI-STEP  KINETIC 


Calculatlonal  techniques  discussed  In  the  main  body 
of  this  thesis  have  assumed  that  the  reaction  rate  of 
burning  gases  Is  limited  by  a single  chemical  reaction 
throughout  the  composition  range  of  Interest.  It  must  be 


recognized  that  this  Is  an  assumption  Justified  only  by 
the  resulting  simplifications  In  the  theory  and  the  hope 
that  It  will  Introduce  only  a small  error  In  the  predicted 
model  performance.  In  fact  It  Is  possible  to  make  rigorous 
allowance  for  multi-step  reaction  although  for  more  than 
two  or  three  steps  the  calculatlonal  work  becomes  prohibitive. 

Calculations  of  multi-step  cases  are  greatly  simplified 
by  the  use  of  graphical  techniques  to  solve  the  stirred 
reactor  equations.  To  Illustrate  the  use  of  such  techniques 
a typical  two-step  mechanism  will  be  considered.  The  two 
steps  postulated  are: 

C3Hq  + 7/2  0a  ► 3 CO  + 4 HaO  B-l 

3 CO  + 3/2  0*  — 3 C08  B-2 

For  this  postulated  mechanism  two  variables  are  required  to 
specify  the  composition,  representing  respectively  the 
progress  of  reactions  B-l  and  B-2.  It  Is  assumed  that 
operation  Is  adiabatic  so  that  the  temperature  Is  fixed 
by  the  composition  and  Initial  temperature  of  the  unreacted 
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expression. 

dN, 


C3H, 


Vdt 


a . /Tk,'  e-*i/BT  c 


CaH, 


8 


B-3 


dN, 


CO 


Vdt 


3 /Tki'  e~El/RT  CCa 


H8 


- /T  ka'  e“Ea/RT  CCQ-  C, 


B-4 


CO  Oa 

Substitution  for  appropriate  concentration  terms  yields: 

(fc/^  - /Tki  e"El/RT  (1  - 6)  P_ 

£ BT  B-5 


^(p-7)  » 3 Ik1  e"El/RT  ^ 

- T k*  e"E*/RT  3(5/»  -7/26-3/27) (6-7)  £_ 


B-6 


Belative  reaction  rate  constants  of  reactions  B-l  and  B-2 
as  well  as  a reduced  feed  rate  are  now  defined: 

BT*  ki  e_El/RTK 


Belative  Beactlon  Bate,j** 


pk,  .-k»/htk 


and 


0 r B*  T, 


« — [ 
ur>8  L 


3/2 


r VPa  L ka  e"Ea'RTK 
Equation  B-5  thus  becomes: 


] 


orap 


e.(E,/BV  (i/ir  -i) 

r l 87 
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fuel-air  mixture.  It  will  be  assumed  that  the  oxygen  Is 

from  air,  and  consequently  specification  of  composition 
and  temperature  requires  that  allowance  be  made  for  the 
presence  of  nitrogen.  l<et  the  two  composition  variables, 

P and  7 be  defined  as: 

P = fraction  of  CaHg  consumed  by  B-l 
y - fraction  of  CO  equivalent  to  C3Hg  consumed  by 
B-2. 

The  composition  of  the  gas  corresponding  to  any  pair  of 
P,7  values  may  now  be  calculated  (let <p  times  the  stolchio- 
metrically  required  fuel  be  used) : 


Component 

mol  pr  mo  CaHg 

C3H8 

1 - p 

0a 

5 (]/<P  ) - 7/2P  - 3/27 

Na 

18.8  (1/<P  ) 

Ha0 

4P 

CO 

3(p-7),  (P<7) 

C03 

37 

Total 

£s  1 + 23.8 (l/<p  )+5/£p  - 3/27 

Temperature 

T = T + -ff-1  p + 7 

0 Cp  Cp 

For  each  reaction  postulated  there  Is  a corresponding  reaction 
rate  expression.  For  the  present  case  assume  that  the 
Initial  cracking  of  propane  Is  the  rate-limiting  step  In 
reaction  B-l  and  that  the  combustion  of  carbon  monoxide  may 
be  represented  by  a simple  second-order  reaction  rate 


The  calculation  of  temperature  for  substitution  into 
B-7  and  B-8  is  simplified  by  the  relation: 

Tr  “ T r + Air  P + *2T  y B-9 

where  Axr  and  &zr  represent  the  corresponding  reduced 
adiabatic  flame  temperature  rises  resulting  from  complete 
reaction  of  reactions  B-l  and  B-2  at  the  prevailing  fuel- 
air  ratio. 

In  calculations  based  on  B-7  and  B-8,  two  cases  may 
be  distinguished;  (l)  where  it  is  desired  to  investigate 
a single  stirred  reactor  with  varying  values  of  % , and  (2) 
where  it  is  desired  to  investigate  staged  reactors  in  each 
member  of  which  the  relative  reaction  rate  parameter  is  the 


same. 


If  a single 


stirred  reactor  is  considered,  then  Gd0  and  Gd7  may  be 
replaced  by  l»  p and  G y respectively.  Corresponding  divi- 
sion of  each  equation  by  0 and  y yields  the  equations: 

-(Ei/RTJ  (1/T  - 15 


E3/RTK(1/T  -1) 
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I 


Each  of  these  equations  may  be  represented  graphically  as 

G vs  0 for  constant  values  of  y (or  vlsa-versa) . B-10 

r 

Is  plotted  for  a £ value  of  1.0.  Superposition  of  the 
two  graphs  permits  solutions  for  the  values  of  0 and  *Y 
corresponding  to  any  given  Gr<  If  these  two  figures  are 
plotted  on  a coordinate  of  log  Qr  then  changes  In  £ are 
represented  by  translation  of  the  graphs  with  respect  to 
each  other  along  the  Gr  axis,  and  for  any  position  desired 
a set  of  0,7  values  may  be  found.  Figures  B-l  and  B-2 
show  such  a graphical  representation  of  equations  B-10 
and  B-ll  for  the  following  constants: 


E!/RTk 

- 5 

Ea/RTK 

* 10 

<P 

* 1.0 

= 0.674 

Aar 

* 0.462 

To 

r 

- 0.25 

tk 

= 2000 °K 

Three  typical  sets  of  solutions  for  £ of  1.0,  0.3  and 
0.1  are  given  In  Figure  B-3.  Since  the  absolute  magnitudes 
of  the  values  of  Gr  in  Figure  B-3  reflect  the  normalization 
procedure  they  are  of  no  significance  In  comparing  these 
curves  with  experimental  stirred  reactor  data.  The  choice 
of  or  relative  reaction  rate,  In  this  two-step  mechanism 
must  be  based  on  the  relative  performances  of  the  stirred 

id 
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reactor  with  respect  to  some  key  point  (such  as  the  rate 
at  the  firing  rate  giving  the  key  temperature) . 'Ihe  correct 
£ having  been  chosen,  the  absolute  values  of  the  reaction 
rate  constants  are  then  found. 

Case  II:  Multi-stage  reactors,  fixed  £ . In  consider- 
ing a multi-stage  reactor  system  It  Is  necessary  to  predict 
the  effect  of  feeding  a stirred  reactor  with  a partially 
burned  feed  when  the  burning  rate  Is  significantly  affected 
by  two  (or  more)  reactions.  So  long  as  both  reactions  are 
significant  In  determining  combustion  rate  the  relative 
rates  of  change  of  0 and  y differ  at  each  point  along  the 
single-stage  stirred  reactor  curve.  Hence  for  each  postu- 
lated degree  of  combustion  of  the  feed  a new  stirred  reactor 
curve  must  be  drawn.  In  the  general  case,  a graphical 
relation  between  Inlet  and  outlet  composition  of  a stirred 
reactor  being  fed  a feed  composition  0Q,  7 is  desired. 

Since  the  rate  of  change  of  0 and  y are  specified 
uniquely  by  0 and  y (i.e.  independent  of  reactor  inlet 
conditions),  one  may  plot  lines  of  constant  d0/d7  on  a 
graph  of  0 vs  7.  On  such  a plot  the  locus  of  all  possible 
reactor  outlet  conditions  for  any  given  0,7  combination 
(or  visa-versa)  may  be  represented.  For  any  combination  of 
0,7  represent  reactor  outlet  conditions  the  (0-0o)/ (7-7Q) 
ratio  across  the  reactor  must  be  equal  to  d0/d7.  From  this, 
loci  of  reactor  outlet  conditions  from  any  given  0Q7O 
combination  may  be  constructed.  Similarly  loci  of  plug-flow 


- 183  - 

reactor  conditions  may  be  found. 

Given  reactor  inlet  and  outlet  conditions  (found  from 
Figure  B-4)  the  feed/volume  ratio  is  easily  found  from 
either  Figures  B-l  or  B-2.  Since  these  figures  were  plotted 
as  fi  (p,*y)/p  and  fa(p,7)/7,  a correction  must  be  made  for 
the  ratio  p/(0-pQ)  or  7/(V7c). 


OF  CONSTANT 


HGURE  B-2:  PERFORMANCE  OF  TWO-STEP  REACTION  IN  A WELL-STIRRED  REACTOR:  STEP 


appendix  c 

ESTIMATED  HEAT  LOSS  PROM  THRKK-IF/p  chamber 


Heat  losses  from  the  experimental  chambers  were  known  to 
be  substantial  because  of  the  high  surface-volume  ratio  and 
because  of  the  opportunity  for  radiant  heat  transmission 
through  and  from  the  Vycor  windows.  Experimental  estimation 
of  this  heat  loss  would  be  very  difficult  because  of  the 
multiple  paths  (radiation,  conduction,  convection)  by  which 
heat  was  lost.  Heat  loss  rates  were  therefore  estimated  by 
calculation.  The  calculations  were  based  on  the  geometry  of 
the  three-jet  combustion  chamber. 

In  calculating  heat  losses  from  this  chamber  several 
parallel  mechanism  must  be  considered: 

(1)  Convection  to  the  firebrick  and  Vycor  walls  followed 

by : 

(a)  direct  or  indirect  radiation  from  the  firebrick 
to  the  ambient  gas 

(b)  transfer  by  convection  and  radiation  from  the 
Vycor  walls  to  the  ambient 

(c)  conduction  through  the  firebrick  from  an  inside 
face  to  an  outside  face  followed  by  radiation- 
conduction  transfer  to  the  ambient. 

(2)  Gas  radiation  from  the  burning  gases  to  the  firebrick 
and  Vycor  walls  followed  by  a,  b,  and/or  c above. 

In  addition  to  these  possible  multiple-mechanisms,  the  problem 
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Is  complicated  by  the  existence  of  strong  temperature  gradients 
in  the  firebrick  and  presumably  also  In  the  Vycor,  especially 
when  the  gases  are  burning  unevenly. 

Because  the  precise  accounting  for  all  of  tnese  losses  Is 
unjustifiably  complex,  simplifying  assumptions  to  the  problem 
were  sought.  These  are  enumerated  as  follows: 

A.  Because  of  the  short  gas  path-lengths  and  low  pressure 
operation,  all  heat  transfer  by  gas  radiation  is  neglected. 

B.  Because  of  the  strong,  temperature  gradients  in  the 
firebrick  simplifications  in  the  radiation  from  and  convection 
to  the  brick  are  necessary.  About  2/3  of  each  side  wall  and 
the  entire  roof  (consisting  of  flues  and  spacer  bricks)  normally 
glowed  visibly  — a bright  cherry  red.  Radiation  from  the  floor 
and  lower  portions  of  the  side-walls  appeared  to  be  much  colder, 
presumably  due  to  colder  gases  and/or  lower  gas  velocities  in 
this  region.  Calculations  were  simplified  by  assuming: 

(1)  uniform  gas  temperature, 

(2)  side-walls  and  roof  of  uniform  temperature,  receiving 
heat  by  convection  from  the  gas  and  radiating  to  the 
floor,  Vycor,  and  ambient, 

(3)  the  roof  has  a heat-transfer  area  equal  to  an  equiva- 
lent plane  (holes  for  the  flues  were  ignored) , 

(4)  radiation  from  the  roof  and  walls  is  black, 

(5)  the  floor  receives  and  looses  heat  only  by  radiation 
and  is  black, 

(6)  heat  loss  by  conduction  through  the  firebrick  is 


- 190 


negligible  because  of  the  low  thermal  conductivity 
of  kaolin  brick. 

C.  In  the  temperature  range  of  Interest,  Vycor  (nearly 
pure  silica)  Is  partly  transparent  and  a partial  reflector. 

The  following  assumptions  were  made: 

(1)  Vycor  Is  transparent  to  radiation  of  wave-length 
shorter  than  3.7p  and  black  to  radiation  of  longer 
wave-length.  Hence  Its  absorptivity  (or  emlsslvlty) 

Is  found  from  the  Stefan-Bolzman  equation  as  a 
function  of  the  temperature  of  the  source  (or  Vycor) . 
Since  all  emmlssors  In  the  chamber  were  postulated 

to  be  black,  complications  arising  from  multiple 
reflect ior.3  are  avoided. 

(2)  Since  the  reflectivity  of  Vycor  Is  low  (10^)  It  Is 
neglected. 

(3)  Thermal  conduction  through  the  Vycor  Is  assumed  to 
be  rapid. 

D.  Heat  transfer  by  convection  has  been  studied  In  various 
geometries,  but  no  data  specifically  adapted  to  this  geometry 

Is  reported  In  the  literature.  Hence  assumed  values  are  neces- 
sary. Since  most  transfer  coefficients  across  gas  films  are 


in  the  range  of  1 to  10  Btu/hr  ft*  °P  (17) , and  are  approxi- 


mately proportional  to  gas  velocity  (16),  the  following  relation 
Is  assumed  for  convection  from  the  burning  gas  to  the  fire- 
brick and  Vycor; 
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hg  - 10  G/Gqq  (Btu/hr  ft3  °P) 

where  G Is  the  firing  rate  and  °B0  Is  the  firing  rate  at 
blowout.  The  transfer  coefficient  from  the  Vycor  to  the 
ambient  air  was  assumed  to  be  hc  * 2 Btu/hr  ft3  °F. 

E.  All  back  radiation  from  the  surroundings  to  the 
chamber  Is  neglected. 

With  these  assumed  conditions  the  heat  loss  from  the 
burning  gas  to  the  ambient  air  may  now  be  calculated.  Radia- 
tive heat  transfer  Is  calculated  by  the  method  of  Hottel  (16). 

For  any  given  gas  temperature  and  firing  rate  all  independent 
parameters  are  fixed  and  the  heat  loss  Is  uniquely  determined. 
However  because  of  the  temperature  dependence  of  the  Vycor 
emisslvity  and  absorptivity  (assumption  C-l  above)  it  is 
necessary  to  evaluate  the  temperatures  of  the  firebrick,  Vycor, 
and  floor.  These  are  found  by  a heat  balance  on  those  respective 
areas. 

A.  Heat  balance  on  walls  and  roof; 

Vw'VV  ■ W <’■ v- < V>  * <V-V) 

C-l 

B.  Heat  balance  on  floor 

<rVwf<V-V)  - a- Vf*  < V- ( -r V>  C-2 

C.  Heat  balance  on  Vycor 

hgVW  + Vf.^vV-^V' 

Vwv<T  ^avwTw  vTv  ^ 

A(r6  vV  + haAv<V  Ta) 


+ 


C-3 
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U 


where : 


‘g 


absorptivity  of  Vycor  for  radiation  as  a function 
of  source  tempera ture.  Subscript  xn  refers  to  wall 
radiation,  x£  refers  to  floor  radiation, 
emlsslvlty  of  the  vycor,  a function  of  Vycor  tempera- 
ture. 

direct  view  factor,  fraction  of  radiation  leaving 
source  arriving  at  sink.  First  subscript  refers 
to  source,  (on  whose  area  the  view  factor  Is  based) 
and  second  subscript  refers  to  sink, 
convection  heat  transfer  coefficient  from  burning 
gas  to  walls  and  Vycor. 

h ■ convection  heat  transfer  coefficient  from  Vycor 
c 

to  ambient  air. 

Simultaneous  solution  of  equations  C-l  to  C-3  for  each  gas 
temperature  and  firing  rate  permits  evaluation  of  T^  and  Ty. 

Heat  loss  from  the  chamber  Is  then  found  by  the  equation: 

c = hgvv  V - Vv(YV  c-4 

Figure  C-l  gives  ay  or  of  Vycor  as  a function  of  Vycor 
or  source  temperature  based  on  the  assumptions  given  in  C above. 
This  plot  must  be  used  in  conjunction  with  a trial  and  error 
solution  of  Equations  C-l  to  C-3.  The  heat  losses  calculated 
from  these  equations  are  shown  as  functions  of  gas  temperature 
and  firing  rate  in  Figure  C-2. 


OR  « OF  VYCOR 
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APPENDIX  D 

DETAILS  OP  EQUIPMENT  AND  APPARATUS 

The  combustion  chamber  used  In  this  work  consisted  of 
three  parallel  two-dimensional  jets  mounted  In  a framework 
holding  a firebrick  lining.  The  front  and  back  of  the  frame 
were  closed  with  Vycor  windows  to  permit  ease  of  observation. 

In  this  chamber  recirculation  was  created  by  the  aspirating 
effect  of  the  Jets.  This  recirculation  could  be  controlled, 

If  desired,  by  varying  the  spacing  of  the  Jets;  however, 

In  these  experiments  only  one  Jet  spacing  was  used.  Sketches 
showing  the  dimensions  of  the  chamber  and  nozzles  are  given  in 
Figures  D-l  and  D-2.  (An  Isometric  sketch  of  the  assembly 
has  been  given  in  the  main  body  of  this  thesis  Figure  IV-i). 

The  feed  to  this  chamber  was  a pro-mixed  fuel-air  mixture. 
Fuel  and  air  were  metered  separately  and  mixed.  This  mixture 
was  then  manifolded  and  the  manifold  mixture  to  each  jet  was 
metered  to  Insure  equal  distribution  of  feed.  The  orifices 
used  in  the  primary  air,  fuel  and  manifold  mixture  lines  were 
standard  square-edge  orifices  and  were  calibrated  to  Insure 
accuracy  of  metering.  Calibration  was  done  with  air  against 
a standard  orifice  set.  Corrections  were  made  to  the  cali- 
bration curve  for  fuel  to  account  for  the  compressibility  and 
higher  molecular  weight  of  propane.  A similar  correction 
should  have  been  made  to  the  calibration  curves  for  the  manifold 


mixture  orifices  (addition  of  propane  changed  the  molecular 
weight  but  did  not  affect  compressibility);  however  this 
correction  was  small,  and  since  it  entered  into  all  three 
equally,  it  was  not  necessary  to  achieve  the  desired  distri- 
bution of  feed.  A flow  diagram  showing  the  orifice  designa- 
tions is  given  in  Figure  E-3.  Calibration  curves  for  each 
of  the  orifices  are  given  in  Figures  D-4  to  D-9.  Compressibil- 
ity factor  for  propane  is  given  in  Figure  E-10.  (Calculated 
from  data  reported  by  Stearns  and  Oeorge  (25) ) . 

A sampling  train  was  constructed  to  facilitate  sampling 
of  the  chamber  gases.  (A  flow  diagram  of  this  train  is  given 
in  Figure  IV-12) . In  sampling  gases  where  the  gases  impinge 
directly  on  the  probe  opening,  the  velocity  of  the  gases  at 
the  probe  should  be  matched  to  the  sampling  velocity.  To  do 
this  it  1 8 necessary  to  meter  the  sampling  rate.  This  metering 
was  accomplished  by  means  of  a capillary  flow  meter.  The 
metering  element  was  a short  piece  (about  4")  of  hypodermic 
tubing  of  l.d.  0.049".  This  was  calibrated  against  a wet-test 
meter.  The  calibration  curve  is  given  in  Figure  D-ll. 


5-40  MACH.  SCREW 
SILVER  SOLDERED 
INTO  NOZZLE  HEAD, 
TWO  EACti  SIDE. 


FIGURE  0—1:  ENLARGED  DETAIL  OF  NOZZLE.  SCALE:  2”*  I 


QUENCH  AIR 


FIGURE  0-5:  CALIBRATION  CURVE,  ORIFICE  0-1  MED.  RANGE 
ORIFICE:  0.500"  IN  ly  " PIPE 
FLUID:  AIR 


FIGURE  D-6:  CALIBRATJJN  CURVE,  ORIFICE  O-l  HIGH  RANGE 
ORIFICE.  0.706*  IN  ly"  PIPE 
FLUID:  AIR 


-•7s  ORIFICE  CALIBRATION  i 
ORIFICE:  0.082"  IN  2"  PIPE 
FLUID:  PROPANE 


FIGURE  0-8 : CALIBRATION  CURVE,  ORIFICE  0*2  HIGH  RANGE 

ORIFICE:  a 1 163"  IN  2*  PIPE 
FLUID:  PROPANE 


FIGURE  D-9:  CALIBRATION  CURVES:  ORIFICES  0-4  TO 

ORIFICES:  0.4375“  IN  $■"  PIPE 
FLUID  AIR 


PRESSURE,  PSIA 

FIGURE  0-10:  COMPRESSIBILITY  OF  PROPANE 
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FIGURE  0- 1 1:  CALIBRATION  CURVE  FOR  CAPILLARY  FLOW 
METER  IN  SAMPLING  TRAIN.  FLUIO-AIR. 
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APPENDIX  E 
SAMPLE  CALCULATIONS 

A.  Calculation  of  blowout  or  loading  rates  (example,  Run  4B) 

TABLE  E-l 

DATA  OBTAINED  IN  A TYPICAL  RUN 


Station 

P 

AP 

T 

0-1 

24.2  psig 

9,22"  H20 

525°R 

(0.706"  orifice) 

0-2 

34.8 

30.35 

513 

(0.1163"  orifice) 

0-4 

-4.21"Hg* 

17.14 

527 

0-5 

-4.23”* 

17.00 

527 

0-6 

-4.26"# 

16,15 

527 

P-7 

19.40"* 

atm 

30.06" 

* Pressure  measured  in  inches  of  mercury  vacuum. 

Negative  pressure  indicates  pressure  above  atmospheric. 

From  these  data  flow  rates,  w,  in  lb/3ec  of  gas  are  calculated. 
Since  the  calibration  charts  for  the  various  orifices  are  given 
in  (PdP)/T  units  of  (lb/ft2)2/°R,  a conversion  factor  must  be 
used  to  account  for  the  units  of  the  measurements.  The  flow 
rates  calculated  from  the  data  of  Table  E-l  are  given  in 
Table  E-2.  It  should  be  noted  from  the  flow  diagram,  Figure 
D-3,  that  the  sum  of  the  flow  rates  through  orifices  0-1  and 
0-2  should  equal  the  flow  through  orifices  0-4  to  0-6.  Loading 
rates,  fuel-air  ratio  and  other  calculated  quantities  are  also 
indicated  in  Table  E-2. 


CALCULATED  QUANTITIES 


Loading  rate  » 0.761 

Temperature  correction  ■ 1/0.973  (Figure  E-2) 
Loading  rate  corrected  to  540°R  * 0.7^2 


CALCULATED  QUANTITIES 
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Temperature  correction  * 1/0.973  (Figure  E-2) 
Loading  rate  corrected  to  540°R  ■ 0.782 


CORRECTION  FACTOR 


0 2 4 6 8 10  12  13 


AT  FROM  540*R 

FIOURE  E*  2;  TEMPERATURE  CORRECTION  FACTORS  (l£).  FOR  TEMPERATURE 
ABOVE  340*R  MULTIPLY  BY  FACTOR.  FOR  TEMPERATURE  BELOW  S40*R 
DIVIDE  BY  FACTOR. 
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B.  Calculation  of  combustion  efficiency  (example  Run  86A) 


TABLE  E-3 


DATA  OBTAINED 

IN  A TYPICAL 

RUN 

a) 

Loading  rate 

P 

AP 

T 

0-1 

56.3  psig 

22.26"  h2o 

537°R 

(0.348"  orifice) 

0-2 

29.7  psig 

14.22"" 

531° 

(0.1163"  orifice) 

0-5 

537° 

0-7 

20.04"  Hg  vac 

- 

0-9 

30.23"  Hg 

i b) 

sampling  data  — 

probe  position,  cold 

side  of  flue 

1 

rate  data 

P 

AP 

T 

P°2 

Impact  10.19"  Hg  4.9"  H20 

Sample  20.10  mm  Hg  33.1  oil  297°K  32.3  ram  Hg 

From  the  above  data  the  following  quantities  are  calculated: 

a)  Loading  rate  (calculated  as  in  previous  example); 

N./VP1*^)  s 0.620,  F = 0.575,  A/F  a 11.53  Ib/lb 
* 540 

b)  Sampling  rate; 

PAP/T  = (20.10) (33. l)/(297)  • • 0.044 
w - 3.67  x 10~5  lb/sec  (Figure  E-3) 

c)  Composition; 

0.7  mm 


Ng  correction*  • 

Mol  dry  products. 
Mol  wet  products 


0.856 


(Figure  E-3) 
(Figure  E-4) 


* A correction  to  magnetic  oxygen  meters  must  be  jaade  to  correct 
for  the  magnetic  properties  of  other  gases  present.  It  is  assumed 
for  this  work  that  the  background  gas  is  only  N2.  Only  a small 
error  is  introduced  by  this  assumption. 
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c)  Composition  — continued 

Moles  O2/  100  mols  wet  products  B 

(33*5) (.856) (100)/(20.10)  s 14.26 

02  req' d/100  mol  wet  prod.  ■ 17.74  (Figure  E-5) 

s 1 - (02  remaining )/(02  req*d)  » 19.7# 

d)  Estimated  gas  temperature; 

Adiabatic  flame  temperature  * 3755°R  (Figure  E-6) 

(Taf  “ 530)  + 530  * 1145°R 

e)  Sample  rate  required  by  impact  measurement; 

PAP/T  s (10.19) (4.9)/(1145)  - 0.044  (Figure  E-7) 

w = 5.79  x 10-5  (Figure  E-7) 

w required  ■ (.656) (5.79  x 10”5)  a 4.9  x 10“5 

(to  be  compared  with  3.67  x 10“5  used) 

C.  Pressure  eixponent  calculation  (example  Run  76A) 


TABLE  E-5 


DATA  OBTAINED 

IN  A TYPICAL 

RUN 

P 

A? 

T 

0-1 

57.0  psig 

2.02"  H20 

544°R 

(0.34S"  orifice) 

0-2 

43.3  psig 

6.20" 

536° 

(0.1163”  orifice) 

0-5 

5400 

P-7 

23.90”  Hg  vac 

atm 

30. OS"  Hg 

From  these  data  the  following  quantites  are  calculated  as  in 
example  A above: 

w of  air  * 0.0129  lb/sec 

Fuel  fraction  F » 0.599 


C.  Pressure  exponent  calculation  — continued 


Vv 

= w.  x 4.26 

A 

* 0.0549 

Temp 

= 540°R 

Temp.  Correction  * 1.00 

F 

- .599 

F correction  c 0.993 

F correction  calculated  from  relative  blowout  rates 
of  actual  fuel-air  mixture  and  the  fuel-air  mixture  desired. 
On  the  rich  side  this  is  F corr  * e 7^r'^  and  on 

the  lean  aide  this  is  F corr  * e*tf  . 

VV)540.  .600  " °-°5« 

P « 6.18"  Hg  • 0.207  atm 


A series  of  these  points  are  taken  for  a given  fuel/air 
ratio  and  plotted  on  log-log  paper.  The  slope  of  the  resulting 
line  is  determined. 


FIGURE  E-4:  MOLES  OF  DRY  COMBUST  ION  PRODUCTS 


FIGURE  E-5:  MOLES  OF  OXYGEN  CONSUMED  AND  REMAINING  AT  COMPLETE  COMBUSTION  (£&) 


•I* 


FIGURE  E-6:  ADIABATIC  FLAME  TEMPERATURE,  To  = 70° F '.?.§) 


- 220 


APPENDIX  F 


SUMMARY  OF  CALCULATED  VALUES 


TABLE- F-l 

BLOWOUT  MEASUREMENTS,  THREE- JET  CHAMBER 


Run 

Air  rat® 

Pres. 

Temp 

Loading,  Nj/VP1*# 

Gen*l 

Zones 

lb/sec 

in  Hg 

°R 

lb  mol/s.ec  ft^  atm^« 
raw  corr'd 

3 fuel 
fr'n 

blown 

out 

4B 

0.0424 

10.61 

527 

0.761 

0.732 

0.524 

1 

4C 

0.0431 

10.62 

527 

0.722 

0.791 

0.493 

1 

4D 

0.0445 

IO.64 

527 

0.792 

0.315 

0.445 

2 

4E 

0.0420 

10.66 

525 

0.741 

0.765 

0.525 

1 

4F 

0.0420 

10.66 

525 

0.741 

0.765 

0.529 

2 

40 

0.0420 

10.65 

522 

0.741 

0.770 

0.469 

1 

4H 

0.0417 

10.64 

525 

0.763 

0.790 

O.442 

2 

41 

0.0445 

10.71 

524 

0.803 

0.332 

0.483 

1 

4J 

0.0473 

10.69 

536 

0.357 

0.393 

0.470 

2 

5A 

0.027$ 

10.46 

0.525 

0.530 

0.558 

1 

5B 

0.0273 

10.43 

537 

0.524 

0.529 

0.569 

2 

6A 

0.0264 

10.06 

533 

0.533 

0.536 

0.554 

1 

6B 

0.0263 

9.95 

540 

0.542 

0.542 

0.607 

2 

6c 

0.0272 

10.14 

542 

0.541 

0.537 

0.457 

1 

6D 

0.0272 

10.07 

542 

0.547 

0.545 

0.428 

2 

6e 

0.0263 

10.14 

540 

0.534 

0.534 

0.552 

1 

6F 

0.0263 

10.11 

541 

0.534 

0.532 

0.602 

2 

60 

0.0270 

10.10 

542 

0.539 

0.537 

0.460 

1 

6H 

0.0270 

10.07 

542 

0.545 

0.542 

0.424 

2 

6T 

0.01365 

10.23 

542 

0.361 

0.359 

0.570 

1 

6J 

0.01365 

10.27 

542 

0.361 

0.360 

0.019 

2 

6K 

0.01380 

10.13 

542 

0.372 

0.370 

0.433 

1 

6L 

0.01330 

10.13 

542 

0.372 

0.370 

0.415 

2 

6m 

0.01900 

10.26 

542 

0.371 

0.369 

0.506 

1 

6N 

0.01900 

10.21 

542 

0.372 

0.370 

0.616 

2 

60 

0.01900 

10.23 

542 

0.372 

0.370 

0.447 

* 

6P 

0.1900 

10.21 

542 

0.376 

0.374 

0.413 

2 

7* 

0.0222 

10.31 

533 

0.423 

0.430 

0.557 

1 

7B 

0.0221 

10.29 

540 

0.423 

0.423 

0.620 

2 

7C 

0.0220 

10.32 

540 

0.424 

0.424 

O.485 

1 
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TABLE  F-l  continued 


Run 

Air  rate 

Pres. 

Temp 

Loading, 

NA/Vpl-S 

Genfl 

Zones 

lb/sec 

in  Hg 

°R 

lb  mol/sec 
raw 

ft^  atm^-* 
corr’d 

8 fuel 
fr’n 

blown 

out 

7D 

0.0221 

10.29 

540 

0.428 

0.428 

0.417 

2 

7F 

0.0219 

10.35 

541 

0.420 

0.419 

0.560 

1 

7G 

0.0216 

10.31 

541 

0.416 

0.415 

0.607 

2 

7H 

0.0214 

10.36 

540 

0.411 

0.411 

0.442 

1 

71 

0.0215 

10.30 

540 

0.417 

0.417 

0.411 

2 

7J 

0.0314 

10.31 

540 

0.605 

0.605 

0.535 

1 

7K 

0.030S 

10.19 

(540) 

0.606 

0.606 

0.596 

2 

7L 

0.0323 

10.33 

541 

0.623 

0.622 

0.528 

1 

7M 

O.O323 

10.24 

543 

0.633 

0.629 

0.590 

2 

7N 

0.032S 

10.44 

540 

0.620 

0.620 

0.470 

1 

70 

0.032S 

10.36 

542 

0.629 

0.626 

0.427 

2 

7P 

0.0329 

10.44 

538 

0.623 

0.626 

0.471 

1 

7Q 

0.032S 

10.42 

538 

0.623 

0.626 

0.427 

2 

*32A 

0.0411 

10.26 

526 

0.803 

0.828 

0.495 

1 

*32B 

0.0450 

9.78 

523 

0.889 

0.919 

0.456 

2 

*3  4 A 

0.C279 

10.56 

533 

0.511 

0.523 

0.626 

2 

*34B 

O.C315 

10.60 

530 

0.5C0 

0.600 

0.633 

2 

*34C 

0.0276 

10.63 

533 

0.504 

0.516 

0.648 

2 

*34D 

0.0246 

10.69 

534 

0.445 

0.455 

0.657 

2 

*35B 

0.0296 

10.49 

532 

0.490 

0.502 

0.620 

1 

*35C 

0.01735 

10.34 

532 

0.332 

0.341 

0.655 

1 

36A 

0.0291 

3.0.34 

534 

0.55S 

0.568 

0.410 

1 

36b 

0.0367 

10.36 

534 

0.738 

0.750 

0.414 

1 

36c 

0.0487 

10.42 

531 

0.918 

0.937 

0.431 

1 

37A 

0.0413 

9.38 

530 

0.944 

0.965 

0.534 

2 

37B 

0.0440 

9.37 

528 

1.009 

1.034 

0.532 

2 

37C 

0.0418 

9.42 

530 

0.951 

0.970 

0.531 

2 

37D 

0.0430 

9.40 

528 

0.982 

1.000 

0.507 

2 

37E 

0.0446 

9.48 

529 

1.001 

1.025 

0.532 

2 

37F 

0.0424 

9.47 

529 

0.9  44 

0.966 

0.533 

2 

*38 

0.0349 

9.13 

533 

0.839 

0.854 

1.018 

0.557 

2 

*39A 

0.0403 

8.93 

534 

1.003 

0.561 

2 

*39B 

0.0402 

8.95 

532 

1.001 

1.023 

0.571 

2 

*39C 

0.0396 

8.SS 

531 

1.000 

1.023 

0.437 

2 

*39D 

0.0395 

8.89 

531 

0.996 

1.017 

0.694 

2 

*39E 

0.0386 

8.96 

531 

0.956 

0.980 

0.592 

2 

*39F 

0.0386 

9.05 

530 

0.940 

0.967 

0.598 

2 

*39G 

0.0385 

8.95 

530 

0.959 

0.986 

0.434 

2 

*391 

0.0470 

8.96 

530 

1.177 

1.202 

0.469 

2 

*39J 

0.0468 

8.92 

528 

1.173 

1.222 

O.484 

2 

* These  data  were  disregarded  in  plotting  as  they  were  unreasonable 
compared  to  the  rest  of  the  data.  These  values  reflect  probably 

the  consequences  of  air  leakage  and  poorly  fitting  firebricks. 


- 222  - 


TABLE  F-l  continued 


Run 

Air  rate 

Pres. 

Temp 

Loading, 

Na/VP1^ 

Gen’l 

Zone3 

lb/sec 

in  Hg 

OR 

lb  mol/sec  f*3  atm^>3 
raw  corr ’ d 

fuel 

frfn 

blovm 

out 

40A 

0.0367 

10,47 

533 

0.690 

0.700 

0.609 

1 

40B 

0.03 6 A 

10.43 

533 

0.681 

0.690 

0.498 

1 

40C 

0.0364 

10.34 

531 

0.638 

0.703 

0.434 

2 

40D 

0.0309 

10.50 

531 

0.579 

0.591 

0.605 

1 

40E 

0,0312 

10.60 

531 

0.574 

0.585 

0,467 

1 

40F 

0,0312 

10.50 

531 

0.578 

0.590 

0.447 

2 

41A 

0.0330 

IO.65 

533 

0.606 

0.617 

O.56I 

1 

41B 

0.0329 

10.69 

532 

0.593 

0.603 

0.474 

1 

41C 

0.0329 

10.61 

531 

0.429 

0.593 

0.607 

2 

41D 

0.0335 

10.75 

531 

0.598 

0.612 

0.559 

1 

41E 

0.0330 

10.75 

531 

0.590 

0.601 

0.467 

1 

41F 

0.0330 

10.70 

531 

0.592 

0.605 

0.442 

2 

41G 

0.0270 

10.43 

530 

0.476 

0.487 

0.582 

1 

41H 

0.0270 

10.37 

530 

0.469 

0.485 

0.432 

2 

411 

0.026 8 

10.80 

530 

0.474 

0.485 

0.465 

1 

41J 

0.0268 

10.73 

530 

0.476 

0.488 

0.428 

2 

42A 

0.0327 

9.23 

532 

0.763 

0.730 

0.437 

1 

42B 

0.0336 

9.23 

532 

0.787 

0.798 

0.506 

1 

42C 

0.0336 

9.24 

534 

0.790 

0.800 

0.598 

1 

*43  A 

0.0317 

10.13 

529 

0.631 

O.648 

0.454 

2 

*43  B 

0.0297 

9.38 

530 

0.560 

0.579 

0.610 

2 

*43  C 

O.O463 

10.00 

530 

0.935 

0.952 

0.469 

2 

*43  D 

0.0304 

9.38 

528 

0.632 

0.655 

0.629 

2 

44A 

0.0377 

10.84 

525 

0.633 

0.691 

0.590 

2 

44B 

0.0602 

10.70 

524 

I.O85 

1.120 

0.475 

2 

*45A 

0.0335 

10.30 

524 

0.742 

0.735 

0.629 

2 

*45B 

0.0480 

9.45 

522 

1.090 

1.150 

0.603 

2 

*45C 

0.0432 

9.48 

523 

1.070 

1.125 

0.593 

2 

*45D 

0.0433 

9.52 

523 

0.972 

1.020 

0.612 

2 

*48A 

0.0243 

10.17 

529 

0.498 

0.515 

0.590 

1 

*48B 

0.0243 

10.19 

528 

0.498 

0.517 

0.629 

2 

*480 

0.0417 

8*46 

525 

1.166 

1.203 

0.540 

1 

*43D 

0.0410 

8.46 

528 

1.100 

1.147 

0.533 

2 

49A 

0.0255 

10.25 

540 

0.498 

0.493 

0.561 

1 

49B 

0.0253 

10.30 

540 

0.490 

0.490 

0.597 

2 

49C 

0.0289 

10.27 

537 

0.560 

0.561 

0.459 

1 

49D 

0.0302 

10.27 

536 

0.587 

C.592 

0.448 

2 

49E 

0.0311 

10.35 

537 

0.597 

0.603 

0.555 

1 

49F 

0.0311 

10.27 

536 

0.605 

0.610 

0.597 

2 

49G 

0.0332 

10.39 

538 

0.633 

O.636 

0.462 

1 

* These  data  were  disregarded  in  plotting  as  they  were  unreasonable 
compared  to  the  rest  of  the  data.  These  values  reflect  probably 
the  consequences  of  air  leakage  and  poorly  fitting  firebricks. 


- 223 


TABLE  P-1  concluded 


Run 

Air  rate 

Pres. 

Temp 

Loading, 

N./VP1**5 

Gen*l 

Zones 

lb/sec 

in  Hg 

°R 

lb  mol/sec 
raw 

ft3  atm1*3 
corr’d 

fuel 

fr’n 

blown 

out 

49H 

0.0363 

10.34 

536 

0.704 

0.711 

0.434 

2 

491 

0.0347 

10.39 

535 

0.662 

0.670 

0.540 

1 

49J 

0.0347 

10. 46 

536 

0.650 

0.657 

0.596 

2 

*51 

0.023d 

10.16 

532 

0.466 

0.461 

0.532 

2 

52A 

0.0317 

11.12 

535 

0.527 

0.534 

O.634 

2 

52B 

0.0376 

12>.22 

535 

0.626 

0.635 

0.622 

2 

52C 

0.0350 

11.20 

533 

0.560 

0.593 

0^540 

2 

52D 

0.0367 

11.21 

533 

0.606 

0.617 

0.535 

2 

52E 

0.0590 

11.2$ 

529 

0.973 

0.995 

0.520 

2 

52P 

0.0311 

11.22 

529 

0.515 

0.531 

0.593 

2 

*53  A 

0.634 

O.656 

O.636 

*53A 

0.0340 

10.49 

531 

2 

*535 

0.0402 

10.61 

533 

0.737 

0.752 

0.637 

2 

*53  C 

O.O646 

10.99 

534 

1.205 

1.22 

0.504 

2 

*53D 

0.0402 

9.46 

535 

0.902 

0.913 

0.604 

2 

*53E 

O.O56O 

9.52 

537 

1.245 

1.25 

0.545 

2 

*54F 

0.0500 

9%62 

537 

1.099 

1.095 

0.569 

2 

* These  data  were  disregarded  in  plotting  as  they  v;ere  unreasonable 
compared  to  the  rest  of  the  data.  These  values  reflect  probably 
the  consequences  of  air  leakage  and  poorly  fitting  firbricks. 
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TABLE  F-2 
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TABLE  Fi-3 

COMBUSTION  EFFICIENCY,  THREE- JET  CHAMBER 


w • 0.0239  - 0.0005  #/sec,  P«  10.5  - 0.4"Hg,  T * 531°  z 2°R 
Loading  rate,  * 0.467  lb  taOl/sec  ft^  atm^*^ 


Run 

Gen»l 

fuel 

frfn 

Left 

flue 

Cntr 

flue 

Right 

flue 

26 

0.549 

66.7# 

62. 1* 

70.0# 

29A 

0.566 

72.1 

64.3 

76.6 

29B 

0.561 

69.6 

60. S 

72.9 

29C 

0i547 

72.6 

60.4 

76.7 

29D 

0.534 

72.6 

62.6 

77.6 
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TABLE  F-A 


BLOWOUT  MEASUREMENTS,  ONE- JET  CHAMBER 


Run 

Air  rate 

Pres. 

Tgmp 

Loading, 

N /ypl.8 

Gen*l 

lb/sec 

in  Hg 

lb  mol/sec 
raw 

ft3  atm**** 
corr*d 

fuel 

fr»n 

60A 

0.00934 

10.18 

530 

0.276 

0.287 

0.661 

60B 

0.00929 

10.19 

530 

0.274 

0.285 

0.661 

60C 

0.00923 

10.18 

530 

0.276 

0.286 

0.651 

60D 

0.00924 

10.17 

531 

0.274 

0.266 

0.286 

0.316 

60E 

0.00907 

10.25 

531 

0.277 

0.316 

60F 

0.01180 

10.26 

531 

0.346 

0.356 

0.630 

600 

0.01175 

10.20 

532 

0.347 

0.356 

0.633 

60H 

0.01175 

10.28 

532 

0.343 

0.353 

0.371 

601 

0.0118 

10.24 

532 

0.347 

0.357 

0.359 

60J 

0.0138 

10.25 

532 

0.405 

0.415 

0.624 

60K 

0.0139 

10.23 

532 

0.410 

0.420 

0.625 

60L 

0.0138 

10.27 

532 

0.403 

0.415 

0.359 

60M 

0.013d 

10.39 

531 

0.395 

0.409 

0.356 

60N 

0.01685 

10.27 

531 

0.493 

0.507 

0.618 

600 

0.01680 

10.34 

532 

0.485 

0.496 

0.608 

60P 

0.01685 

10.27 

532 

0.493 

0.505 

0.617 

60Q 

0.01685 

10.30 

532 

0.491 

0.503 

0.371 

60R 

0.01685 

10.20 

532 

0.499 

0.531 

0.371 

61A 

0.01505 

9.30 

535 

0.524 

0.530 

0.625 

61B 

0.01495 

9.19 

536 

0.533 

0.540 

0.368 

61C 

0.01840 

9.44 

536 

0.626 

0.633 

0.607 

610 

0.01840 

9.32 

538 

0.639 

0.641 

0.590 

61E 

0.01835 

9.33 

53$ 

0.637 

0.640 

0.606 

61F 

0.01840 

9.37 

538 

0.633 

0.636 

0.392 

61G 

0.01840 

9.37 

538 

0.634 

0.637 

0.383 

61H 

0.0225 

9.1$ 

538 

0.805 

0.809 

0.581 

611 

0.0226 

9.18 

538 

0.808 

0.809 

0.585 

6lJ 

0.0226 

9.21 

53$ 

0.802 

0.807 

0.410 

61K 

0.0225 

9.26 

53$ 

0.794 

0.799 

0.405 
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**  At  this  high  firing  rate  the  flame  is  not  stable  when  the  probfe  is  inserted 
across  the  interior  of  the  chamber  this  distance. 


PRESSURE  EXPONENT  DATA,  ONE-JET  C: 
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APPENDIX  G 

SCHLIEREN  PHOTOGRAPHY  OF  THE  THREE-JET  CHAMBER 

In  developing  conceptual  models  of  combustion  chambers 
knowledge  of  the  flow  pattern  as  well  as  of  the  combustion 
efficiency  pattern  Is  useful  In  defining  regions  which  may  be 
considered  as  one  zone.  Since  Schlieren  photography  permits 
visualization  of  density  gradients  In  a gas  it  was  hoped  that 
either  still  or  motion  pictures  taken  via  a Schlieren  apparatus 
might  help  in  defining  appropriate  zones.  Still  pictures  would 
help  show  hot  and  cold  regions  and  motion  pictures  would  help 
In  visualizing  flow  patterns. 

The  principles  of  Schlieren  photography  have  been  widely 
discussed  in  the  literature.  A summary  has  been  given  by 
Shipman  (20)  in  discussing  the  design  and  operation  of  a small 
Schlieren  apparatus  for  the  ft.  I.  T.  Fuels  Research  -Laboratory. 

This  unit  has  two  6"  diameter,  42"  focal-length  mirrors  as  the 
principal  optics.  These  are  mounted  on  a heavy  steel  frame  to 
maintain  their  alignment.  The  frame  was  designed  to  permit 
assembly  around  moderate  size  apparatus.  The  vacuum  tank  used 
in  this  work  was  fitted  with  two  opposing  windows. 

Photographs  were  taken  using  a xenon  discharge  tube 
supplied  through  the  courtesy  of  Dr.  H.  E.  Edgerton  (7).  This 
tube  gave  a flash  duration  of  2 to  3 msec  of  sufficient  intensity 
to  easily  expose  a panchromatic  film  with  an  ASA  rating  of 
100. 
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The  presence  of  the  turbulence  found  with  these  short- 
duration-flash  pictures  suggested  that  the  regions  which 
were  considered  as  well-stirred  zones  were  In  fact  somewhat 
Inhomogeneous.  The  presence  of  this  inhomogeniety  also 
suggested  that  a single  parcel  of  gas  could  be  Identified 
adequately  to  permit  one  to  follow  Its  movement  over  a short 
timed  Interval.  It  was  proposed  to  follow  this  movement  by 
means  of  a double-flash  picture.  This  was  achieved  with  a 
special  flash  control  circuit  supplied  by  hr.  tdgerton.  The 
flash  was  a spark  discharge  of  about  0.1  sec  duration  with 
intervals  ranging  between  flashes  of  5 to  100  secs.  The 
intensity  of  this  flash  was  low  enough  that  special  developer 
was  required  on  high-speed  panchromatic  film  to  achieve  an 
effective  film  speed  of  600  to  800  ASA.  These  pictures  did 
not  show  the  turbulence  with  sufficient  definition  to  follow 
a single  gas  parcel  ; however,  the  shorter  flash  durations 


[. 


[ 


did  give  substantially  clearer  resolution  of  the  Schlieren 
patterns. 

Finally  high-speed  motion  pictures  (2-3000  frames  per  sec) 
were  taken  of  the  flame  using  an  Eastman  High  Speed  camera 
equipped  with  a synchronizing  device  to  time  a strobe-light 
discharge  with  each  frame  exposed,  again  supplied  through  the 
courtesy  of  br.  Edgerton.  These  pictures  showed  a substantial 
gas  movement;  however  the  movement  suggested  that  it  was  hot 
gases  external  to  the  combustion  zone  which  were  being  observed 
rather  than  flow  pattern  within  the  chamber  Itself.  Typical 


FIGURE  G-l:  TYPICAL  SCHLIEREN  PHOTOGRAPHS.  LEFT,  TWO  ZONES 
ACTIVE;  RIGHT,  ONE  ZONE  ACTIVE.  FLASH  DURATION  O.l  pSEC. 
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appeml»ix  I 


TABLE  GP  NOMENCLATURE 


a Proportionality  constant,  Equation  III-8. 

C Concentration,  moles  per  unit  volume.  C^,  fuel 
concentration  cq3»  oxygen  concentration. 

E Effective  activation  energy  of  overall  reaction, 
stirred  reactor  equations,  Cti,  II  and  III. 

F Generalized  fuel  fraction,  F = 0/(1  + 0)  . Appendix  0; 

F is  direct  view  factor  in  radiation  calculations. 

G Firing  rate  to  stirred  reactor  or  combustion  chamber, 
moles  of  fuel  per  unit  time. 

h Heat  transfer  coefficient,  Appendix  L;  h , gas  film 
coefficient  from  burning  gas  to  Vycor  and  walls; 

hc  convection  coefficient  from  Vycor  walls  to  ambient  air. 

J Fractional  contribution  of  a reactor  unit  in  combustor 
model  to  overall  combustion  rate. 

k Proportionality  constants,  ki  and  k2  in  equations 
III-l  and  2 respectively. 

k1  Pseudo  collision  factor  in  overall  reaction.  Well-stirred 
reactor  equations,  Ch.  II  and  III. 

L Length  of  combustor,  Equation  III-4,  ff. 

m Momentum  flux,  Equation  III-ll,  ff. 

n Overall  order  of  reaction,  well-stirred  reactor  relations. 

N.  Air  firing  rate  to  well-stirred  reactor  and  experimental 
chamber,  lb  moles/sec. 

p Number  of  zones  in  extended  recirculation  model. 

P Pressure,  atmospheres. 

i Gas  constant,  aeclrculation  ratio  in  recirculation  models. 

T Temperature,  degrees  ranklne.  T f,  temperature  of  an 
adiabatic  flame  at  complete  1 combustion,  T^, 


: 


nozzle  temperature,  T , inlet  temperature  of  fuel/air 
mixture,  TK,  key  ° temperature  in  normalization 
of  well-  stirred  reactor  equation,  Appendix  A. 

U Velocity,  Chapter  III,  U.  Inward  velocity  (toward  Jet 
axis),  U , outward  velocity  (away  from  Jet  axis), 

nozzle  velocity  of  jet.  Chapter  VI;  overall 
heat  transfer  coefficient,  defined  by  Equation  VI-1. 

V Volume,  cubic  feet. 

x Distance  coordinate  parallel  to  Jet  axis. 

y Distance  coordinate  perpendicular  to  Jet  axis. 

Y Fraction  of  fresh  fuel/air  mixture  reaching  a combustion 
zone. 

2 Fraction  of  fuel  fed  to  recirculation  zone  in  split-fuel 
feed  model. 

a Reaction  order  exponent  on  fuel,  well-stirrea  reactor 
equations . 

0 tfurnedness,  fractional  fuel  consumption. 

y Variables  describing  the  composition  of  a burning  mixture. 

8 Variables  describing  the  composition  of  a Duming  mixture. 

6 Emissivity  m radiation  calculations,  Appendix  1>. 

^ Oxygen  consumption  efficiency. 

£ Relative  reaction  rates,  defined  by  equation  D-  , App.  D. 

P Density. 

T*  Stefan-Bolzman  Constant,  0.171  x 10 . 

^ Summation  of  moles  present  per  mole  of  fuel. 


0 Fuel/air  ratio;  fraction  of  stoichiometric  fuel. 
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